Prace | 
Towarzystwa Naukowego Warszawskiego. 


ИЕ, — Wydział nauk matematycznych | przyrodniczych. 


Travaux de la Société des Sciences de Varsovie. 
HL Classe des sciences mathématiques et naturelles. 


N 12. 
JAN KRASSOWSKI. 


0 RUCHU PLANETOID TYPU 3/4 (THULE). 


WYZNACZENIE WYRAZÓW TYPOWYCH ORAZ PRZERWY. 
SUR LE MOUVEMENT DES PETITES PLANÈTES 

DU TYPE DE THULE (8/4). ў 

N 1144 27 Na 


DÉTERMINATION DES TERMES ÉLÉMENTAIRES ET DE LA LACUNE 
POUR LE TYPE THULE. 


WYDANE 7 ZAPOMOGI KASY POMOCY DLA OSÓB PRACUJACYCH 
NA POLU NAUKOWEM IMIENIA D-ra JÓZEFA MIANOWSKIEGO. 


WARSZAWA. 
NAKŁADEM TOWARZYSTWA NAUKOWEGO WARSZAWSKIEGO 
Skład główny w księgarni E. Wende i S-ki. 


1916. 
Dodatek do „Sprawozdań z posiedzeń* T. N. W. 


À 


| { 
http://rcin.org.pl 


Trace И м, 
Towarzystwa Naukowego Warszawskiego. 


Ill. — Wydział nauk matematycznych i przyrodniczych. A 


Travaux de la Société des Sciences de Varsovie. 
HI. Classe des sciences mathématiques et naturelles. 


№ 12. 
JAN KRASSOWSKI. 


0 RUCHU PLANETOID TYPU *4 (THULE), 


WYZNACZENIE WYRAZÓW TYPOWYCH ORAZ PRZERWY. 


SUR LE MOUVEMENT DES PETITES PLANÈTES 
DU TYPE DE THULE (3/4). 


DÉTERMINATION DES TERMES ÉLÉMENTAIRES ET DE LA LACUNE 
POUR LE TYPE THULE. 


ЫИ Pafstu. ЫРШ re 


WYDANE Z ZAPOMOGI KASY POMOCY DLA OSÓB PRACUJACYCH 
NA POLU NAUKOWEM IMIENIA D-ra JÓZEFA MIANOWSKIEGO. 


WARSZAWA 
NAKLADEM TOWARZYSTWA NAUKOWEGO WARSZAWSKIEGO 
Sklad glówny w ksiegarni E. Wende i S-ki. 


1916. 


Dodatek do ,Sprawozdañ z posiedzeñ* T. N. W. 


http://rcin.org.pl 


Geprüft und freigegeben durch die Kais. Deutsche Presseabteilung. 
Warschau deu 14. 7. 1916. T. M 1908. Dr. X 103. 


хем Ic 


Druk Rubieszewskiego i Wrotnowskiego w Warszawie, 


http://rcin.org.pl 


O ruchu planetoid typu *, (Thule). 


Wyznaczenie wyrazów typowych oraz przerwy. 


Wstęp. 


Wiadomo, że za charakterystykę danej grupy planetoid służy 
wymierność stosunku ruchu dziennego Jowisza, Marsa lub Saturna 
do ruchu dziennego płanetoidy. Gdy dana planetoida posiada ruch 
dzienny taki, iż stosunek ruchu dziennego planety głównej do ru- 
chu dziennego planetoidy wyraża się w przybliżeniu stosunkiem 
dwóch małych liczb całkowitych, jak np. $, 4,4, 2, ... etc., to zabu- 
rzenia w ruchu tej planetoidy, spowodowane planetą główną, są 
tem większe, im bardziej się ów stosunek ruchów dziennych zbliża 
do stosunku takich dwóch liczb. Największe zakłócenia w ruchach 
planetoid powoduje Jowisz, dzięki swej olbrzymiej masie; niektóre 
tylko planetoidy, których ruch dzienny jest współmierny z ruchem 
Marsa, podlegają dość znacznym zaburzeniom z jego strony, np. 
Hebe (6), Vala (131), Asia (67), Nysa (44), Lutetia (21) o dość 
znacznych ruchach dziennych i mimośrodach. Zaburzenia Jowi- 
szowe wywierane na tę grupę planetoid nie są duże, a toz powodu 
znacznej odległości Jowisza od tej grupy, Dopiero planetoidy gru- 
py š (z ruchem dziennym około 900”), o małym mimośrodzie, pla- 
netoidy grupy $ (z ruchem dziennym około 600") właściwie rozpo- 
czynają grupę planetoid silnie podlegających zaburzeniom ze stro- 
ny Jowisza. O ile stosunek ruchu planetoidy i planety głównej nie 
jest bardzo blizki stosunku dwóch małych liczb, to perturbacye 


О ruchu planetoid 1 
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nie beda zbyt wielkie, i wtedy zwykle metody obliczania perturba- 
cyj stosuja sie bez trudnosci. 

Gdy jednak ów stosunek staje sie coraz biizszy takiego pro- 
stego ulamka, jak np. dla przedstawicieli kazdej grupy planetoid, 
a wiec Hestii, Hekuby, Hildy etc., których ruchy dzienne najbar- 
dziej sa zbliżone do pewnej niewielkiej wielokrotności ruchu Jo- 
wisza—trudności przy obliczaniu perturbacyj stają się coraz więk- 
sze; zwykłe metody zawodzą, gdyż w wyrażeniach na perturbacye 
występują bardzo małe dzielniki, które przy całkowaniu jeszcze 
bardziej się zmniejszają. 

Przy pierwszem przybliżonem obliczeniu perturbacyj takich 
planetoid należałoby już uwzględnić wyrazy rzędu kwadratu masy 
zakłócającej, co, jak wiemy, niepomiernie komplikuje rachunki 
i często nawet uniemożliwia ich wykonanie o ile chcemy postępo- 
wać według metod klasycznych. 

Gyldćn podał właściwie jednę z metod przybliżonych, która, 
jak paru uczniów jego wykazało (Brendel, Backlund, Harzer 
i in.) ‘nie zawodzi nawet i wtedy, gdy stosunek ruchów jest bardzo 
blizki do prostego ułamka. Metodę Gyldćn'a, dość skompliko- 
waną i niezbyt dobrze przystosowaną do rachunków praktycznych, 
przerobił i uzupełnił Brendel. Metodę tę uczniowie Gyldćn'a 
zastosowali już do planetoid grupy Hestii, Hekuby i Hiidy — i za- 
wsze z dobrym rezultatem. Brendel i Kramer podali obszer- 
ne tablice, pozwalające w niedługim czasie obliczyć najważniejsze 
perturbacye planetoid typu Hestii i Hekuby i ich pozycye geocen- 
tryczne, tak że prawdopodobnie z łatwością będzie można zawsze 
odnaleźć te planetoidy. Wyniki dla ruchu Hekuby, otrzymane 
przez Harzer'a'), zostały do pewnego stopnia skontrolowane na 
innej drodze przez Simonin'a*). Przypadek planetoid typu Hil- 


dy 5) trudny z powodu bardzo blizkiego sasiedztwa Jowisza z te- 


mi planetoidami, zbadał Buchholz (obliczył on charakterystycz- 
ne wyrazy, ale nie podat jeszcze wyrazeñ samych perturbacyj). 
Planetoida Thule (279) nalezy do oddzielnej klasy planetoid: 


4 i jest, jak dotad, jedynym tej klasy przedstawicielem. Thule je- 


1) Harzer P. Untersuchungen über einen speziellen Fall... etc. 
2) Simonin M. Sur l'orbite d'Hécube, Annales de l'observatoire de 
Nice, t. IV. 
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szcze bardziej niz planetoidy typu Hildy moze sie zbliz yé do Jowi- 
sza; takich zbliżeń, jak wspomina Bidschof!), bywa po trzy na sto 
lat. Podczas takiego zbliżenia z Jowiszem, ponieważ Thule przez 
dłuższy czas znajduje się pod jego siinym wpływem zakłócajacym, 
droga jej ulega bardzo znacznym zmianom; np. dzięki zbliżeniu 
do Jowisza, które miało miejsce w 1912 r., średni ruch dzienny 
Thule znacznie się zmniejszył, jak to widzimy z rachunków Vilje- 
va?) Kobb?) przypuszczał nawet na zasadzie wzorów S.H. 
Darwin'a*), odnoszących się do orbit okresowych, że ruch Thule 
może z czasem stać się niestałym i droga jej może się nawet kiedyś 
przeciąć z drogą Jowisza. 

Dotąd nie badano ogólnie ruchu tej ciekawej planetoidy. 
O ile wiem, tylko Bidschof i Wedemeyer?) obliczyli lub po- 
prawili pierwiastki drogi Thule, a Viljev obliczył tylko zaburze- 
nia ruchu, odnoszące się do niedługiego okresu czasu, i to opie- 
rając się na danych orbity, otrzymanych przez Wedemeyer'a. 

Wydało mi się więc rzeczą inreresującą przystąpić na zasa- 
dzie metod, podanych przez H. Gyldćn'a, które z tak wielkiem 
powodzeniem dopełnił M. Brendel, do zbadania ruchu tej cie- 
kawej grupy planetoid typu Thule, których ruch średni dzienny 
jest bardzo blizki */, ruchu średniego dziennego Jowisza. Punk- 
tem wyjścia naszych rozważań będą dane liczbowe, odnoszące się 
do planetoidy Thule, podane!przez Wedemeyera. 

W niniejszej pracy zamierzam streścić ogólny zarys teoryi 
Gyldćn'a, obliczenia perturbacyj małych planet w modyfikacyi, 
wprowadzonej przez Brendel'a, i następnie wyznaczyć charakte- 
rystyczne wyrazy, odnoszące się do ruchu planetoidy Thule, oraz 
szerokość przerwy *) w pierścieniu planetoidalnym, która dla typu 
Thule dotąd nie została jeszcze oznaczoną. 


1) Sitzungsberichte d. Akad. d. Wissenschaften in Wien. Bd. 100 
Math.-Naturwiss. р. 937. J. Bidsc hof, Bestimmung der Bahn des Planeten 
(279) Thule. 

3 Viljew A. М. Bd. 195 

3) Kobb G. Sur un cas d'instabilité possible, Bull. Astr. 1902. 

') Darwin G. H. On periodic orbits, Acta Mathematica t. XV. 

5) Wedemeyer. Archiv der Seewarte XXXI Jhrg. 1909. 

5) Brendel. Theorie d. kl. Pl. II, p. 5. 
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I. Wyprowadzenie równań ruchu zakłóconego w spół- 
rzędnych Gyldénowskich. 


1. Wiadomo, że ogólne równania ruchu trzech ciał, które 
charakteryzują ruch zakłócony jakiejś planety dokoła słońca, są 
następujące: 


dx © 2Q 
qe Cm) = +1) ar 


Ë dy ? Y ys 29 
dË + k?(1-- m) уз =k*(m --1) oy 


2 ә 
= -- k*(14-m) 5 — k*(m 4-1) 4. 
gdzie, według przyjętego zwyczaju, przez k oznaczamy stałą Gaus- 
sa, przez m masę planety zakłóconej (przyjmując masę słońca za 
jedność), r promień wodzący tej planety; т’, 7' takie same wiel- 
kości dla planety zaktócajacej; v, y, 2 i z', у’, 2' odpowiednio dla 
każdej за to spółrzędne prostokątne, odniesione do pewnego stale- 
go kierunku w przestrzeni i których początek jest w środku słońca: 
t czas liczony w dniach średnich słonecznych, który upłynął od ja- 
kiejś epoki początkowej. Funkcya © jest to funkcya perturbacyjna: 
_ m jl ae +yy+z'| m |1 r 
й=ү M" ra | rr là moen 
gdzie 
А — (z — x -(y —y') (2 — 2) 2r*--r* —2rr'. cos Н 
zaś H jest to kat w środku słońca, pomiędzy promieniami wodzą- 
cymi r i r' obydwóch planet. 

Równania I należy przedewszystkiem odpowiednio prze- 
kształcić. 

H. Gyldćn analogicznie domyśli Hansen'a rozpatruje ruch 
planetoidy zakłóconej na swej chwilowej drodze oddzielnie od ru- 
chu samej płaszczyzny drogi w przestrzeni. 

Hansen okazał”), że istotnie można znaleźć takie wyrażenie 
analityczne dla ruchu zakłóconego planetoidy, iż elementy chara- 
kteryzujące położenie orbity w przestrzeni, a więc 9 i à i z nimi 
związane stałe, wystąpią oddzielnie od pozostałych elementów. 
Taki rozdział tych dwóch rodzajów charakterystycznych wielkości 


' Hansen Р. A. Auseinandersetzung etc. I. $ 3. 
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można osiągnąć przez wprowadzenie pewnego nowego układu 
spółrzędnych z,, y;, 2, mającego — jak zresztą i poprzedni układ 
xyz — swój początek w środku ciężkości słońca. 

Całkowanie układu równań różniczkowych I, w których pra- 
we strony byłyby zerami, dawałoby nam, jak wiemy, ruch po elip- 
sie Keplerowskiej. 

Równania te są następujące: 


2 ` 
zi + AP +m)=0 


ŻY 4. pt о 
(2) dp ^ prd m-—0 


dz, z 
de À rk -m)=0. 

Stałych całkowania w tym wypadku jest 6; są one dwóch 
rodzajów: jedne określają położenie drogi planetoidy w przestrze- 
ni, pozostałe zaś kształt i rozmiary drogi. Metoda waryacyi stałych 
dowolnych, jak wiemy, upoważnia nas, zamiast całek równania I 
(ruchu niekeplerowskiego, zakłóconego) wziąć całki równania (2) 
dla ruchu Keplerowskiego— pod warunkiem atoli, abyśmy stałe 
dowolne, które są w całkach równań ruchu keplerowskiego znajdu- 
ją, uważali jako funkcye czasu, które muszą spełniać równania |. 
ах dy dz 
dt * dt” dt 
obliczone na zasadzie układów równ. 2 i I mieć będą i te samą po- 
stać analityczną. Wszelkie spółrzędne, spełniające powyższy waru- 
nek (t. j. mające taki sam kształt analityczny w ruchu Keplerow- 
skim i w ruchu zakłóconym) Hansen nazwał spółrzędnemi ideal- 
nemi !). 

Niech te spółrzędne idealne beda z,, Yı 2,5 dostawy kierun- 
kowe osi z, oznaczymy przez a, В, т, Osi y, przez a,, f,, т, a osi 2, 
przez m, В», Ya — wówczas, jak wiemy; zachodzić będą związki 
następujące: 


(3) zy =a ++ By +yz; У Fc BY FET: à — arg + BY | T 
i odwrotnie 
(4) s= 0%, + 6444 + 921; Y= PtH PY, + 821: 

2— 123 + Yi + 1%, 


Wtedy wartości spółrzędnych x, y, z i ich pochodnych 


1) Hansen P. A. Auseinandersetzung etc. 1. $ 3: 
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aa, |- PB, т! =0 ++ r=l 
ча) aa, + BR TS =0 98814-11 1 
+ Babe te = 0 aż +B | 13 = 1. 
Aby nowe spółrzędne x,y,2, były „idealne“, wystarcza aby 
były spełnione związki następujące: 


da de її à Ei df, d п 
(5) ww dt ү rot +y dt Teg =" 
“Ж A 13 
Bu # -F dio A: 


Zwiazki (5) Mon sie do dwóch tylko warunków niezalez- 
nych; jakoz uwzgledniajac zwiazki (3) i podstawiajac w równ. (5): 


gdzie dla krötkosci oznaczyliśmy : 


Adi — N (4$ — — Y par 
Bdt = Y ad= — тав 
Cdt = Y pda=— Y ady 


(wielkosci te sa to dostawy kierunkowe chwilowej osi obrotu). 

Zrównania (6) widzimy przeto, ze we wszelkim ukladzie 
spólrzednych idealnych, odniesionym do osi ruchomych, o$ chwi- 
lowa obrotu ukladu schodzi sie z promieniem wodzacym planetoi- 
dy. Jeśli założymy, że stała 2, —0, wtedy płaszczyzna z, y stale 
zawierać będzie w sobie promień wodzący planetoidy i nowy układ 
spółrzędnych ruchomych x,, y,, 2, będzie najzupełniej określonym 
w przestrzeni. 

Ponieważ wprowadziliśmy warunek 1) 2, =0, mieć więc bę- 
dziemy 


(7) at + Bzy + 122 = 0; 
1) Nale2y przypomniec, 2e Hoene Wronski, podajac swoje nowe metod y 
obliczenia perturbacyj w „Réforme absolue du savoir humain*, rozpatruje per- 


turbacye, odnoszace sie do plaszczyzny drogi planetoidy oddzielnie od pertur- 
bacyj elementów jej orbity. 
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to równanie, jak to widzimy od razu, jest równaniem plaszczyzny 
przechodzącej przez promień wodzacy, a zmiennej z czasem f; 
płaszczyzna, zawierająca w sobie w czasie £ nietylko promień wo- 
dzący planetoidy dla tej chwili, ale jeszcze i styczną do jej drogi, 
spełniać będzie równanie, w którym zamiast x, y, 2 podstawimy 
z+ ах, ydy, z+ de. 

Wówczas muszą być spełnione związki nastepujace: 


(8) a dz + fady + tydz 0 
oraz 
(9) Eda + ydf, + 2d, = 0. 


Równania (7) i (8) dają nam warunek, który musi być speł- 
niony, aby płaszczyzna x, y, znajdowała się w płaszczyznie chwi- 
lowej drogi planetoidy; przy pomocy zaś równania: 

(10) ada, +8df, + (d, = 0 
wyznaczymy па tej płaszczyznie położenie osi 0z,. 

Pomnóżmy równania (4) raz przez da, df, dy, a drugi raz 
przez da, df, ‚dr, dodajmy i uwzględnijmy związki: (4a), 2, =0, 
oraz (10); otrzymamy: 

ado +ydf +12d =0 

хао, -ydg + 2d, — 0. 
Równania (3) zrózniczkujmy i uwzglednijmy równ. (9) i (11), 
wtedy otrzymamy, dołączając równanie 2, =0, równania nastę- 


pujące : 
dz,— ade +fxy --1рг 
(12) d= аа Вау таг 
dz, = 0 = adz + Bad у + 1,d2. 

Równania (12) pomnózmy przez а, a,, %,, później drugi raz 
przez В, Ви, fe, i Wreszcie trzeci raz przez 7, 11, Y, i dodajmy, to po 
uwzględnieniu związków (4a) otrzymamy : 

da = ada, + a,dy, 
dy = вах, + Вау, 
dz = таз, + т.001. 

Z tych równañ ро uwzglednieniu (4a) otrzymamy: 

dada + didy + dydz =0 
da dæ + ddy + айг = 0 


(11) 


igneus rein. "org. ToD 
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po uwzglednieniu których z (12) otrzymamy równania nastepujace: 
dx, =ad*x --gd^y +yde, 
(13) а?у, = s, dx +B,d?y+ rd’, 
d?z, —0, oraz: 
| dx = ат, + a,d*y, + dadu, + da, dz, 
dy = Ва?х, -- 8,d*», + dgdz, + dp dy, 
фа = Me + ydy, + ут, + di dy, . 


Funkcya perturbacyjna 9 da nam zwiazki nastepujace: 


29 29 2Q 29 
аа арта 
9 29 20 29 
(14) ou o эс TP" ay |^ 32 
го 
— =0. 
| 82, 


Jeżeli wreszcie równania I odpowiednio pomnozymy przez 
a, B, y 1 dodamy do siebie po uwzględnieniu dopiero co otrzyma- 
nych związków, to mieć będziemy !: 


2 ; ә 
o | erem eram 
2 
| w + k?(1-Em) a =hk?*(14m) у 
Układ równań II jest najzupełniej analogiczny do równania I 
i wyraża ruch planetoidy po chwilowej drodze?). W układzie rów- 
nań Il zamiast spółrzędnych prostokątnych wprowadźmy spółrzę- 
dne biegunowe r iv. Ponieważ płaszczyzna drogi planetoidy leży 
w płaszczyznie 2, у;, więc wskutek perturbacyj droga ta bynaj- 
mniej nie jest zamkniętą, i dla tego musimy liczyć kąty v nie tak, 
jak zwykle w elipsach od 0% do 360%, lecz od — oo do + co. 
Spółrzędne biegunowe, które wprowadzamy do równań II, są 
nam dane przez równanie następujące: 
(15) T, =7COS V; Yar Sin v; 


zaś rów. II przejdą w następujące: 


') Hansen. Auseinandersetzung I. p. 69. 
*; Harzer P. Untersuchungen etc. p. 8. 
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k? ier + Em) 


dr F 29 
(ш) | a) = (1+, 
dv dr ر‎ PRE (К 
| Tu. a tap (rm r % 


W tych równaniach v oznacza prawdziwą długość planety, 
zaś r jest to jej promień wodzący. Te równania są to równania 
zasadnicze Haesen'a, ale w spółrzędnych biegunowych. 

2. Równ. ШН. Gyldén przekształca przez wprowadzenie pe- 
wnych nowych zmiennych zależnych р, n, S oraz zmiennej niezale2- 
nej v „prawdziwej długości“; te wielkości są to tak zwane spółrzęd- 
ne Gyldénowskie. Zanim jednak to przekształcenie równania 
Ш wykonamy, przypomnimy w ogólnych zarysach bieg myśli G y l- 
dćna. Jeżeli w równaniach Ш przypuścimy, że > i A są równe 
0, to wtedy te równania przedstawiać będą ruch Une 
i, jak to wiemy, otrzymamy calki nastepujace: 


(1) = 1e cos (v — z) 


oraz stala predkosé wycinkowa 
, do 

"dt 

Poniewaz mase planetoidy zawsze mozemy zaniedbaé wobec 


masy slofica, która przyjmujemy za 1, wiec w ruchu Keplerow- 
skim jest zawsze: 


(3) a(1—e) <r <= a(1--e). 


(2) =} Va(1— e?). 


Jezeli а 01 =+0, to mieé bedziemy do czynienia z ruchem 


SOR GUN i M prawdziwa droga planetoidy znacznie sie róż- 
nić będzie od elipsy Keplerowskiej, jednak, jak to nam wyka- 
zują spostrzeżenia, dla czasu skończonego promienie wodzące pla- 
netoidy wahać się będą w pewnych skończonych granicach. Otóż 
Gyldćn wprowadził tak zwane krzywe periplegmatyczne 1), zawar- 
te między dwiema współśrodkowemi kulami, skończenie od siebie 
oddalonemi. Krzywe te są to poszczególne drogi planety, które się 
tworzą wskutek ciągłych zmian wielkości połowy osi jej orbity. 


1) Gyldén H. Orbites absolues r. I, 1. 1; id. Undersökningar... I. 
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Zmiany te, wedtug Gyldén'a sa nieograniczone w czasie, jednak 
zawarte są w skończonych granicach; wielkość а Gyldćn nazywa 
protometrem. Ponieważ w dalszych rozważaniach naszych 
a nie będzie absolutnym) elementem, przeto, idąc za Brendel'em, 
wielkość a nazywać będziemy wprost połową osi drogi (nie 
jest to jednak półoś chwilowej elipsy Keplerowskiej). 

Przypuśćmy, że planeta zakłócona porusza się po drodze, le- 
żącej w pewnej stałej płaszczyznie; niech 7 będzie jej promień wo- 
dzący, v—prawdziwa długość, liczona od pewnego stałego kierunku. 
Wskutek perturbacyj elipsa Keplerowska stale się odkształca, 
irównanie drogi, opisanej przez planetę zakłóconą, napisać moż- 
na, według Gyldén'a, w postaci następującej: 


a(1 — 1) 
4 Mom 
w I+p 
gdzie a jest to pewna $rednia wartošé promienia wodzacego, anip 
sa to funkcye prawdziwej długości, które przedstawić możemy 
przy pomocy szeregów trygonometrycznych; ч jest to pewna wiel- 
kość, bardzo powoli zmienna z czasem i porównywalna z mimośro- 
dem eliptycznym. W ruchu Keplerowskim prędkość wycinkowa 
była stałą, w ruchu zakłóconym, według Gyldé n’a, już tonie ma 
miejsca, gdyż ta prędkość wycinkowa waha się dokoła pewnej $re- 
dniej wartości. Przez analogię z ruchem Keplerowskim, mo- 
żemy więc napisać, że w ruchu zakłóconym zachodzi dla prędko- 
ści wycinkowej równanie następujące: 
dv KVall — nš) 

PIE — mx Ы 
(5) E aA +5 
gdzie S jest to pewna mała wielkość. Wprowadzmy funkcye S, р, т, 


w równania Ill. 
Drugie z röwnan Ш mozemy napisac: 


а | 430). 4 eQ 
dt (r dt =} (14m) Ten 

Połóżmy m —0; wtedy będziemy mogli napisać przez analo- 
gię, ze wzorem na ruch Keplerowski, równanie następujące: 


1) Absolutna drogą Gyldćn nazywa pewną orbite, w której elementy są 
to pewne zupełne stałe, i która różni się od rzeczywiście opisywanej przez da- 
ne ciało drogi, o wiełkości rzędu siły zakłócającej. 
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d ŻE: — 12) j=: 
dt 14-8 
wykonajmy wskazane rózniczkowanie, wtedy otrzymamy: 


a DE мас dS 1kVa 1 ТЕ 
` dv 


? (+9 (IFS "dv 2148 Y]? 
czyli, ješli wraz z Gyldén'em i Brendel'em!) oznaczymy: 
7? 99 
bis = —. DEM 
A 9 = TET u’ 


powy2szy wzór przejdzie w równanie nastepujace, które nam okre- 
sli S: 

1 din? 
qu WEERTUWA 94-5 1 LA 


Zwróćmy się teraz do wzorów (4) i (5); bez trudności znaj- 
dziemy: 
1 
dr _  kVa(l—w?) é 
d — 1-8 ` dv 
oraz 


1 1 
АР d; 
dr ka? (1 — n? 
d?" pm acp [^ qs т 09/93 1l 


Uwazajmy teraz równania Ill; w tych równaniach podstawmy 


; y dr dry | 
znalezioną dopiero co wartość na д: 9132 na ( d A , to otrzymamy: 


d? : 


(2) ай — 19) | gęs + (+8) 


— (1 +8) = — (+ 


ROSE | 
Równanie r= => ) da nam: 


') Brendel. Om Anvándningen р. 5; Gyldén, Undersókningar II. 
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ar! | 

r 1 СЬ. Ni Po, 1--р = l+p d^ 
dv? al — 1} lav: IP?" dv ` E — p?) dv 1—1 do 

co podstawione w równanie (a) da nam w rezultacie: 

d?p 2 
(09 Syp E ro ras c — 
1 2 dr?\? 
— (+8)P — fi — y A LR z) 


+8) , dr? | 

т 1—12 qa u Fre 2 
gdzie analogicznie do poprzedniego przez P oznaczylismy funkcye 
następującą: 


(Vois) P=r 80 


ey ` 

Równanie (V) daje nam właśnie wyrażenie na p w funkcyi v. 
Poniewaz P i (J sa to pewne szeregi trygonometryczne, wiec 

równania IV i V mieé beda postaé nastepujaca: 


dS А 
(1) de = У, An Sin (Av — Ban) 


d?p zł AA 
(2) dat — Be bx cos (A,v — B,). 


Zgodnie z oznaczeniami Brendel'a przyjmujemy, Ze w röwn. 
(1) i (2) wielkości @,, 6», Am Bn, В за to pewne stałe wielkości; 
w wielkościach a, , bw, В, (co jest istotne) masa ciała, wywołujące- 
go zaburzenia (w tym przypadku masa Jowisza), wystepuje jako 
czynnik. Te trzy wieikosci sg przeto porównywalne z masa zaktó- 
cajacego ciala. 

Calkowanie równania (1) da nam odrazu rezultat nastepujacy: 


(3) 8—a— Y i cos (Anv — B), 


gdzie przez a, oznaczyli$my stala calkowania. 
Z równania za$ (2) otrzymamy: 


1) Brendel. Th. d. kl. Planeten I. 
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b 
j + 2 : à 5 ` n A 
(4) p = k cos [1 —@—1 ] + N ] 8 r3 А COS (A, B,). 


Wyrażenia (3) i (4) zawierają jako dzielnik wielkość zwykle 
małą A, rzędu masy zakłócającej; wskutek całkowania we wszyst- 
kich wyrazach zależnych od X, wystąpią wielkości, w których masa 
zakłócająca zachodzi w mianowniku w wyższej potędze, wskutek 
czego wyrazy zależne od X, zostaną bardzo powiększone. 

W równ. (1) istnieją wyrazy pomnożone przez masę ciała za- 
kłócającego w pierwszej potędze; te wyrazy w całce (3) już nie bę- 
dą zawierać tej masy jako czynnika. Podobne wyrazy Gyldćn 
nazywa elementarnymi; ponieważ w przypadku tych wyrazów 
elementarnych spółczynniki A, będą bardzo małe (masa Jowisza 
145) przeto te wyrazy będą bardzo powolnie zmienne z czasem 
i według określenia Harzer'a!) są to wyrazy elementarne 
długookresowe. 

Jeżeli ruch planetoidy i planety zakłócającej jest tego rodza- 
ju, że stosunek ich ruchów średnich dziennych jest prawie równy 


p 


prostemu ułamkowi 0” to wtedy tez w równ. (3) wystapia wyrazy 


będące pierwszego rzędu w stosunku do masy planety zakłócają- 
cej; w niektórych przypadkach specyalnych mogą nawet być rzędu 
niższego od pierwszego, a mianowicie takiego, jak pierwiastek 
trzeciego stopnia z masy zakłócającej. Wyrazy, których rząd jest 
zawarty pomiędzy 4 i | —Gyldćn nazywa charakterystyczny- 
mi; ponieważ jednak w tych wyrazach czynnik 4, jest bardzo ma- 
ły, przeto i te wyrazy będą charakterystyczne długookre- 
sowe. Wreszcie w równ. (1) mogą występować takie wyrazy, dla 
których czynniki A, nie będą małemi wielkościami, wskutek czego 
po zcałkowaniu w (3) wyrazy te nie będą zwiększone; podobne 
wyrazy Są to wyrazy zwykłe. . 

Oznaczmy przez s, spólczynniki przy v w wyrazeniach wyra- 
zów elementarnych długookresowych (s jest rzędu тазу zakłóca- 
jącej); przez à, tę samą wielkość w wyrazach charakterystycznych 
długookresowych (à, jest wielkością małą, jednak znacznie więk- 
szą, niż masa zakłócająca); wtedy będziemy mogli napisać równ. 
(3) w postaci następującej: *) 


1) Harzer Р. Untersuchugen, p. 3. 
23) Brendel. Theor. d. kl. Planoten 1. p. 18 et sqq. 
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+ vum M 3 A, RYN ç 
(5) 8=и— v 5, COS (ant — A) — > z. cos (2,0 — On) 
as! А 
— x cos (А, — B,) 
x 


i w tem wyrazeniu À, juz nie jest mala wielkošcia. 

Zbadajmy obecnie w analogiczny sposób równ. (2). Przede- 
wszystkiem należy tutaj zauważyć, ze V1 — 8 == A, wskutek czego 
w (4) nigdy nie wystąpią dzielniki zerowe. Z pośród różnych wy- 
razów 2 М, wchodzących w równ. (4), musimy odróżnić przede- 
wszystkiem wyrazy, w których A, różni się od jedności tylko o wiel- 
kości rzędu masy ciała zakłócającego; w tych wyrazach À, będzie 
mieć kształt 1—5,, gdzie s,, jak i poprzednio, oznacza wielkość 
rzędu masy zakłócającej. W równ. (4) w odnośnych wyrazach wy- 
stąpią dzielniki tego samego rzędu, wskutek czego te wyrazy będą 
wyrazami elementarnymi; jednak będą one elementarnymi 
krótkookresowymi. Dalej w równ. (2) wystąpią także wyrazy, 
w których À, różni się od 1 o wielkości rzędu 2,, tak iż możemy 
napisać „== 1 — 2, w tych wyrazach występują dzielniki rzędu 2,; 
wyrazy te są to wyrazy charakterystyczne krótkookre- 
sowe. 

Wreszcie w równ. (2) są i takie А,, które bardzo się różnią od 
jedności; wyrazy, zawierające te wielkości, przy całkowaniu nie 
uzyskują małych dzielników, nie są więc one powiększone: są to 
wyrazy zwykłe 

Podobnie, jak to zrobiliśmy dla S, możemy napisać równanie 
na p w ten sposób, że każdy rodzaj wyrazów charakterystycznych 
występować będzie oddzielnie; p zatem mieć będzie postać nastę- 
pującą: 

f p=kcos [(1—«) v — T] + Zk, cos [(1 —0,)v — F,] + 
| +38, cos [(1 —2,) v—D] +20, cos [Anv — B,) 
gdzie, idąc za Brendel'em, oznaczylismy : 

Da ere b, А 
2e,— 2, --B —2(s,— €)— (o3 —?)' 

B = — E^ = bn | р" —. b,” 
" 28, — 8,2—8 21) <)’ 18 


5 Harzer. Untersuchungen, p. 3. 


(6) 


(7) 1—5=V1—$;k,= 
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Wyrazy charakterystyczne długookresowe i krótkookresowe, 
wyrazy elementarne obydwu kategoryj odgrywają zasadniczą rolę 
w rozwinięciach wielkości P i Q, które wchodzą w wyrażenia samej 
funkcyi perturbacyjnej. 

Harzer ') przejrzyście oznaczył w następujący sposób te 
wszystkie wyrazy, wypisując dla każdego rodzaju charakterystycz- 
ny kąt, którego sinus lub cosinus zachodzi w naszych wzorach. 
Harzer rozróżnia cztery postacie A, B, C, D wyrazów, charaktery- 
zujące się następującymi argumentami: 


(8) | А) 6.0 — Au; C) 8,0 ` Cz; 
| В) (1—0,)0—T y; D) (1—28,)v— Du. 


Wyrazy А) і В) są to wyrazy elementarne; C) i D) — chara- 
kterystyczne. Wyrazy kształtu A)i C) są długookresowe; zaś B) 
i D) krótkookresowe; wszelkie wyrazy, których argumenty nie są 
kształtu (8), są wyrazami zwykłymi. Należy zauważyć, że najniż- 
szy rząd względem mimośrodów drogi ciała zakłóconego i zakłóca- 
jącego dla wyrazów kształtu (A) jest rząd drugi, zaś dla wyra- 
zów (В) już pierwszy. Najniższy rząd wyrazów kształtu (C) i (D) 
względem mimośrodów, zależy od stosunku ruchów średnich y. pla- 
nety zakłóconej i planety zakłócającej. Niech stosunek ten mało 
się różni od stosunku dwóch liczb całkowitych p i q, wtedy ów 
rząd najniższy będzie | p — д | albo) | (p—q)—1 | |; w przypad- 


ku planetoid typu Thule, dla których : malo sie rózni od */,, wy- 


razy (B), (C) i (D) beda ju2 pierwszego rzedu, a wiec bardzo znaczne. 

Przy calkowaniu, które dalo poczatek spólczynnikom 5,, mo- 
że się zdarzyć, że niektóre z tych b, wypadną niezupełnie elemen- 
tarne, t. j. takie, które nie są rzędu 0 względem m”, lecz wyższego. 
Te wyrazy Gyldén nazywa wyrazami podelementarnymi. 
Wyrazy elementarne zawierają w sobie perturbacye wiekowe, w któ- 
rych według klasycznych teoryj czas występuje przed znakami 
funkcyj sinus i cosinus. Gyldén przez wprowadzenie wyrazów 
elementarnych uniknął tej niedogodności. 

Wielkością stopnia n-go nazywać będziemy wyrażenie, które 
zawiera jako czynniki niektóre z wielkości v, 4 Kn, kn" przyczem 
suma wykładników przy tych czynnikach jest równa n. Brendel 


1) Harzer. Untersuchungen, p. 3. 
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do teoryi malych planet wprowadza jeszcze okreslenie nastepujace. 
Wiadomo, że wielkością rzędu #-g0 nazywamy wielkość, która za- 
wiera jako czynnik masę zakłócającą w potędzę n-tej. Może się je- 
dnak zdarzyć, że podobny wyraz oprócz masy ciała zakłócającego 
zawiera jeszcze mały dzielnik rzędu 2,; otóż Brendel!) nazywa 
wyrazami rzędu istotnie n-go takie wyrazy, których wartość bez- 
względna może być porównywalna z n-tą potęgą masy m. Np. 
w równ. (7) zwykłe wyrazy (a także ich współczynniki 5,) są isto- 
tnie 1-go rzędu; charakterystyczne wyrazy (także i spólczynniki ß,) 
nie są istotnie pierwszego rzędu; natomiast iloczyn 5, X ё, jest 
formalnie drugiego rzędu, ale tylko istotnie pierwszego rzędu; 
idąc za Brendel'em oznaczać będziemy przez litery łacińskiego 
alfabetu wielkości istotnie pierwszego rzędu, zaś przez litery grec- 
kiego alfabetu wszystkie te wielkości, w których występują dziel- 
niki rzędu 0. . 

3. Przejdziemy teraz do zbadania ruchu planetoidy w jej 
chwilowej płaszczyźnie drogi. Gyldćn zakłada, że protometr or- 
bity a (w naszym przypadku jest to półoś orbity) jest zdefiniowany 
przez takie samo równanie, które daje nam związek w ruchu Ke- 
plerowskim pomiędzy a i n, a więc, 


_ kV] +m 
(1) = À 


w ruchu Keplerowskim. Według Gyldén'a a musi też takie 
samo równanie formalnie spełnić, tylko w tem równaniu n nie od- 
powiada już ruchowi średniemu po elipsie Keplerowskiej, jest 
raczej pewną stałą całkowania, którą Brendel nazywa stałą ruchu 
planetoidy, i która, jak zobaczymy, pozostaje w związku z wiel- 
kością n,, oznaczającą ruch dzienny średni. 


à dv Va(l—m) aq #7 (1 +5) 
2 М = =. 
Poniewaz równ. r R” "TAB da nam Дө ud 
więc możemy napisać: | 
dt „(b= 12} 
(2) "as U Xo? (1-5) 


równanie to będziemy musieli całkować. 
Wielkość p Gyldén rozbija na dwie części (p) i R: 


(3) p= (e) +R 


> Brendel. Theorie d. kl. Pl. I, p. 21. 
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w ten sposób, że w (p) są objęte wszystkie wyrazy typu (B), a R 
jest funkcyą pierwszego rzędu i zawiera wyrazy charakterysty- 
czne oraz największe zwyczajne. Równanie na (p) możemy 
napisać w postaci następującej : 
(p) =k cos [(— =) е —Г|-- ХК, cos [(1—о,) е —T„] 
albo 
(4) (p) = k cos (v — w) -+ Ў, cos (v — on) 
gdzie 
w= 4р0; e, =I, Sn t). 
Wprowadzmy do równ. па (р) funkcye q i П przy pomocy 
związków następujących: 
n cos I[ = k cos e + УК, cos w, 
n sin П==Ё sin e + Zk, sin w, 


z tych związków z łatwością znajdziemy 


(5) (p) = cos (v — II), 

gdzie q i II są to funkcye zawierające wyrazy elementarne długo- 
okresowe. Gdy funkcye 4, Il, R są już znane, to możemy przy po- 
mocy następujących związków obliczyć promień wodzący r, mia- 
nowicie równanie str. 4 jest: 


(4) s; 809) .. a(1— 12) = a(l — 12) 
I+p — 1+(p)+R l1+ncosv4+ В’ 
gdzie v — v — II. 


Z równania ostatniego wydzielmy znów pewna czesé r, ana- 

logicznie do tego cosmy zrobili dla p, mianowicie 
all a) 

(6) Оаа permet th 
Wielkość (r) w ten sposób zdefiniowana jestto według Gyldén'a 
„absolutny promień wodzący*. 

Równanie (2) z poprzedniego paragrafu daje nam związek 
pomiędzy czasem ź i miejscem plauetoidy na jej drodze; to równanie 
przekształca Gyldćn, posługując się analogią wzorów dla zmien- 


1) р jest rzędu masy zakłócającej i charakteryzuje ruch apsyd. 
Y 
O ruchu planetoid. 2 
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nych ruchu eliptycznego i tych nowych zmiennych, które on wpro- 
wadził jako charakteryzujące ruch planetoidy. W tym celu Gyl- 
dćn wprowadził nową zmienną e, analogiczną do anomalii mimo- 
środowej w ruchu eliptycznym, która jednakże występuje nie w wy- 
rażeniach promienia wodzącego 7, lecz absolutnego promienia (7). 
Posługując się wspomnianą analogią kształtu wzorów dla ruchu 
Keplerowskiego i zakłóconego, znajdziemy 

(7) (r)=a (1 — y cos s) 

stąd 


à 1—1? . 1-l-cos v 
a) sine — or sin v; b) cos шура бтр 


б и ES 
c) в =) {С во. 
Podobnie oznaczajac przez М zmienna analogiczna do апо- 
malii średniej w ruchu eliptycznym, mieć będziemy : 
(9) M=e—1 sin e. 
Postugujac sig formalnie wzorami ruchu eliptycznego, rozwi- 
niemy anomalie M w szereg !): 


M =з — qsine=v ++ YD B, sinn v, gdzie 
B, =— 2; B= + iT. В,——}— m 


А 2 
| B,=2 1; B,= x A anni. 


Oraz 


(8) 


(10) 
32 n 1--V1— 7 
Rózniczkujemy pierwsze równanie z (10), uwazajac, iz n jest 
stałą : 
1 de—n cos sde —dv--XnBncosnvdv, czyli 
(11) (1—1 cose) = 1-2, B,cos nv: 
1 — n° 


se STA © 
rózniczkujmy ten związek względem у: 
1— 12 VI — т> 


> MDC Ua D DD 
sin s z тову) sin v, lecz sine [cos Y 


') Brendel. Th. d. kl. PL I, p. 30. 
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podstawiajac te wartosé otrzymamy : 
de 1—12 de И1— 1? 


Бакр? | эшен E Darah 
LES hg (1 + qcos v) — dv 1—+ncos v ` 
Obie strony (11) (a) pomnózmy stronami przez otrzymane równa- 


V1 


A de 4 
nie dla de otrzymamy: 
(1— cos yA UP 
dv (14-37 co: -+7 cos v?) 
co znów podstawiajac w równanie (11) da nam rozwiniecie naste- 
puiace: 


3 
(12) бе er 
(1+-ncos у)? : 
Równanie (12) jest bardzo wazne a to dlatego, Ze dzieki samej defi- 
пісуі О yIdén osiągnął to, iz po prawej stronie (12) niema żadne- 
go wyrazu diugookresowego, podczas gdy z lewej strony znajduje 
sie funkcya o długim peryodzie. 

Zwróémy sie obecnie do zwiazku pomiedzy anomalia v i cza- 
sem {: mamy równanie (2), p. 16; w tem równaniu n jest to tak 
zwana stała ruchu planetoidy. Niech Z oznacza długość średnią, 
analogicznie do definicyi w ruchu eliptycznym, będzie ona okre- 
śloną w sposób następujący: 

(13) L=n+A, 


gdzie A oznacza długość średnią w chwili gdy ¿=0; wtedy będzie- 
my mogli zdefiniować nową zmienną czas zredukowany (t). 
„Czas zredukowany* spełnia związek następujący: 


(14) n(t) +A =v + ZB,sin nv. 


Na zasadzie znanego wzoru z trygonometryi możemy napisać 
uwzględniając wartość v: 


B,sinnv = B,sin nvcosn П — B,cosnvsinnIl, 
skąd otrzymamy: 


(14a) n Ó=1+-2nB, cos ne. © dv 
. dŠ _ «xd (B, cos n Il) d ( B, sinn IT) 
gdzie „, = m jp sin — F- 3 cos nv 
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jest to zresztą mała wielkość, którą przeważnie można zaniedbać 
(w naszym przypadku nie możemy tego zrobić). Ostatecznie więc, 
uwzględniając poprzednie związki, otrzymamy: 


| BL Ne p d 
dE = (12-3 cos v? "do 
Niech W oznacza róznice pomiedzy katem, odpowiadajacym 
prawdziwemu czasowi nt i kątem odpowiadającym czasowi zredu- 
kowanemu n(7)—mieć będziemy : 


(16) W =nt —n(t), 
a ponieważ 
L=nt--A 
zaś 
, A=v—n(t) +- XB, sin nv, 
wiec 
L-nt — n(t) +v+ XB, sin nv, 
czyli 


L-v-rFXB, sin nv + W. 
Wynika z (16) 


HV e n At) 
do "do "dv 


ai 


ale dla n A mamy wyrażenie (2), zaś dla n ¡y wyrażenie (14a), 
co uwzględniając otrzymamy: 
(17) ¥ (1 = E LER 2 — 1 _ ds 
— (1 +7 cos у)? (+= >) dv 
Rozwijając wyrażenie (17) według potęg R i S, otrzymamy: 
dW * : 
P —S—2R—2RS--3R? -- 3Sg*— AR* +.. 
+ (68 — 28 — 12R*-+-6R8S--...)q cos v 
= — зв (> S— 68+... nt cos 2v +... 
19 | 
+ 6R cosv+-(3 R — 8) y? cos Зу 
gdzie v — v— П. 
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Cala metoda Gyld&n’a jest zawarta we wzorach IV, IV*is, V, 
Vois i VL Zajmiemy się teraz zcalkowaniem tych równań i wyzna- 
czeniem wyrazów charakterystych i długookresowych w przypad- 
ku ruchu planetoidy Thule. Nim jednak do tego przystąpimy, 
musimy się zapoznać z rozwinięciem funkcyi perturbacyjnej 
w spółrzędnych Gyldénowskich. 


II. Rozwinięcie funkcyi perturbacyjnej. 


W równaniach otrzymanych poprzednio dla S i p wystepnja 
wielkości P i Q, które zawierają pochodne cząstkowe funkcyi per- 
turbacyjnej 9. W funkcyi perturbacyjnej występuje część 


a = (7272 — 2rr' cos H у i 


wskutek tego nie mozemy calkowaé powy2szych wyrazów bezpo- 
$rednio. Aby catkowanie wykonaé jestesmy zmuszeni rozwinaé 


A na taki szereg, którego wyrazy można całkować. Gyldén roz- 


wija funkcyę perturbacyjną na szereg trygonometryczny według 
wielokrotności kąta Н i potęg wielkości 7 i r”. Gyldén zakłada że: 


K = R,--2R, cos HR, cos 2H +... 


W tem rozwinięciu spółczynnik R, wyrazi się przy pomocy 
wzorów następujących : 


(1) 


a Г cos n 4 d à _4 1 Ji cos nds 
z. Vre —2rr'cosp T "hf uus m 


Oznaczajac a = s, i wprowadzając nową wielkość zmienną A, 


(2) R,= 


bardzo małą, określoną przez równanie: 


o as 
otrzymamy 
(4) = (i gne 


i związek (2) przejdzie w następujący: 
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LR Ik cosy ` 
Jezeli do tego ostatniego wzoru zastosujemy znane prze- 
ksztalcenie z teoryi calek eliptycznych, to w rezultacie otrzy- 
mamy: 


а 1 б cos nodo 
al 
0 


а! 


3 
2 x sin” odo 
== —;@”+1(1— А)? > 
(5) R z y” (1 af = 


ѕіп2 +02 sin?g * 
W calce (5) wyrazenie pod pierwiastkiem mozna napisaé 
w sposób nastepujacy: 


3 
2 ut f sin’ ede > 
т J V1 — atsin? g + afisintg — 
(6) z 
2 "T sin?" e d$ 
т : 2) sin? 
' 1 — o? sin? DET 
i J Te 


Gyldén, majac juz wartosé kazdego R z wyrazeniem (5) 
rozwija je według potęg A. W ten sposób znajdziemy: 


2 ? sin?"pde 
7 art! | == s == Te; À LN M. 
Wo” Miał MOT iren 


0 


gdzie za Gyldén'em !), oznaczajac przez 
(8) 8,09 = | 


T, 
0 


_ sin"ede 
$1 —adsin? | 5 
otrzymamy 

(1) 1 (3) 3 (5) 
| Tho = "+1; Yna == y t? EM 
(9) 


_ 1.3.5...(23—1) „m (224-1) 
|== N erte: №. 


1) Gyldén H. Undersökningar etc. II art. 38. 
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Spólczynniki Y»; zależą więc tylko od а. Wprowadzmy 
otrzymane wyrazenia z (7) i (9) do wzoru (5), wtedy znajdziemy: 


(00 B= AA maed aa" + 
Rozwiniecie (10) jest jeszcze zbiezne i w przypadku planeto- 
idy Thule. Podstawiajac wzór (10) w równ. (1) otrzymamy zadane 
rozwinięcie tylko części А 
wiem, Ze cala funkcya ERA 1) ma kształt następujący: 


funkcyi perturbacyjnej —wiadomo bo- 


at a > 
(11) a(g) mt —^ MES Hj=m EG a3(1— Xj cos H] 
Rozwiñmy w szereg wzór (11); w tym celu uwzgledniajac 


T otrzymamy ze wzoru (11) wzór nastepujacy : 


rozwiniecie dla A 


(12) a9) — mn, +2 n. fi E (Лу е |е H--2R, cos2H-L... 


mozemy ogólnie więc napisać, że szukane przez nas rozwinięcie 
funkcyi perturbacyjnej mieć będzie kształt następujący: 


(13) a (9) =2,-+-22, cos H+ 28, cos 2H+...2Q, cos nH = 


i =Q, + 9, cos n H, 
gdzie 


U n 
(14) 9, =m (1 — А) rio — Tni À + 1823 — Таз... | 


Тло Oznacza, że dla wszystkich wartości n, z wyjątkiem п —1, mieć 
będziemy 
a? 
tao = fo zaś dla n=1  (10— 110 — 9° 


Otrzymaliśmy więc rozwinięcie funkcyi perturbarcyjnej,w;po- 
staci, podanej przez Gyldéna. Jest tojednak jeszcze postać 
nie ostateczna, gdyż w rozwinięcie (13) zamiast wielkości A musi- 
my wprowadzić wielkości p, р’, n°, 12. 

Związki pomiędzy r, r' in, q' są następujące: 

kB a, Al z ДЭ. 


1+0 1--p' 


1) Z pominięciem zaburzeń samej płaszczyzny drogi. 
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a wiec 
A JM тта 1—70 1+p'\? 
0% к= (2) (7) A "17 | (=з) j 
Podnoszac À do potegi od 1-ej do 3-ej i uwzgledniajac wyra- 
zy trzeciego rzędu włącznie (a więc w krórych mamy př, р”, p?n, 
р'2 т," mp, etc.) znajdziemy: 
À = 2p — 257 — 3p?4 4pp!— p 212 — 2112 + 4p! — 6p*p" + 
+ 2pp°?— 4рт> + Apr? + 4012 — 4p. 
(16) + M = 4p— 8pp' ]- 4" — 129?-- 28p?p'— 20pp"?+ 4p4-Bpr? + 
+ 80172 — 8p'n? — 8p". 
№3 = 8p! — 24р%р' — 24pp'?— 8р. 
skad na а, (1 —à otrzymamy związek następujący: 
(pet 1 — py 
+" da 


OED N an 


(17) MT Eu =lm-+-(n+- 1)p' + 


2 
pd) p OE ay— REA рр" + 


n =” ED anna pt nn рт 


4 

n(n + D*g^q?. 

Otrzymany zwiazek (17) wprowadzmy do wzoru (13), wtedy 
otrzymamy ogólnie wzór nastepujacy: 


(18) a (Q) —2m' > Q (п, s, Sw pp" 127 12" cos n H. 

W tym wzorze należy pamietać, że dla п =0 mamy opuścić 
czynnik 2. W ten sposób ostatecznie otrzymaliśmy rozwinięcie 
postępujące według potęg n’, 1, p i p' dla całej funkcyi perturba- 
cyjuej. 

Spólczynniki 9, wyrażają się przy pomocy związku następu- 
jącego: 

(19) Q, = m! > Озу» ^p n mm. 


n, 5, #' v, v' 
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Te spólczynniki sa to funkcye liczbowe znanej nam wielko- 
ści a, a więc nie trudne do obliczenia. Harzer!) obliczył licz- 
bowe spółczynniki w wyrażeniu Q,, zaś Gyldén?) podał warto- 
ści tych spółczynników dla wielkich planet aż do wyrażeń zawiera- 
jących a’ a także i liczbowe wartości stałych niezależnych od a, 
które wchodzą w wyrażenia na ©. 

Ou... w funkcyi spółczynników liczbowych q,. wyrażają 
się w sposób następujący 3): 


(9 „оо = {n.0 
Quio = — N{n0 — 11.0 ; 
1 
g. = ED v E (nF Lima 


Q,4111— — (n+ 1) (0-2) ти — 2(2n4-5ruri1 — S 4122 
2,111 F— (n—1)n Ta—10 — 2(2n--lyt 14 — 814 12 
$4119 = — (Nn 1) тана — 21e 
Q, ii = — (n— 1-11 —2Ta—1a 
(20) ¿Rozo = 34 +4102 
2 

9,121 = Ся у сон = Теа — 10111 — 8112 + 

+4(3n + yt 23-2413 
9.1.20 = Ter ERE 2) Yni H (4-5) tia + 4186412 


n? (n — 1) 


Mn = 9 Tn—1o + (3n°+4n-|-2) Тина + 
FHAN + 5y(,—12 +) 24415 
(e. 120 = es Sio +(2n + 1) 'In—14 + 4-12 etc. 


W równ. różniczkowych na S i p wchodzą wyrażenia pocho- 
dnych cząstkowych funkcyi perturbacyjnej względem r i v—-a nie 
bezposrednio rozwiniecie samej funkcyi — przeto obecnie, zna- 
jac wyrazenia rozwinietej juz funkcyi perturbacyjnej mozemy zna- 


1) Harzer P. Untersuchungen pg. 22. 

з) Gyldén. Orbites absolues t. 1. 

3) Brendel. Th. d. kl. РІ. I, pg. 52; Harzer. Untersuch., pg. 29; 
Kramer. Th, d. kl. Planeten, pg. 9. 
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lezé i związek pomiędzy P i Q oraz funkcyami Q,,,,,5; W ten spo- 
sób wielkości P i Q staną sie znanemi funkcyami liczbowemi 7,2. 
co zatem idzie z, a wiec ich obliczenie, choé klopotliwe, dla danej 
planetoidy'nie przedstawia zasadniczych trudnosci. 

Mieli$my równ. nastepujace: 


2Q 7? 29 
=p: = з . 
P=r =; oraz () (1—1) 20 
Poniewaz 
_ a(1—1') 
1--p 


TN r 0 1—1р 20 
— a(fl—n?) 83» (14-р)? w 
W tem wyrazeniu pierwszy czynnik rozwiümy w szereg po 
tegowy, potem zrózniczkujmy wzgledem v; otrzymamy wtedy 


to znajdziemy: 


© 
Q= — (1 — 2р 4 302 — ...) ZEN nas, pp" 1 y?" sin nH E 
(bo wiemy, ze H zalezy od v) co mozemy napisac: 
ә 
(21) О = — 27 п Que pen?" sin п Н = 


(dla n = 0 należy opuścić czynnik 2) — gdzie przez (,,,.,, Ozna- 
czylismy pewna funkcye Queso а Więc funkcyę wielkości Yno, któ- 
rą można adi na drodze rachunkowej: postać tej funkcyi jest: 


(22) ЯЕ Миа = О ES A 20,4— Ls. y d ЗӨ» oa v! 


bye A analogicznie znajdziemy wyrazenie 
Р= т = — (1 _) a 
Po wykonaniu dziala i podstawien wskazanych, jak to 
zrobilismy dla Q, otrzymamy rezultat nastepujacy: 
(23) P=232 Pay pr: р? т." cos n H, 
gdzie przez P,,,,, oznaczyliśmy wyrażenie następujące: 


1 


(24) m 


Fa sa! vy! = — (8- r 1)R scis + (8 + Mans 1..»—1.v — 
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Wielkości P są tylko funkcyami liczbowemi wielkości Q,...,, 
а więc Tns, przez co także te wielkości sa znanemi funkcyami a. 


Brendel?) Harzer?) podali wartości Pass i Que: 


Pavo =— 9,10 Quoo = noo 
Paio = — 29,20 Qn.10 = 9.10 — 29.00 
Pazo —=—3Qus0 Quao = Luzo — 201.0 +200 
Ро 1— Qua нод = Q, toa 
Р, 210. — — Q, 14 Qi Q, 104 
Pu-i00= >= Q, 100 Qut00 = Q,1100 
Роло = —2Q20 Quoi = P101 

(4)! Pa =— Buu Фоо = 9100 
Pioo =— 910 Q. 444 — 29 
Pr 111 =—2 Q, 21 Qn—1.1 1 = Q, 111— 29, 101 
Pai Z — ZB»; 21 On+ı 1.02 4110 — 29-4100 
Pa+11.0 == Q, {120 Qn—1.10 =Q —1.1.0 — 20, 100 
Pu = — 29. 120 Quia = 911 — 29:0, 
Розо == 3 20.3.0 Qı 10 = Qi10 23 291.00 
Pisi = —29,s; 

(Pao = — 2 Rizo 


W wyrazeniach na P i Q wystepuje kat H. Musimy obecnie 
ten kąt wyrazić przy pomocy prawdziwych długości v Го’ plane- 
toidy i planety zakłócającej; jak dotąd tak i teraz rozpatrujemy 
tylko ruch planetoidy w płaszczyznie drogi chwilowej, pozostawia- 
jac do dalszej części niniejszych rozważań, wyznaczenie „wyrazów 
funkcyi perturbacyjnej, zależnych od nachylenia i zmian położenia 
płaszczyzny samej drogi. W ten sposób znajdziemy po zaniedba- 
niu drugiej potęgi wstawy kąta nachylenia pomiędzy drogą plane- 
toidy i planety zakłócającej: 

(25) ?) H=v—v' 


Po wprowadzemiu tej wartości H w wyrażenia па Рі Q 
otrzymamy wzory, w których występować będą obecnie długości 
prawdziwe v i v'. W równaniach, z których możemy otrzymać 


1) Brendel, 1. cit, pg. 56. Harzer, l. cit, pg. 30, 31. 
*) Masal H. Formeln u. Tafeln etc., pg. 5. 
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S, p, W—zmienną niezależną jest prawdziwa długość v — więc aby 
módz zcałkować równ. dla Р iQ musimy wyrazić v' przez v. 
W tym celu zwróćmy sie, idąc za Brendel'em, do równania (13) 
pg. 19, z którego otrzymamy związek pomiędzy położeniem plane- 
toidy na jej drodze i czasem. Znane wzory ruchu Keplerow- 
skiego dadzą nam dla planetoidy: 


(a) nt+A=v— 2 sin v + in sin 2v —...-- W 

i analogicznie dla planety zakłócającej: 

(b) n't 4- A'— v'— 2y' sin v'--*/, т sin 2v'—...4- W” 
gdzie у=0—П, zas v'—v'— IT. 


, 
Mnozac wyrazenie (a) przez E=" otrzymamy: 


(26) n't=y0— u A—2p7 sin v+ ip v sin 29v+-...W y. 
Z równañ (b) i (26) przez porównanie znajdziemy: 
v'= (po — pA — 291 sin V--3/, pr? sin 2 v... Wu.) 
+ A'+ 21 sin v'’+...+W". 
Wprowadźmy oznaczenie następujące: 
(27) B=N—pA i G—yXB,sin nv — ХВ, sin nv' 


to otrzymamy: 


(28) v=w + B+ G -yW — М’ 
a wiec: 
(29) v—v= (1—p)0— B— G — узу =w —– С. 


(Możemy opuścić W' jako małą wielkość, gdyż w tym wyra- 
zie jest zawarty wpływ Saturna na ruch Jowisza, zresztą ta wiel- 
kość jest nam znaną z teoryi Jowisza), gdzie: 


(30) w=(1—p)v— B—pW. 
, Polózmy dalej 

(31) V'=—w, + G+ v, 5 

wtedy 


(32) G = —2y«sin v +27 sin v,'+iun?sin 2v-]- 27"? sin 2v, +... 


1) Brendel. Th. d. kl. PL I, pg. 63. 
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Aby przedstawić cos n(v —v') jako funkcyą v, rozwińmy 
przedewszystkiem, uwzględniając związek (31), wyrażenie 4/ sin v, 
w szereg według potęg тү: i 
(33) y sin v,'—-—"' sin [w, — G — v,j =—7 sin (w; --V;) + 

+16 cos (W, — v1) —... 
gdy w (32) uwzglednimy wyrazy tylko do pierwszego stopnia 
włącznie, to otrzymamy: 
G= ——P2ku sin v — 21 sin (w, — vi) 
Podstawiajac te wartosé w (33) znajdziemy: 
т sin v,'— — sin (w,— У) — 29.74 sin v cos (w,—v;,) — 
— 28.12 sin (W, — vj) cos (W, —v,) =— sin (w, — У) — 
— pnn sin (v— V, +W) -Henn sin (w, — v — vi) 
— 1" sin (2w, —2v,) +... 

Analogicznie postepujac otrzymany dalej: 

— in"sin2v',— 11" sin 2(w, — v, + G) - ix? sin 2(w, — Vi) +... 

Podstawiajac te wyrazenia we wzór na G, otrzymamy w re- 
zultacie: 

G=—2p7 sin v—27 sin (w, — V4) + Fp? sin 2v — 
(34) 1 — 2p nn sin (w, + V—V,) 3-217 sin (W; — v — v)— 

— łn? sin (2w, — 2v +... 

Poniewaz ze wzoru (29) mamy 

v—v—w,—G, 
wiec 
A > т? С? 

cos n (v — v') = cos n W,+nG sin n w,— ^ COS RW; +... 

Możemy więc, uwzględniając wyrazy do III stopnia wiącznie, 
wstawić wyrażenie (34) w to ostatnie rozwinięcie. Ostatecznie 
otrzymujemy cos n(v—v') w postaci szeregu, w którym zmienną 
jest tylko v. Dla sin n(v — v') analogiczny wzór mieć będziemy 
z tą różnicą, że zamiast cos wystąpi sin. 
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Bedzie wicc: 


[= n(v—v') = in? Wid- DOM i т утә 7 
‚ Cos 
+k" sin [(n—Dw, +y] +=, 7 sin 95 [m--1)w,—v.] 
cos cos 
TES) sin пу, + ke n? sin (nw,--2v) 
(85); x Co (nw,—2v) + km S i D) wv v] 


Th. me Se 04D) w, vv Hi itt АЕ 


cos 


, 
Heun sin 


[m--1)w,—v м]? © 


sin ” 


l TÉ ү = [(n — 2)w; +-2v,] Ha Si sn [(04-2)w, —2v,] 


gdzie: 

k,—np; kv=—nt; k,— —n; k=--"n; Ку np; 
mo MT RECE 2) ais M S N: 

Msn 9 g e? AE) 2 Tg" 

k,— —п(п — 1); k,— 4-n(n4-1)p; k= +n(n —1)ц.; 

k,,=—nln--1)p; k= n’; n=l LÉ "i 
ELI ME 

ме 9 т 8 5j. 


Gdy juz wykonalismy rozwiniecie En - n H według potęg v, 


to wtedy wyrażenia (21) i (23) pg. 26 napisać będzie można w po- 
staci następującej: 


(36a) p=2% Passw pop? n°?" cos п (0 — v'), 
oraz 
(37) Q=—2Š Q, sn pr po nn?" sin n(v—v'). 


(š oznacza, ze dla n=0 nalezy opušcié czynnik 9). 
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W tych wzorach wydzielmy z wielkošci p czešé elementarna 
ksztaltu B oraz zwykla wraz z charakterystyczna; pomiedzy temi 
róznemi czešciami p mamy, na zasadzie poprzedniego, zwiazki na- 
stepujace: 
(38) p=ncosv+-R oraz p=(p)--R, 
gdzie, jak zawsze, (p) oznacza część elementarną p postaci B; zaś 
w R są zawarte wyrazy zwykłe i charakterystyczne. 

Ponieważ: 


и + Rs). 8--(5)e) . +... Be, 


więc gdy w tym wyrażeniu na p uwzględnimy tylko część charakte- 
rystyczną i podstawimy w (36a), to otrzymamy rozwinięcie postaci 
następującej: 
P=P,-+- P,R+-P,R*-+- P,R° +... P, Р +... 
(39) | кале á 
| P i= 22 Passw (о)* . (py > cos (v— v'). 

Dla Q, mieć będziemy analogiczne wzory, tylko w tych wzo- 
rach zamiast cos (0—v') wszędzie znajdować się będzie sin (v—v'). 
Zatem według drugiego wzoru ze związku (39) mieć będziemy: 


R, [Y P (0.9900 TE „aj TR P n.0.100 Ch 


Qíu.0.0) 0.0 m.1.0).0.0 Zi (Jin.0.1)00 
P 1 7 Р 102500 

ЕЕ У (n.2. e p)? > Pas o. p' (n.0.2)0, 12 
Re + Qi RS. qd Qn.0.20.0 P m 
R P COS 

+ y ("w.0.0)1.0 (m.0.0)0.1 y? LAD 2 n(v —w) 
Qo 011.0 Oe sin 

m | P (1000 ый ju M * Pio 

En Y ы (u.2.0)0.0 لے‎ ©0(n.1.1)0.0 


x RT pu ` n(v—wv) i t. d. 


Aby otrzymać ostateczną postać rozwinięcia funkcyi P i Q — 
musimy wykonaé jeszcze jedno przeksztalcenie. Kazda funkcya 
S, R, W, zawiera bowiem czesé zwykla, krótkookresowa, oraz diu- 
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gookresową, wyrazy długo i krótkookresowe elementarne, takież 
wyrazy charakterystyczne i wreszcie wyrazy zwykłe. 

W ten sposób funkcyę W możemy napisać w postaci nastę- 
pującej: 
(40) W —1v-4- Wi+ W,-4- W, 

Część yv jest to część wiekowa. Wu oznacza wyrazy długo- 


okresowe, W;,—wyrazy krótkookresowe, W. zaś wyrazy zwykłe. 
Według Brendel'a oznaczmy 


(41) W;+W,=K, 
to wtedy 
(42) W=:0-+W¿+K, 


ale we wzorze (30) pg. 28 mamy 
(43) mw,—n(l—p)v —n B—nuW —n(1—p)v—nB—nyW ,—npK, 
gdzie dla któtko$ci oznaczylismy 
(1 + (= р. 
Jezeli jeszcze oznaczymy 
(1 —p)0— В — рМ. 20, 


to wtedy !) 
w,—w-—K. 
Wtedy dla 
sin nw, 
otrzymamy wzór: р 
2 
(44) sin nw, = sin nw— ny cos nw - ч: EE sin nw... 


Jeżeli teraz uwzględnimy wszystkie wartości Po, P,, Py... 
i złączymy razem wyrazy zależne od tego samego argumentu, 
oraz uwzględnimy związki (44), to otrzymamy w rezultacie rozwi- 
nięcie w szereg nieskończony według potęg т, тү, R, K cząstko- 
wej pochodnej P funkcyi perturbacyjnej, w którym występować 
będą tylko argumenty v, v,, oraz w, i pewne liczbowe spółczynni- 


ki (n.s.s')v.v' zależne jedynie od s, = a, które dla każdego po- 
szczególnego przypadku zawsze będziemy mogli obliczyć. 


1) Brendel, I c. pg. 69. 
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W rezultacie więc otrzymamy dla P rozwinięcie?) podane 
przez Brendel'a w postaci następującej: 


(45) P=EB, оо cos mw-|-E BC ncos(nw+v)+ EB) ncos(nw —V) 
+ УВ) cos(nw--v,) ® ВС) y! cos (nw — vi) 
+ Z В, >01? cos nw 
-Е®В y cos (nw--2v) -+ ХВС 2 n° cos (nw —2v) 
+ XB "hyn'cos(nw-|- v-v) EBE u’ cos (mw +v—w, ) 


“li 


EB, m cos(nw—v+v,)H2B mcos (nw—v—v,) 


n.1.1 


4-Х В, оз cos n 10 
HE Ba cos (т--2у)--ЕВ s 1 "cos (1102 v,) 
des 


+ В УВ, cosnw+ EBE? cos (nw + v) 
+YB- «cos (nw — v) J- X BT qf cos (nw + v) 


n.0.1 


+ SB; cos (mo — v) +...) 
+pK 4 Xn B,,, sin nw+-Zn BE” y sin (nw-- v) 


n.1.0 


+ Xn Bi «sin (nw —v)-- En BU" « sin (nw-+-v,) 


n.1.0 n.0.1 
+n BO) w sin (nw —v,) +... ) 
+ R?{ EB}, cos w+}... Fp RK {En В! sin me^...) 


— wet v В, cos moy °} 


') Brendel, 1. c. pg. 70. Harzer. Untersuchungen etc., pg. 30. 

2) Wskaźniki przy spółczynnikach A i B są tak dobrane, że u dołu 
pierwszy wskaźnik odpowiada czynnikowi przy w, następny daje nam wykła- 
dnik +, trzeci zaś wykładnik * — górne wskaźniki dają czynniki przy v oraz 
wykładnik przy E (confr. Brendel 1. c, pg. 71). Dla wartości n=0 niektó- 
re ze spółczynników À i B są zerami (confr. Brendel, 1. c., pg. 71 i 79). 


O ruchu planetoid. 3 
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Dla Q za$ mieé bedziemy analogicznie: 


(46) Q=2 Anoo sin nw--X 4907 q sin(mo+v)+2 AC «sin (mw v) 


#.1.0 


+ ХАС q sin (nw + v,) FIA 4! sin (mw — v,) 


n.0.1 
+A, ,, "^ Sin nw 
ХА nf sin (nw-+-2v) LAC? y? sin (nw — 2v) 


n.2.0 
+EA qw! sin (mov + vi) + X ACC qu sin (nw +v — vj) 
+ЕЎА M sin (nw—v +v) J-E АС, am sin (mo v-- vi) 
+A Lou 7? sin nw 
-- 2455, 1? sin (mw-]-2v,) — --E AC q* sin (mw—2v,)+-... 


no2 1 
+ RATS, sin nw + 2 4715 nsin (mwv) X А1210 qsin(nw —v) 
+24, W sin (nw+ vj) FIA 1! cos (mo —v,)-+... $ 
—pK À Xn Anoo COS nw+ Xn An cos (nw-- v) 


+ En AC2 q cos (nw —v) + Xn AU" чү cos (nw+-v,) 


n,10 n.0.1 
+ En AT cos (nw—v,) +-...) 


+ R? { ХА?) sin nw...) — RK (n AY, соз тш} 


"0.0 2.0.0 


2 dy i LY 
Ея 5 Anoo sin зе 


Wzory (45) i (46 rozwinięte do 3-go rzędu włącznie posłużą 
nam do wynalezienia wyrażeń analitycznych dla ruchu planetoidy 
Thule (279). Wyrażenia (45) i (46) będą wystarczające dla nasze- 
go celu, ponieważ nachylenie orbity Thule względem orbity Jowi- 
sza jest bardzo małe. 

We wzorach (45) i (46) występują spółczynniki À i B opa- 
trzone wskaźnikami; te spółczynniki są w dość prostym związku 
z wielkościami Pyss i Que, które, jak wiemy, są to funkcye Ques w, 
przeto spółczynniki Ai B ze wskaźnikami dadzą się z łatwością 


1) l. c. pg. 33. 
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obliczyć јако funkcye liczbowe n.s, czyli wielkości z. znanej nam 
z elementów drogi uważanej planetoidy i Jowisza. 
Związki pomiędzy À i B oraz P i Q są następujące”): 


Anoo = — 2nQ«oo 
Ан, = — M Quio 
A00 = — Qaso 
АФ) = – пі. np Q.cof 
Aj = — NiQuio Sm Оно! 


A = (а-да DO! 
Ai = — (n—1) |Qs_101 +2 (n— 1) Qu—1.00 


AT —=- 1120.20 + 2ng Qn. 1.0) 
AN 120030 — 24 Quo! 
Ad = (н) Qaa — 20 +1) ало 


А129 — —(и—1) (Queras + Ua lQ 


В, оо = 2Р„о„ 
н; = 2P.40 
i = 2Pn20 


ВО = Pa1o2ny Pavo 

Hp = P, — 2 Poo 

BU = Р„нол— 2(n+1)P,+ioo 
BO = Pa-101 +2(n—1)P, 100 
Bii" = 2Р,зо-- ть Pro 
В = 2Pu20—2np Рило 

Bi” = Paris — 2(n--1) Prae 


Bo = Р, лал + 2(n—1) Р. 110. 
i) Brendel. Th. d. kl. Pl. I Кар. V sparsim. 
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ПІ. Wyznaczenie wyrazów charakterystycznych 
i elementarnych w rozwinięciach odnoszących się 
do ruchu planetoidy Thule. 


Nastr. 13i 14 podaliśmy ogólny kształt równań różniczkowych 
wyznaczających wielkości S i p; równanie dające W jest analogicz- 


ers 
ne do równ. na S, mianowicie W — C, + X E sin (v — Bu). 


W powyższych równaniach wielkości а, b, c, są to bardzo 
małe stałe wielkości; À,. jak to wiemy, jest mały dzielnik. Mo- 
że się zdarzyć, iż pomimo, że wielkości a, 5, c, są małe, je- 
dnak wskutek całkowania wystąpią jako dzielniki A„, które w nie- 
których przypadkach mogą mieć takie wartości, że dany wyraz, 


w którym wystąpią po całkowaniu czynniki : albo 1 moze 


E 
sie stać albo bardzo małym blizkim zera, albo też bardzo dużym +). 

Gdy À, jest bardzo małe, blizkie zera; wtedy argument cosi- 
musa jest bardzo małym kątem, z czego wynika, że okres takiego 
wyrazu, zawierającego w sobie jako dzielnik bardzo małe A, jest 
bardzo znacznym. — Wyrazy zależne od takiego X, będą wyra- 
zami długookresowemi. 

Gdy zaś x, jest blizkie 1, wtedy okres takiego wyrazu jest 
niewiele różny od okresu zmienności kąta v, a więc od czasu obie- 
gu planetoidy. Te wyrazy nazywać będziemy wyrazami krótko- 
okresowemi. 

Jest rzecza oczywista, iz okresy krótkookresowych wyrazów 
są blizkie okresu obiegu planetoidy. Jak to już poprzednio wspo- 
minalismy w wyrażeniach odnoszących się do S, wskutek całko- 
wania tylko wyrazy długookresowe zostaną powiększone, podczas 
gdy wyrazy o krótkim okresie nie ulegną zmianie; odwrotnie zaś 
dziać się będzie, gdy całkować będziemy wyrażenia dające namp— 
tam długookresowe wyrazy nie zmienią się, natomiast krótkookre- 
sowe będą powiększone. Brendel zauważył*), że w przypadku 
współmierności bardzo blizkiej pomiędzy ruchem planetoidy i pla- 
nety zakłócającej, wyrazy nawet drugiego rzędu bywają jeszcze 
bardzo znaczne. — Z drugiej strony zasadnicza myśl Gyldén'a 


) Brendel. Om amvändningen pg. 6. Ludendorff. Die Jupiter- 


stórungen d. kl. Planeten von Hecuba-Typus pg. 4. 
2) Brendel. Th. d. kl. Pl. I Kap. VII sparsim. 
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polega właśnie na tem, ze w wyrażeniach przyblizonych, jakie 
otrzymujemy na P i Q porządkujemy wyrazy nie według rzędu ani 
stopnia, lecz jedynie według rzeczywistej wartości tych wyrażeń, 
które one osiągają w wyrażeniach całkowych. 

Dla tego też przedewszystkiem musimy w rozwinięciach na 
S, p, R, W zająć się wyznaczeniem krótko i długookresowych wy- 
razów dla specyalnego przypadku ruchu planetoidy Thule (279). 
W wielkościach S, р, W, R, zgodnie z Brendel'em, oddzielmy 
część długookresową, część zwykłą i część krótkookresową w ten 
sposób, że pod odnośną liczbą dopisujemy wskaźniki d, 2, k; przez 
(p), zgodnie zresztą z tem co poprzednio pisaliśmy, oznaczymy 
w wielkości p wszystkie wyrazy elementarne postaci В; te wyrazy 
zresztą, jak to już zauważył Brendel, są krótkookresowe !) (są 
one oznaczone przez (p)). Mieć więc będziemy: 


(1) S= Sa + Se +- 8,; R= R. + К, В. 
p= (e) + Е; w—=10+WitK, 


gdzie К = W,+ W., zaś yv oznacza wyraz wiekowy. 

Przejdźmy do naszego przypadku specyalnego; stosunek ru- 
chu dziennego n' Jowisza do ruchu dziennego średniego n, plane- 
toidy (279) Thule wynosi: 


IT > _ 3-3, 
(2) AA z” 
gdzie à, jest bardzo mała wielkość. 
Zatem 
, 14-2 
(3) {м=р 


a ponieważ ogólna postać argumentu, wchodzącego w nasze roz- 
winięcie jest 

(4) nw=n(l —p)v— nB—unW;i=n(l — в.) —nB—nywW,;, 
to w wyrazach krótkookresowych spółczynnik przy v w tym argu- 


mencie jest blizkim I, w długookresowych zaś różni się on mało 
od zera. Jeżeli więc chcemy znaleźć w rozwinięciach na Pi Q 


5 Brendel, I. c, pg. 19. 
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wyrazy krótkookresowe i dtugookresowe, to musimy tak dobrac 
czynnik п, mnoznik przy v, aby argument cosinusa, wzglednie si- 
nusa odpowiedniego kata byl w pierwszym przypadku blizki je- 
dnosci, w drugim zaś blizkim zera. Uwazajmy zatem nasze rozwi- 
nięcie (45) i (46). 

Ponieważ w obu tych wzorach wskaźniki odnośnych spół- 
czynników są te same, zaś same spółczynniki są inaczej oznaczo- 
ne, przeto chwilowo wystarczy nam w zupełności zbadać tylko 
spółczynniki zawarte w rozwinięciu dla P, gdyż spółczynniki roz- 
winięcia dla Q otrzymamy ze spółczynników P, zamieniając tylko 
symbol À na symbol В. 

W rozwinięciu (45) następujące wyrazy l-go rzędu 0, 1 i 2 
stopnia będą długookresowe. W 1, 6, 13 wyrazie składnik argu- 
mentu zależny od v jest blizki 0 dla n —0, a więc długookresowe 
funkcye będą przy Ву, Boos: В}, Ву; dla n=4 3-ci i 5-ty 
wyrazy są długookresowe. Odnośne tedy spółczynniki dla tych 
wyrazów będą: В, у, Ви,. Podobnie 8, 12, 15 wyraz dla n= 8 
będzie blizki zera; spółczynniki tych wyrazów będą odpowiednio 
Bio» Bua) Po. Wyrazy krótkookresowe tego samego rzędu 
` co poprzednie, a więc I-go, są te, dla których argument zależny od 
v będzie blizki 1; te wyrazy są następujące: drugi i czwarty wy- 
raz dlan=0. Spółczynniki tych wyrazów będą; Bi m. ; da- 
lej З i 5-ty które nam dadzą B,, i i m. Biorąc 1 —4 otrzyma: 
my krótkookresowe wyrazenia dla 1, 6, 9, 10, 12 wyrazu, dla któ- 
rych spólczynniki będą B,,,, Biso Bios B. |o Dae 
Bii, Bios; dla n=8 argumenty powodujące wyrazy krótkookre- 
sowe występują w wyrazach 3 i 5-ym, dla których spółczynniki be- 
dą Bus i в wreszcie dla n=12 takie same wyrazy wystepuja 
w 8, 12, 15 wyrazie rozwinięcia, skąd otrzymamy: Bso, Bau» 
ns . Wyczerpalismy tedy wśród wyrazów |-go rzędu już 
wszystkie, które dla Thule (279) będą krótko i długookresowe. 

Wypiszmy te wyrazy: 

Wyrazy o długim okresie: Wyrazy o krótkim okresie: 
Bi 1 cos v 


2 +1 
Bo20" Boo 1 СОЗУ, 
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Baoa т? 
B® mn! cos (v — v.) 


0.1.1 


Ву тпр cos (v, — v) 


.1.1 


ВС 4 cos (4w — v) 


4.1.0 
Ву, cos (4w —v,) 


—2) 
8.2.0 


n? cos (8w—2v) 


°В m cos (8w — v — v4) 


8.11 


BU" ү? cos (8w— 2v,) 


8,0.2 


Bio 1 COS V 


В! 1 cos v, 


В, oo COS 4 w 


B 


4 


В 


4 


В? үпү! cos (4w—-v — v,) 


4.1.1 
BĘ) mn! cos (4w — v + v4) 
BC? n° cos (4w — 2v) 


4.2.0 


зо № COS 4 w 


оз 12 COS 4 w 


ВС) 4/2 cos (4w — 2v,) 


4.0.2 


ВС? qu cos (4w —v — vi) 


B! 1 cos (8w — v) 


8.1.0 


ВС? g cos (8w — v,) 


8.0.1 


ВӘ q? cos (1210 — 2v) 


12.2,0 


po 


1244 À 


BĘ), чү? cos (12w — 2v,) 


12.0.2 


m! cos (120 —v —v,) 


39 


Analogiczne wyrażenia znaleźlibyśmy dla spółczynników A. 
Zajmijmy się teraz wyszukaniem wyrazów o długim i o krótkim 


okresach wśród wyrazów drugiego rzędu. W tym celu, jak to już 


zrobiliśmy poprzednio, wielkość R rozbijmy na dwie części, dłu- 
gookresową Ва i krótkookresową А, tak ze R=R,+-R,. 


Ale R jest to szereg trygonometryczny, musimy więc rozwi- 


nąć oddzielnie część wyrazów drugiego rzędu, zależną od długo- 
okresowych argumentów wchodzących w R, t.j. zawartych w R; 
i część w której występuje Rr. 


Postać tedy tego rozwinięcia, 


w którym zawarte będą wyrazy drugiego rzędu, będzie następująca: 


Е, = (В +- Ra) | Bi), cos nw- Вг n cos (nw tv) +... 
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Wyrażenie w nawiasie }...| oczywista zawiera wyrazy o krót- 
kim i długim okresie, a także wyrazy zwykłe. Jak to powiedzie- 
liśmy wielkość R jest to szereg trygonometryczny, w którym wy- 
stępować będą argumenty kątów, analogiczne do tych, które są 
w wyrażeniu w) |. Przez wykonanie mnożenia wskazanego 
powstaną nowe argumenty, będące sumami lub różnicami po- 
szczególnych pierwotnych argumentów. Buchholz podał 
praktyczne wskazówki, kiedy w wyniku ostatecznym, przy mno- 
żeniu dwóch funkcyi różnych (krótko lub długookresowych) otrzy- 
mamy wyrazy o długim lub krótkim okresie lub wreszcie wyrazy 
zwykłe. Oznaczymy według Buchholz'a przez у i $ funkcye 
długo lub krótkookresowe, występujące w wyrazach szeregu, które 
mnożymy przez siebie i niech rezultat otrzymany będzie W, wtedy 
mieć będziemy: 


fr da Vas ga. G= Vi; Xe. = Va ya. =, 
l Хе. weht; X: - ‚=. Ч W,. 


Zastosujmy te wzory Buchholz'a do naszego przypadku 
i wyszukajmy wyrazy drugiego rzedu w rozwinieciach (45). 

Przedewszystkiem rozważajmy tę część iloczynu R,,, w któ- 
rym wchodzi jako czynnik wielkość Ва. W tej części wynajdziemy 
następujące wyrazy II rzędu: 


(2) 


Długookresowe Krótkookresowe 
Boss ph. VT UT 
DES y cos (4 w — v) guo 9 eot y 
BiU y cos (wv) BO cos v 


(3) Бы Y cos v 


Bu * cos v, 


В 11? n cos (81 — v) 


8.1.0 


Bos" 4! cos (810 —v,) 


Dalej w iloczynie R,, musimy uwzglednié takze wyrazy kry- 
tyczne powstałe wskutek mnożenia przez R,. Wyrazy długo 
i krótkookresowe będą przedewszystkiem wszystkie te, które już 
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wyznaczylismy dla poprzedniej czesci (w której jako czynnik wy- 
stępowało Ra); oprócz tych wyrazów wystąpią jeszcze wyrazy po- 
chodne od wyrazów zwykłych, szeregu | | dla których argument 
jest równy 2v. Te wyrazy musimy uwzględnić, gdyż po pomno- 
żeniu przez Ry one dadzą nam nowe, wśród których będą wy- 
razy o krótkim okresie. Oprócz więc powyżej już wypisanych 
wyrazów (3) w tej drugiej części rozwinięcia dla R,, wystąpią je- 
szcze wyrazy zależne od następujących zwyczajnych wyrazów, mia- 
nowicie: 
Вуд» COS 8w 


BL? y cos (4w +) v) 
В} ү” тү cos (4w+- v) 
Bi; 7 cos (12w— v) 


Вб 1/ cos (1210 — v!) 
Wyrazy drugiego rzedu jeszcze powstana z iloczynu 
p K iEn Banoo sin nw +...{. 

Rozlózmy znów K na części składowe K = Ka K,; na za- 
sadzie poprzedniego wiemy, ze K;— 0, a więc pozostanie rozpa- 
trzeć tylko część 

LK | Ул Вьоо Sin w+...!. 


W wyrażeniu w | | następujące wyrazy będą te, które na- 
leży uwzględnić: 


Wyrazy długookresowe: Wyrazy o krótkim okresie: 
| 4y BI) q sin (410 — v) 4y B,,, sin 4w 
4p B, 5, лү sin (4w — v) 8u BU y sin (8w— v) 


8u BG) 4! sin (8w — vi) 
Wyrazy zwykłe: 
7 Зв В, оо sin 8w 
4p В? q sin (4w + v) 
4p Bin, т sin (4w + vy) 


12 В" q sin (12w—v) 


12.1.0 


| 12 В) 4 sin (12w—v) 
3* 
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W ten sposób wyczerpalismy juz wszystkie wyrazy w rozwi- 
nieciu funkcyi Р az do drugiego rzędu włącznie. Wyrazy trzeciego 
rzedu zwykle sa juz tak male, Ze mozna ich nie uwzgledniaé — je- 
dnak w przypadku ruchu niektórych planet krytycznych bliskich 
Jowisza będziemy musieli uwzględnić jeszcze wyższe wyrazy 
(trzeciego rzędu) w tych rozwinięciach, a to dla tego, że one nie 
będą o wiele mniejsze od wyrazów niższych rzędów, które uwzglę- 
dniamy. Z pośród wyrazów trzeciego rzędu należy uwzględnić 
oczywista tylko te, które są zerowego stopnia. 


š | B;4, cos 4w; 4 Bj, , sin 4w; 8 B, cos 4w 
(9) 
| Bi cos 12w; 1227, sin 1240; 62 B,,,, cos 12w 


te wyrazy przemnożone przez Ri К dadzą nam wyrazy o długim 
io krótkim okresie. 

W rezultacie będziemy mogli ostatecznie napisać następujące 
rozwinięcie dla Р. W tym rozwinięciu zostały uwzględnione 
wszystkie wyrazy szczególne aż do lll-go rzędu, zachodzące 
w przypadku ruchu małej planety typu Thule — t. j. takiej, której 


ruch dzienny średni jest bliski ы ruchu $redniego dziennego Jo- 


wisza. 


P= Booo + B.oocos 4w + В, 1 соз у + Ву, w cos vi + 


+В) 4 cos (4w — v) -- В у 11 cos (4w—v,) + 
+ BG) q cos (810 —v) + ВС т cos (w — v,)+ 

+ Boo т? + Boga 7? + В mm! cos (v—v,)+ B,,, n? соз 4ш + 
+В, о, m? cos 4w + ВС) 1? cos (4 — 2у)+ 


+ ВСЭ чү? cos (4w—2 v,) + ВЕ!) 4 cos (8w—v) + 


4.0.2 


FRE? " cos (8w — v,)+ pte? nm cos (4w+v— v,) + 


B.0.1 44.1 


+ ВС) m cos (40 — v - vi) + ВС) mn! cos (4 — v — v) + 


+ BC? qx! cos (8:0 —v — v) + ВСЭ my! cos (1220 —v —v,) + 
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+ ВС) 1? cos (8 w — 2v) + ВЕ n°? cos (8w — 2%;) + 
+В» 11? cos (12 — 2v,) + ВС» т? cos (12w — 2v)+ 

+4 Bzy, cos 8w+ В} n cos (4w + v) + Bi cos (4104 vi) + 
+В n cos (1220 — v) + Bzz; W cos (12w—v,) } Re+ 

+ { By + Bin cos 4w + Bras «cog v+ В үзү cos v + 


1.0.1 
+B n cos v+ Ву W cos у, + Ви 1 cos (80 — v)+ 
+Bqa cos (82 — v.) } (В+ R)+ 


+ (4B,,, sin 4w + 8B,,, sin 8w + 4B 7) q sin (45 —v) + 


4.1.0 


+ ABC g sin (4w — vj)--ABG y sin (4w + v) + 


11.0 


+ ABUS лу sin (4w + vı) +8 BG) q sin (8ш — v)+ 


8.1.0 


+8BGY g sin (81 — v,)+12 В), y sin (120 — v) + 


12.1.0 


+1280 v sin (12:0 — v)} p (K,+ K, +) 
cos 12w }R? + { 4В;" 


oo Sin 40 + 


+ (Bio, cos 4w+- В?? 


12.0.0 


+12Bl},, sin 12w } в К, R, — 4 8B,,, cos 40+ 


12.0.0 
+62B,,,, cos 12w u K; . 


12.0.0 


Analogiczne wyrazenie otrzymamy dla wielkosci Q, tylko od- 
powiednio zamieniając spółczynnik liczbowy B ze wskaźnikami na 
odpowiadający spółczynnik À. 

Otrzymamy tedy: 


@=А, oo Sin 40-4 1 sin v АС y sin (4w—v)+ 
HAUT nsin (8w— v) + АС sin (4W0—v,) +46) 4' sin(8w—v,)+- 
FA m! sin (v—v,)-+-4,„, q? sin 404445, 1? sin (4w — 2v)+ 
LAG) т? sin (8ш — 2v)+ 472, т? sin (12w — 2v) + 


HAED m sin (4w--v — у) AGA mn! sin (4 —v—v,)-- 
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+4, m sin (8w—v—v,) + AG), my sin (1210 —v — v) + 
+4) m sin (4 — v-4-vi) + AC? т? sin (4w — 2v) + 
Taba n° sin (8w— 2v)-- 469) 7° sin (12w — 2v) + 

+ (4%, sin 81--4717 s sin (4--v)--45, q sin (120—v)+- 
+ Aja" sin (4ш + v H- Az P^ sin (120 — v) ) Bi 

+À 44, Sin 404714 nsin (4w—v) + A; 17^ y sin (8t0—v)-I- 
+ Ama sin v, + Аш 7! Sin (dw—v,)+ 

+ Ал” W sin (8ш—у,) ) (Ra + Ri) — (44,0, cos 4w + 
-H84,,, COS 8w +4 А02) n cos (40-4v) +44} cos (4w-]-v,)4- 
+444,,, т cos (4w — v) + 8 АС! q cos (8 0 — v) + 

+ 12 45157 cos (120 — v)-I-4 AL, 7! cos (4100 — vi) + 

-+Е8А , W cos (8W—v,)--12 АС? у cos (12w—v,) ) p (Ki +K)+ 
+ 1459, Sin 4w+ А sin 120) R2— (4415, cos 4w + 

+ 124,50, cos 12 w jp Ro К, — 484,,, Sin 4w+- 

+ 624,00 sin 12w} p?K 244 A, oa cos 8104-47 cos (4104-v) + 
+ 451, 1 cos (1210 — v) + ACE Y] cos (4w4-v,) + 
+ Az, cos (120 — v) ). 


Abyśmy jednak mogli ostatecznie obliczyć wielkości P і (Y 
musimy mieć możność obliczenia 


(*) Sa, Sk, Ra, By, Ky. * 


Wielko$ci te sa to znów szeregi trygonometryczne, w których 
wystepowaé beda wyrazy o któtkich lub dtugich okresach. Musimy 
więc obecnie, mając rozwinięcie dla P i Q, wynaleźć jeszcze wartości 
dla wielkości (*), które podstawimy w nasze wzory. Przedewszy- 
stkiem należy zauważyć, że we wszystkich wielkościach (*) wystę- 
pować będą wyrazy krótko i długookresowe, zależne od tych sa- 
mych argumentów, które już występowały w rozwinięciach dla P 
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10. Mieć więc będziemy w rozwinieciach dla S, В, K wyrażenia 


nastepujace: 
Dtugookresowe. Krótkookresowe. 
Wyraz stały 
í т? n COS V 
y? n’ COS vi 
тїї! cos (у — v4) cos 4w 


(9) 


n cos (4w — v) 

" cos (4w — v4) 
12 cos (8ш — 2v) 
12 cos (810 — 2v,) 


y cos (810 — v) 
т! cos (816 — vi) 
12 cos 4w 
mn cos (4 w-Hv—v,) 
т! cos (4 w—v--v,) 


т! cos (8w —v—v,) y? cos 4w 


n? cos (4w— 2v) 
n’? cos (4w—2v,) 
nn cos (4w—v— Vi) 
a COS (12w—2v) 
1? cos (12w—2v,) 
ттр cos (12w—v—v,) 


W tych wyrazeniach (9) sa zawarte wszystkie kategorye argu- 
mentów Gyldén'owskich, o których wspominali$my na str. 15 
wzór (8) mianowicie argumenty typu A, B, Ci D. 

Z tego, co$my teraz powiedzieli widzimy, ze do typu À na- 
leżeć będą 4 pierwsze wyrazy długookresowe, będą one elementar- 
ne, gdyż te wyrazy są drugiego stopnia względem n — wszystkie 
zaś pozostałe długookresowe wyrazy należeć będą do wyrażeń 
charakterystycznych kształtu C. Podobnie, do wyrazów elemen- 
tarnych krótkookresowych należeć będą dwa pierwsze wyrazy 
o krótkim okresie, reszta zaś należeć będzie do wyrazów kształ- 
tu D 

Aby móc lepiej zbadać istotę wyrazów (9) postarajmy się 
wyrazić wszystkie argumenty w postaci tej, jaką Harzer otrzy- 
mał dla każdego typu charakterystycznego i który wypisaliśmy już 
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na str. 15 (8). W tym celu przypomnijmy sobie, ze mieliśmy 
związki następujące: 


3-58, 
m | 
skąd: 
4(1 —p)=1-+8; 4w=(1--8,)0—4B— 4uW*, 
oraz 
Vy +o—Il; ') vy =v П, 
czyli 


у=0(1— 4) — Ibo; v,—v(1— sı) —Ho,, 


(gdzie П» stała, zaś Il,, znamy z teoryi Jowisza, $ ruch apsyd). 


Wprowadźmy te wielkości do wzorów (9) wtedy otrzymamy: 
Const. 


n cos (1 —p )v— Пе] 


(| т (Bi 
n cos [(1 —9,) v—II,,] 


xm cos [cv + (IT, — П,)) 


y cos [ 2,0 — 4B — 4p Wa + Ш 
| Y cos | 8,0 —4B — WELE] 
n? cos [22,0 — 8B — 8. Wa 2II] 
| 1" cos [28,0 — 8B — 8r W 4-1- 2IT,] 
mn cos [22,0 — 8B — 84 W a + TE 4- TT, ] 


cos [(1+ 3,)v —4.B — 4» Ма 
(10) q cos [(1--28,)0—8B—8uW;--II] 
n! cos [(1--28,)v0 — 8 B— 8u Wa+- IL] 
12 cos [(1-- à,)v—4 B—4yW.] 
тт cos [(1- 2,) v —4 B —4&W;—II— Ш] 
тт’ cos [(1-- 2) v —4 B — 4uW; +II—IL] 
1 cos [1 + 3,)v —4 B— 4p Wa] 
y? cos [(1-+ 8,)0 +4 B-- 4p W 4 — 211] 
1/* cos [(1 — 2,) 04-4 B+4pW¿—211,] 
mn cos [(1 — 2,)v4-4 B 4-4 W, — I—II] 
r? cos [(1 -Е38„) v—12 B — 12 Wa + 211] 
y? cos [(1-1- 32) — 12 B— 12. Wa 2II, | 
Um! cos [(1 2-32, v — 128—194 WII] 


(C) 


(D) 


1) Brendel. L с., pg. 20. 
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W powyzszem zestawieniu wyrazów tylko 5 nie zawiera 
wielkości 8,, jak wiemy, związanej z ruchem pewnej określonej 
planetoidy; te wyrazy odnosić się więc będą do ruchu dowolnej 
planety!), podczas, gdy reszta wyrazów wypisanych w (10) musi 
być obliczona oddzielnie w każdym przypadku szczególnym, dla- 
tego też Gyldén nazwał te wyrazy wyrazami charakterystyczne- 
mi, jako zawierające w sobie cechy ruchu pewnej badanej plane- 
toidy. Zwróćmy się do zestawienia tych szczególnych postaci 
wyrazów podanych na str. 15, wtedy z łatwością w (10) będziemy 
mogli wydzielić wyrazy, odnoszące się do każdego rodzaju, co już 
zresztą uczyniliśmy, wydzielając klamrami każdą grupę argu- 
mentów. 

W wyrażeniach kształtu (C) i (D) występuje Wa, które jest 
bardzo mało zmienne z czasem, gdyż jest to długookresowa część 
wielkości W — redukcyi czasu. Obecność tej wielkości utrudnia 
znacznie całkowanie naszych równań różniczkowych. W wyraże- 
niach kształtu A wskutek całkowania wystąpi mały dzielnik ç któ- 
ry jest rzędu masy zakłócającej. W wyrażeniach kształtu Di C 
przez całkowanie w mianownikach wystąpi wielkość 8, która, acz- 
kolwiek mała, nie jest jednak rzędu masy zakłócającej, lecz zawie- 
ra masę z wykładnikiem mniejszym od jedności: wyrazy więc te- 
go kształtu stają się zerem wraz z masą. Zgodnie z tym, co już 
poprzednio mówiliśmy, musimy zauważyć, że gdy całkujemy rów- 
nanie dające S, to przez całkowanie wyrazy kształtu 4 i C powię- 
kszają się, podczas gdy przez całkowanie równania dającego p — 
są zwiększone tylko wyrazy typu Bi D. Wskutek tej uwagi skon- 


statujemy, że S, będzie rzędu m’, zaś S, rzędu 5 ; Ra, Rx, W, sa 


, 
rzedu + ‚Wa i К, sa tegoż samego rzędu. Ма tej samej zasa- 


dzie można z łatwością stwierdzić, że zwykłe wyrazy w S, В, К sa 
rzędu тазу zakłócającej ?). 

Po tych uwagach, będziemy mogli napisać następujące wy- 
rażenia 3) dla S (w S wyrażenia kształtu (A) będą conajmniej dru- 
giego stopnia). 


1) Brendel. 1. c. pg. 106. 

? Kramer J. Theorie der kl. Planeten von Hecuba-Typus, pg. 22. 

3) Ludendorff. Jupiterstórungen etc. kap. I u. H sparsim.; 
Kramer J., L с., pg. 43. 
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{ S= SaF S, — a, + a, cos 4w 4-а cos v 
+2, 4 cos (4w—v) +a, т cos (8w—v) 


+ ау COS V, 
as 1’ cos (4w —v,) +4, 1 cos (8w —v,) 
+ aqq? cos 40 -|- a,n? cos (4w — 2v) 
Haa n? cos (81w—2v) +a,, 12 cos (12w— 2v) 


qui + авт cos (4w--v—v,) {+ а cos(4w—v— v,) 
a, nn’ cos (8w—v—v,)+ а, mn cos (12)? -v v4) 
Hasni cos (4w—v+-v,) + ау? cos (4w - 2v,) 
a 11? cos (8w—2v,)+a,, 1/2 cos (12w 2у,) 
+a,n'? cos 4w + 
| T 7) 


w tem rozwinieciu jest: 
Sa=4, +2 cos (4W—v )J- an? cos (80 — 2v) 
ayy! cos (4w — vy) H ay; nn! cos (80 — v — v4) 
— an’? cos (Bw — 2v) 

za$ S, oznacza wszystkie wyrazy w S ksztaltu 4. 

Analogicznie mieé bedziemy rozwiniecie (2) dla R, zwazy- 
wszy, Ze wyrazy ksztaltu (B) за juz zawarte w (p): 
R=b,+Bß, cos 4w + Bar, cos (4w—-v) +B,1 cos (Bw — v) + 


+ Byn’ cos (4w — V,)-+-8, 1" cos (8W—V,) + 
+ 8.1? cos 4w - 8,7? cos(4to —2v)--8,,1? cos(8w - 2v) + 
--8,:5* cos (12w—2v) 

(12)4 +87! cos (404 v — v) 4-6, mn cos (40 — v — v) + 
+B, 1" cos (Bw —v—v,)H-f, тд] cos (1200 —v— v) 

f, 10 cos (410 —v-v,)+B1311? cos (410—2v,) + 
+B? cos (8w —2v,) +6, 11? cos (12w — 2v,) 


gdzie R, oznacza czesé R, w której sa ugrupowane wyrazy ksztal- 


tu (A). 


1) Kramer Ј,. 1. c., pg. 28. 
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Wyrazy diugookresowe w R czyli R; sa: 

Rı=b,+B;n cos (4w —v )+ pn? cos (8ш — 2v) 
+Ë, f cos (4w — vi) Hs mn! соз (80 — v — у) 
+Ë? cos (8w — 2v). 

Jak to juz poprzednio zaznaczylismy, po zcalkowaniu w tych 
wyrazach nie będzie występować mały dzielnik ©. 

Ponieważ K¿=0, przeto wielkość K=Ki+-K,. Niektóre wy- 
razy zwyczajne będziemy musieli uwzględnić, gdyż one w naszym 
rozwinięciu będą już nawet pierwszego rzędu. Mieć więc bę- 
dziemy : 


( K —1, sin 4w +4, 1 sin (8w—v) +7, sin (4w — v) 
+7, 1! sin (8 w—v,) 
+7 1? sin 4w +7, 4 sin (4w — 2v) 


+Y 12 sin (12w — 2v) 
+1, m sin (41! ]-v—vj) Ha sin (40 — v — v.) 
(13) 4-1 sin (120 — v — v|) 
+7, m’ sin (4dw—v-+v,) +7, 1?sin (40 — 2v, 
+“ 9 sin (12w — 2v,) 
+, 7° Sin 4w Yo n? sin (4ш + 2v) 
т: Y sin 8w 
+7, т sin (80 +-v—v;) 
gdzie jak widac: 
K,=çt sin (4w — v) 
-Ет sin (4w + 2v) 472,7? sin 8w +120 5 sin (Bw v — v). 
dw . [dw 
qo ene p 
Dla tej samej wielko$ci miec bedziemy rozwiniecie naste- 


pujace: 
[#)= тец cos (4w — v) + 1,41? cos (8w — 2v) 


Wreszcie część o długim okresie wyrażenia 


(14) + 1,7? cos (40—v;) + 115 mn! cos (80 —v—v,) 
| +11 2 cos (810—2v;) +), 


O ruchu planetoid. 4 
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gdzie znów przez (99, oznaczymy wyrazy ksztaltu (А). 


Należy zauważyć, ze spółczynniki a są rzędu masy zaklóca- 
jacej, zaś spółczynniki a if są rzędu 5 (masa zaktócajaca po- 
dzielona przez maty dzielnik). 

Z wyrazeñ (12) (13) (14) odrazu widzimy, które rodzaje wy- 
razów zachodzą w każdej z uważanych wielkości. Napisane szere- 
gi musimy wstawić w wyrażenia wyżej otrzymane dla P i Q i wy- 
konać mnożenie, grupując razem wyrazy zależne od jednego argu- 
mentu. Przy mnożeniu musimy uważać aby w rezultacie, zgodnie 
z zasadniczą myślą Gyldén'a uwzględniać tylko wyrazy o dłu- 
gim lub krótkim okresie, zaś opuszczać wszystkie zwykłe wyrazy. 
które później w otrzymanym rozwinięciu na P i Q sprowadzimy, 
przez dodanie sumarycznego symbolu. Jednakże, ponieważ w ilo- 
czynie zawierającym K, wystąpią, jako rezultat mnożenia wyrazów 
zwykłych w K przez wyrazy krótkookresowe w nawiasie, nietylko 
wyrazy zwykłe, które zwykle pomijamy, lecz także wyrazy krótko- 
okresowe, więc musieliśmy także te wyrazy zwykłe uwzględnić 
w rozwinięciu K. 

Wykonajmy więc wskazane podstawienie i mnożenie, korzy- 
stając z elementarnych związków goniometrycznych dwuch kątów 
À ip: 


sin À . sin y= — ; cos 603 cos (А — в) 


1 


(а) sin А. cos в = 9 sin à +u)+ ; sin (A —&) 


cos À . cos в, = ; cos (A Hp) + 1 cos (A — p). 
Otrzymujemy wzór nastepujacy: 
P= Bono + Воо cos 4w-]- В}! у cos v+ В cos (4w—v) + 
+ Ву cos (8w —v) + Biz, q' cos v, + Вул cos (4w —v,) + 
+ Bzy, 1 cos (80 —v,)-- В, Boog FB; cos (v— v) + 
4- В, osn? cos 4w-4-B,,, qq! cos (404 v — vi)+ 
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+ B 1? cos (4w—2v) + B,5, 1? cos (8v — 2v) + 

+ B5, cos (12:0 —2v) + В, ү, qq cos (4w — v + v.) + 

+ Ву? m cos (&v — v — v) + ВЕЙ m cos (80 —v — v;) + 
+ Biz, gm! cos (12w—v— у) + Bios 1? COS 4w + 

+ Bio. 1? cos (4w—2v,)2-B,5, 11? cos (8w — 2v,) + 

+ Bs, 1? cos (12w—2v,) + 


+É By, cos 8w + В} q cos (4w4-v) + В ncos(12w- v)+ 


12.1.0 
+ Bai" cos (4w+-v,) + В w cos (12w — vj) ? 
4 В, cos 4w+f, y cos (8w—v)+8, n? cos 4w+B,, 12 cos (4w—2v)+- 
+ f,, 1? cos (12w—2v)+6, y! cos (Sw—v,) +8, qr cos (4w--v— v) 
+ Bia if cos (40 — v — vj) + Bum cos (12w — v — v) + 
+ Bari" cos (40 — v + vi) + 61? cos (dw — 2v) + 
+ 81912 cos (12w — 2v,) + Bion? cos 4w $ 
+ À Bio + Biog COS 4w + Baia” т cos v + 
+ ВГ” 4 cos (40—v) + Ву” 4 cos (Bu — v) + 
+ BZ 4! cos VFB, cos (4шф—у,)-ЕВу cos (8w—v,); 
(0, + 6, cos 4w + 8,1 cos (4w—v) + B, 1 cos (8w—v) + 
+ 8, 4' cos (40 — v4) + Вт cos (8w — vi) + Baj? cos 4w -+ 
+ 6,172 cos (4w—2v)+-5,, 1? cos (810—2v) 4,7 1? cos (12!0—2v)4- 
+ Bs t! cos (4w--v— у) + Gi, mn cos (40 — v — v) + 
- Bis n! cos (810 — v — v4) + Bra | cos (126 — v — vp + 
+ Bora! cos (40 — v 4- vi) + Bus? cos (5w — 2v) + 
+ Ben’? cos (80— 2у;) + Ban? cos (12!w— 2v,) 4- лү? cos 4v } 


Lu (4B,,, sin 4w + 8B,,, sin & + 4 Biro sin (10+V)+ 
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+4B$, 1 sin (4w + v,) + 4B,,, n sin (40—v) + 
+8B;,, 4 sin (8w—v) + 12 B, , 4 sin (1210 — v) + 
+-4B,,, sin (4w ——у,)-Е8В ү sin(8w—v,) + 
+128, 1! sin (12w—v,)} 

} f sin 4w- y, 4 sin (8w — v) + т, y sin (40 +- v) + т, + 
+71 sin (80 — v,) + 7,4? sin 4w + Tn 47 sin (4w — 2v) + 
+ ттт? sin (12w — 2v) + ть т sin (4w + v — v) + 
Hram sin (40 —V — vi) + nann sin (1220 — v + v4) + 
Ето’ sin (40 — v 4- v4) + 1,7? sin (40 — 2v,) + 
+04 sin (1200 — 2v) + 1,47? sin 4w + 
Чот? Sin (4w-1- 2v) + a, sin 8w + ү, тр sin (Sw HV — v,)| 
+ Liga cos 4w + В, cos 12w| 1461? 4-3 B cos Bw 
+ (114819, sin 4w + 12 By,, sin 1210 | 

148,7, sin 8w} — p 18B,00 cos 4101-62 Bis00 cos 12w ; 

iy’ — 411? cos 8 |] №. 
Z powyiszego wyrazenia otrzymujemy ostatecznie: 


(15) Р=р, + pı cos 4w--pzq cos v +p,n cos (4ш —v )+ 
+ p n cos (Bie — v) + pa? cos (1210— 2v) 
T РТ cos v, +-p, cos (4w — v,) - p, т cos (80 — vj) + 
-- 2,7" cos (1220 — v — vi) 


+p: 7 +p 12 cos (40—2v)+ 
-Ери л cos (810—2v)-|-py, rf? cos (12w—2v,) 
+p, ч} cos 4w pis m! cos (40—vV—v,) + 


pis m! cos (Bw— v— v4) Pag лү! COS (V— V,) 


Pro mm COS (40-|-у—у,) + 11? cos (4w— 2v,)+ 
р, 11 cos (Bw — 2v) 
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pu mm! cos (4w—v+4-v,) 
HP 1? cos 4w 
+P т? +6, 


Przez G, oznaczyli$my sume wyrazów zwyklych. 
Dla krötkosci oznaczyliśmy : 
D=B, + $ Bino By + 28 Biooti + $ Booo tt 
Pi=B,00 + Boo + BI UB; ao ti t 
+4 Bist Вы Вало 
+u (Bros + 3Bi ls n He |— 2B,,,+% Bios! 
РЕВ +В B + 4 Вало Bi +4 Bio, Bu ñ + 
+28 Boote + 49 Bano 14 + 2v BE n + Że Br Tı 
Bros F4 В, F4 Bios Bi + 4 Big b + 4 Boi h + 
+44 Boots + 29 Bhoi + 28 Вт, 
Pi=B o + Booo 8,4-4 Bio tb Bo Bi + 4 Вю В 
+26 Boots + 46 В, оте + 4e BZ 
ь›=Вү-ЕВ Bs +4 Bino Bs tb Bigo Bs +3 Baar b +4 Ban fı 
-F2y Boots + Au Вот, 
Bao + Bio Bit Booo fat t Bino fa + Bao BiH ORBo т, 
— 28 Bo 
Pr= Bon +4 Bios Въ Вь +4 Boo b + 4 Bon h 
+6p Broi Ti — 24 Bios 
p=B,,, + $ Bigo fr +B, Bath Вы b + 2вВ т; 
+ 408, tn + 2 В + 49 Bas T 
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Р.В, + $ В, F4 Вто Boast +4 Bono Ps + 
+ ¿Bo Ba HA Bzyy BaH 2e В, та +48 Bund 
+ 6p Bro 14 + 28 Вт, 
Фо ВИЛ В, Po +4 Biron Bet Boo fe +4 Biso в + 
Bas bs + 4 Bu & + 2e В, ото + A Выль + 
+ 64 Въ + 2% Poe Ts 
ри==В Bzy, bat i В, ts + Bono bo + à Boon 8 + Bior Ba 
+4 Hn Bs 4uB, оо т + 6t ө; T; + 29 Da Ta 
Pa=B os + 4 Ba fio Ві. +4 Bios ñ +4 Bons В + 
+ Ве, Bst Boo Pio + 44 Bano то + бр Bizo Ts + 
+2 Bio, Ys 
P Boos + } Ba bo В" Ps Ва Ps +20 Bosse 
TAB 8,017, 
Ри B.) Вызовы + 1 Bio Bet Bono bu FA Bisou Bona fs 
+46В, „т. 49 Bao, to 24 Bh Ta + 4v Biot 
PB Вью b +) Bia 6, +3 Bor 8, + Booo Bia + 
+4 Bo fis + BE, Au Bago tie + 
+ Ap Bots +R, + Bao, Te 
Ри Bios Hi Bios hai Вол FF Boan Pist 4 Во + 
+4} Bon Bote Bago To + 29 Bio, Ts 
ри=В, + B br Воо 4 Boo Brt В b + 
THEY № — 21 В, ostii + 28 В, от: +64 Bite 
iB, 11 + Bn Pist 4 Воо Bio b Во Bs + 
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HEBR Ba YB. Bs + 4 Bros 8-28 Boats 
+28B, 0.0 Tie + "is Те 

Pis = Вуз + Водо + $ Bu Bi Воб à Bios в — 

— 15228, — 2&B, Та + 28 Boo Vis 

Pro = Bso tt Bano Bu + Bo biri Broo tt Buo B + 
Bio 66— Bano — 28 Bos tà 

ри = Bi $ Bra Pn + Bono ñu tt Bad +4 Brio B + 
+480193, + Bae Ba + 3 Bot Bot — 
— 2& Bio, 1, — 2&В 1» 

Pa = Во: $ Beno bis + Boao Piet $ BI B В Ps + 
TB B— B, T—2 Boots: 

Pas = В} YB Be + 4 Boo Bs + 4 810068, + Bu + 
+4 BG B +4 Baa BB oo AB tt Butt 
+ 29 Bio, te + 49 Boote Bao, Ta- 


Analogicznie postępując z Q otrzymamy wyrażenia następu- 
jące : 


Q= Asoo sin 4w + Аф, т Sin Vi+ Альт Sin (4w—v) + 

+45, y sin (Sw —v)+ А rin” sin (v—v )+ 

+ Ajo 1/ sin (40 —v,) + 475, W sin (8w —v,)+ 

+ А, от? sin 4w + 4,5, n° sin (40—2v) + А72 12 sin (8w —2v) 
+ Aus 9 sin (1210 2v) + АН, my sin (40++v—v,)+ 

+А T sin (40—v— vj) + AGE, md sin (8%—v—v,)+ 

+ Ag sin (120—v—v,) 3- 477, my sin (80—v+v,)+ 

+A 1? sin (410 — 2v) + Аз, 1? sin (8w — 2v) + 
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+ AG, 102 sin (120 —2v,) + 


+ {А10 Sin 8w+ A777" «sin (4 + у) Ai; у sin (120—v) + 
+ Agga Sin (120 — vi) }. 
j B, cos 4t» 4- т cos (Bie — v) 4-3, 7 cos (810 — v4) + 
+B; n? cos 4w + f, 1? cos (4w — 2v) + 44: 1? cos (12w — 2v) + 
4-8. т" cos (4w + v — v) +82 m cos (40 — v — v) + 
4-8, |! cos (12w — v — v4) + Bin cos (4w — v + v4) + 
Bis n° cos (4w—2v,) + 8,4»? cos (12w — 2v,) +- Bj, 1°? cos 4v} + 
+4 445, Sin 4w + 4,5" 7 sin (4w —v) + Az ут sin (Bie — v) + 
--ATH q' sin v, + Ауу sin (410 —v,)+ 
Ago W Sin (8w — v) } 
| bo 8,7 cos (4w—v) + B, 1? cos (8w—2v)+ Br cos (410—v)-1- 
+B, nt cos (Bw — v — v4) +B, 1? cos (8w — 2v) + 
+B m cos (4w — v — v) + fa cos (12w —v — у) + 
+ 8. yq. cos (4w — v + v) + B, 17? cos (40 — 2v,) + 
-+ 8,012 cos (1200 — 2v,) {8.172 cos 4w-- 8, cos 4w + 
+ Bin cos (80 —v) - л cos (80 — v,) +8, 1? cos 4w + 
+ Bain! cos (4w-1- v — vj) НВ» тт’ cos (410 — v —v,) + 
+11 cos (40 —2у){ — 
— {44,,, COS 4w + 8A,,, cos 80 + 44, 7 cos (Aw-tv) + 
+ 4AŻ), n cos (4wv,) --44,1, y cos (40 — v) + 
+ 84;,, 1 cos (Bin —v) +12 А), q cos (120 — v)+ 
+ 4470, 1 соз (4w — v) + 8 Aus, 11 cos (8w — v,) + 
+12 45,077 cos (1200 v,) bt. 
f y, sin dw т, т sin (810 — v) H- 3,7, sin (4w + v) + 
Hys +1, sin (8w — v) H- yrn? sin 4w + т, 12 sin (4w — 2v) + 
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Hin? sin (1280 — 2v) +7, 11! sin (40--у — v) + 

Fa qt Sin (410 — у —v4) Ри! Sin (120 — v + v) + 

+ лү! sin (4,0 — v+ v) +H nsn? sin (44 — 2v) + 

+77 sin (12w — 2v4) + 7,07"? sin 4w- tson? sin (4w+2v) + 

Han? sin 80 +13 77 sin (Bw - v — v;)| + 

+ < Af, sin 4w-I- A12, sin 12w) {1° +3 81, cos ёш ) — 

— {44100 cos 4w + 12412, cos 12w } uf 48,7, sin 8w ) — 
— e? 4 8 A,,, sin 404-62 A,,,, sin 12:0 } (41? — 41? cos 8w } 
--4,,, COS 8w + AŻ; т cos (4w--v) + А т cos (12w —v)+ 
АН, т cos (4w + v) АС ч cos (1900 — vi). 


12.0.1 
Te wyrażenia uporzadkujemy, uwzględniając związek gonio- 
metryczny (a) w ten sposób, abyśmy otrzymali wyrazy jednako- 
wych argumentów dla wyrazów charakterystycznych oraz elemen- 
tarnych. Zachować będziemy musieli także i niektóre wyrazy 
zwykłe, zależne od argumentu 2v, gdyż podobne wyrazy wystąpią 
póżniej w wyrażeniach dla W. 
Gdy wykonamy wskazane podstawienie oraz przekształcenie, 
otrzymamy wtedy ostatecznie wzór następujący: 
Q=4, sin 4w + 9. sin v+ q, sin (40 — v) 
+ gq | m sin (8w — v) 
+4, sin v, + 454' sin (40 — v) 
qr sin (Bw — v,) 
4-9 72 sin 4w + qun? sin (4w—2v)+-q4,4? sin (8w—2v) 
4-9,1? sin (12w—2v) 
(16): + qai! sin (41 - v — v,) +g sin (4w — v — v, + 
+ qu W! Sin (8w—v—v,)4-919 11! sin (1200 —v—v,) 
Haoni sin (410—v4-v,) + 944? sin (dw — 2v,) + 
ант? sin (8ш — 2v,) +95? sin (120 — 2v) 


9,117? sin 4w + Qu! sin (v — vy) 
g, sin 8w +9, 4 sin (40+v) +g,n sin (12w—v) 
| +91 sin (40+v,) + g; 4! sin (12w— vj) + Ga. 
O ruchu planetoid. 4* 
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Przez g oznaczylismy spólczynniki wyrazów charakterystycz- 
nych i elementarnych, za$ g oznacza spólczynniki najwiekszych 
wyrazów zwyktych. 


Wartości dla poszczególnych spółczynników są następujące : 
ФА, + 44500 В + 48А, Ао 400, Ano Bi” 
Aa An 1 11 1-62 | — 44400 24400 У A20011 

Ф = 1445" Bi + A 844-14 008 — RY ve BL — 264,007 + 
F4? Aano Te + 264 т —2 4,1, 1 | 
=A НЗ doi h + 445006 оо — #4 оу + 
+48 450075 + 28 Agoi — 2F Aoa 
U=A Аро 8 В — 26 doo, + 4F Азот, + 
+ 4®А лет, 
0s=4 70134400 — Años” B MA — 28 Asso + 


+ 4 Azot 
ЧА Ави в +4 44006 +4 Ao Bi + бы Ао a —- 
— 4p Aio 
Az --4 Аер, 4-441, B, + FAL B, FOR А; — 
— 2p Аа 1 


Ф.А, + $ Anio B+4 ino Br — 44510 Bath Ао Ë — 
— 2& А, otn + At Auot + бр, 14 — 2F А, o Ta 
qu = Al 44505" Ba +4 Ао An Bs + Аа 8 — 
—2 В, Tas + jów + 6p Az, 14 — 20441905 
du Aca FHA b + 4800 b EA Fa EA h + 


+44 B + 4p Loos + 6p Amots — 2y Ai Ya 
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(Qu A ug ti Ano Fs F-4 4008104001 вв — Aion Bs + 
+4 Ax "+ 4p Asoo Tio F 6p 3, T; — Że 47 ол Ts 
Q12—4 420— KAT BH Ао оби — Año Bu — 4р Asoo Tar + 
+48 400120 — 2p Aj ota T4 Arte 
ФА, mL a Ë. — 4” 5 — $ А8001 — 4 Ао В = 
— 4401 B— 4h Ао — 2 Ат Aol 4A0 to 
AQ t Alor Bs — 3 Are Pio — Азове — FA В — 
—4&A,,, Tio — 29 A Ys 
MA IA o $ А Br А Aal — 
—2ywA,,,1:- 6 A0 6 
=A LFA, b12- dA o0 бів A, b + à Ano Bt 
+ АД! 8,—2p. 4,55 Tia — ŻE A „ти + 60 4175, to 
qu=A à А b — 4 44.0 Pis bio ВН Аз Ps — 
— 2&A,oo Ta — 24 А, 19 
Arno + ABF 4 Atos Buc Ano В hdi Pa 
АА 0011 — AB, Ya 
Фо = А52 +4 Arco Ао fis Ano" Bs +4 Bios BaH 
44451198, + 4 Ago &— 464, ти — ŻE Ао Ти Альт 
в = Argos Hr Aso Bis + Arco bed AG, Bet 4450) Bs — 
— 4А 1з — 24у Ts 
Qu = Aja Ао P Hb Año — Ао Bs НА Ba 
nr. & +4 Aot? 8,— 2584,05 Ts 2 A oos — 
— Ay Anoo Taa — 2t Ато — 445196 + 48 Ат, 
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9145,03 А В, + Азоо — 24,00 Ti Ө A, ooi 
»=А FB Ato Ba А Bd soob АА ли T бА, от, 
а= ау Аа fh Ао А h Ass + 

+ 4р Años 'i T 99 A155 Ts 
0041.-440 +4 Agro" b +4 450,0 hA от – 


—1.10 © 


— 4р dż a+ 84,600 1e + 8845610 T 3 Авто Ё 
ga = Ало ti 4506 + 4 4go3 B E 40 — A Sos — 
—4y A т. + 8g A Ts: 4 Ao В. 


IV. Całkowanie równań różniczkowych, 
odnoszących się do ruchu planetoidy (279) Thule. 


W rozdziale 1, równ. IV, V, VI dają nam wyrażenia różnicz- 
kowe spółrzędnych Gyldćnowskich; w poprzednim rozdziale 
zdołaliśmy już oddzielić wyrazy typowe odnoszące się do ruchu pla- 
netoid. Brendel') i Kramer podzielili planetoidy wedle cech 
ich ruchu na pewne klasy. Według tego podziału planetoidę Thu- 
le (279) zalicza się do planetoid pierwszej klasy. W ruchu tej kla- 
sy planetoid wyrazy stopnia 0 w Ri K są znaczne w porównaniu 
z masą Jowisza. W tej klasie planetoid wystąpi jeden bardzo znacz- 
ny wyraz z w rozwinięciu funkcyi R, i ten wyraz, jak to wykazał 
Brendel*) będzie kształtu D, według oznaczeń Harzera czyli 
charakterystycznym krótkookresowym. Do tej samej 
klasy co Thule należy także i planetoida Hekuba, której ruch, jak 
wiemy, był niejednokrotnie badany przy pomocy metod Gyldć- 
nowskich?), lub ich modyfikacyi. Zasadnicza jednak różnica po- 
między temi dwiema planetoidami polega na tem, że Thule przybli- 
żając się znacznie do drogi Jowisza, z jednej strony ulega silnym 
perturbacyom, z drugiej zaś dla tej planetoidy wskutek znacznej 

1) Brendel Th. d. kl. PI. I, pg. 92; Kramer. Th. d, kl. PL, pg. 20, 
str. 199. 

*? Brendel. Om amvändningen, str. 26. 

3) Patrz prace Harzera, Kramera, Backlund'a. 
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blizkosci drogi Jowisza i drogi Thule rozwiniecia wedle poteg sto- 
sunku osi lub promieni wodzacych sa bardzo powoli zbiezne. Aby 
otrzymać wystarczająco zbieżne rozwinięcia perturbacyi Thule 
(tak, aby położenie planetoidy było określone z dokładnością oko- 
ło I’ na przeciąg jednego stulecia) możemy w rozwinieciach spół- 
rzędnych i w równ. IV, V, IV wyżej wymienionego rozdziału — po- 
minąć wyrażenia rzędu trzeciej potęgi masy Jowisza, zachować 
zaś w wyrażeniach 0 stopnia wielkości rzędu trzeciej potęgi masy 
Jowisza, podzielonej przez mały czynnik ç. Natomiast w wyraże- 
niach I-go i I-go stopnia wystarczy uwzględnić wyrażenia rzędu 
kwadratu masy zakłócającej, podzielonej przez mały czynnik c, 
a to dlatego, że w tych wyrazach występuje jeszcze bardzo mały 
czynnik y lub y! w pierwszej i drugiej potędze. W rozdziale УП 
pierwszej części swojej „Theorie der kleinen Planeten* Brendel 
5 š š LUD. 

podał równania różniczkowe, wyrażające duż do 

Ponieważ w zadaniu naszem postanowiliśmy się ograniczyć 
do wielkości rzędu do trzeciej potęgi masy ciała zakłócającego, 
dlatego też musimy rozwinięcie Brendel'a jeszcze dopełnić 
w sposób analogiczny, jak to zrobił Kramer dla Hekuby!). Gdy 
dodamy do wzorów Brendel'a, odnoszących się do ruchów He- 
kuby, wyrazy dopełniające, otrzymamy wtedy, z wystarczającą zu- 
pełnie dla naszego celu dokładnością, wzory różniczkowe na związ- 
ki pomiędzy spółrzędnemi Gyldćnowskiemi dla planetoidy 
I klasy. 

Ogólnie każdą funkcyę zależną od v to znaczy Q, Р, R, S, 
К, W it. p. napiszemy w ten sposób, że wydzielimy w oddzielną 
grupę wyrazy 0, 1, 2 it.d. stopnia i każdą grupę opatrzymy 
wskaźnikiem, który oznaczać nam będzie stopień uważanego wy- 


razu, np. 
= €; FQ... 


znaczy, że funkcya () składa się z wyrazów stopnia 0, 1, 2i t d. 
Na stopień i rząd dwu wyrazów musimy zwracać baczną uwagę 
przy mnożeniu odnośnych wielkości, aby w wartościach P, Q, R... 
nie otrzymać błędu większego od n Zajmijmy sie wyznacze- 
niem tych wyrazów dopetniajacych. 


1) Kramer. Th. d. kl. Planeten, I, pg. 44, 46 sqq. sparsim; Bren- 
del, 1. c. pg. 139. 
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Mielismy 
p=(p) +B, 


gdzie (p) zawiera w sobie wszystkie wyrazy elementarne kształtu B, 
zaś w À są zgrupowane wszystkie wyrazy elementarne kształtu А, 
oraz charakterystyczne kształtu C i D. a otrzymamy: 


de _ 20 
dv dv + = 


dR А š 
= możemy napisać w postaci 
ав (dR dR dR 
dv la 
gdzie wskaźniki O, 1, 2 oznaczają stopień odnośnych części w roz- 
łożeniu LA 
dv 
Analogicznie będziemy mogli napisać: 
Q=Q + Q, TQ i 8=8,4+ 8, TS. 
W AE S R wystepuja spólczynniki B, które, jak wiemy, 
sa rzędu 7 * ` (pg. 50) а więc wyrazy zależne od (Q, + Q, + Q») Fa 
wejdą w nasze rozwinięcie. Wskutek różniczkowania wyrażeń 


długookresowych R, i R,, gdy tworzyć będziemy wielkości = 
1 


i d ‚ spółczynniki В w tych wyrazach będą pomnożone przez mały 
2 


czynnik à,, przez co te wyrażenia staną się rzędu istotnie masy 
zakłócającej. Gdy zaś otrzymane wyrażenia zostaną pomnożone 
przez funkcyę Q, która jest rzędu masy zakłócającej, to w rezulta- 
cie otrzymamy wielkości rzędu kwadratu masy zakłócającej, które 
będziemy musieli jeszcze uwzględnić. Zatem mieć będziemy: 


ro [ie 


gdzie wskaźnik p oznacza, że w > uwzględnimy tylko niektóre 


wyrazy, mianowicie te, które są kształtu C. Podobnie przybędzie 
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jeszcze i wyraz — Q, (55. Poniewaz S jest rzedu m', to iloczy- 
ny nastepujace beda musialy byé uwzglednione: 


m. 


—28,0, (д); — 28,Po i —28,P,; —28,8, 9, 


+28/40Q, т sin v, + (8,4 ^y, 


gdzie przez wskaźnik (A, C) zaznaczam, że tylko tego kształtu wy- 
razy w odnośnem rozwinięciu są uwzględnione. Ponieważ dalej 


dn? . ! 
A Jest stopnia 2-go ale rzędu masy zakłócającej, przeto będziemy 
2 
musieli uwzględnić jeszcze wyrazy — 2 E (E bedace jeszcze 
0 


drugiego stopnia. W ten sposób w wyrazeniu rózniczkowem da- 
jacem p jako funkcye v zostaną uwzględnione wszystkie wyrazy, 
aż do rzędu trzeciej potęgi masy zakłócającej włącznie. 


Gdy do wyrazów podanych przez Brendel'a (!. c. kap. VII) 
dodamy jeszcze wyżej wymienione wyrazy (porówn. Kramer I, 
pg. 42 et sqq.), to otrzymamy poszukiwane rozwinięcie postaci 
następującej: 


(d?p ha dR ав 
dż HP =2 8 — Р, 95] — 28%, (5) + 
+28, — P,— 284 0 P, +-2 8, 8/40 + 
+01 sin v— Q (7 i). — Q, dob. + 
+28, — P,—28,40 Р, —28,40 P, + 
+28, (8,40 + 5,40 * - Q, sin v - 
- 4 do), — 0, 0), e (6 i). + 
dż 2. 


= ———————— ———— 


(40) i ر‎ 
-F28,49 Q,n sin v—2 de (do 
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Ograniczajac się do takiego samego stopnia dokładności, 
otrzymamiy analogicznie drugie równanie następujące: !) 


dS 2 
(2) ЗЕ" Q,— Di — 9,— 38,00 Q,— 38,0 Q, — 38,00 Q,— 9 dw 


Röwn. (1) i (2) są właśnie równaniami szukanemi. Zajmie- 
my sie obecnie calkowaniem. Zanim jednak to zrobimy, musimy 
równ. (1) i (2) jeszcze cokolwiek przekształcić. Wielkości Pi Q 
składają się z szeregów, w których występują spółczynniki nie- 
oznaczone. Np. w R mieliśmy różne spółczynniki B. w drugich 
zaś funkcyach spólczynnikl a i y. Aby sobie ułatwić zadanie, po- 
staramy zgodnie z myślą przewodnią Brendel'a?) wyrazić wszy- 
stkie spólczynniki 2 i y przez spółczynniki В; w ten sposób uzy- 
skamy wielkie uproszczenie naszych wzorów. Ponieważ w tej 
pracy ograniczamy się do rozwinięć liczbowych nie idących dalej, 
niż drugi stopień (Thule posiada bardzo mały mimośród — przeto 
w pierwszym przybliżeniu; pomijając wyrazy wyższych rzędów, 
niż pierwszy względem mimośrodu, nie popełniamy wielkiej nie- 
dokładności), możemy napisać następujące przybliżone wyrażenia, 
które zawierać będą wyrazy do I-go stopnia włącznie, mianowicie: 


(3) S=a,n cos (40 — v) -- a, 1 cos (4w — v,) 
(w tej funkcyi S są znaczne wyrazy typu С) 
oraz 
(4) R=Bß, cos 40-2, т cos (4w — у) B, «' cos (4 w — v, + 
+f, y cos (810 — v) +B; 1 cos (8w—v,). 
(Widzimy, ze w E wystepuja wyrazy charakterystyczne typu 
Ci D). Alez (1) uwzgledniajac definicye dla В, otrzymamy: 
d: R 
a ponieważ 
ФЕ 
dv? 
1) Ludendorff, 1. c., Kap. H sparsim; Kramer, l. c., pg. 43, 51. 


2) Brendel. Th. d. kl. Planeten 1. p. 134; Ludendorff, 1. c., pg. 8; 
Kramer id. p. 24; Brendel: Om amvändningen, pg. 29. 
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jest co najwyżej rzędu iloczynu тазу zakłócającej przez mały czyn- 
nik, więc możemy je opuścić, wtedy otrzymamy 

R98. 
czyli 
(5) R =2a,7 cos (4w — v) +24, q' cos (4w — v). 
Uwzględniając (4) i (5) jednocześnie otrzymamy z łatwością !) 


az Baj Ь=-у-. 


Przystąpmy teraz do wyznaczenia spółczynników 


їп (п= 1, 2, 3, 4, nhe 6). 
PAR š Я aw 
Mieli$my na str. 20 wyrazenie VI, dajace nam A" fun- 
kcyi R i S mianowicie: 
ай 
(6) dy =—2В- 8, =2 R, +6 Ry cos v+... 


Z drugiej za$ strony mieli$my na str. 37: 
W =< + У. К, 
gdzie yv oznacza wyrazy wiekowe. 
Po zrózniczkowaniu otrzymujemy: 


dw dw dK 
(a) rate ` 


ale 
(b)  K=r, sin 40 + 1,1, sin (80 — v) + тт) sin (4w + v) 
+, sin (8w— v), 
gdzie 
1=% +1; 
oraz czesé dlugookresowa 


(С) 9,1 cos (40 — v) + sn! cos (4w—v) —q,... 


dv 


1) Brendel. Th. d. kl. Planeten, p. 134. 


O ruchu planetoid. 5 
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Wyznaczymy te spółczynniki wyrazów z dokładnością więk- 
szą, niż a i B, mianowicie będziemy musieli pominąć dopiero 
wyrażenia 3-go rzędu, a z wyrażeń 2-go rzędu te, które są istot- 
nie rzędu pierwszego lub wyższego, przeto wyrażenie (6) może- 
my także napisać w postaci następującej: 

aW- : 

(7) de 2h cos 4w + (a; — 28, + 38,8, cos (4w—v)-- 
+ (х, — 28, + 38,8)" cos (4w— у) — 28,1 cos (8w— ү) -- 
—2,1 cos (8w—v,) + 3f44 cos (4w + v). 

Z drugiej zaś strony uwzględniając (a), (b) i (c) oraz (10) 
otrzymamy: 


W 
(8) Auer + (1 +3,)7, cos 40 4- 1,7, cos (dw—v) + 


+ (1++2д,)т,л cos (810 -v тул] cos (4w—v,) + 
+ (1+28,)5;4 cos (8w - v - (2--2,) 1, cos (4w- v). 


Z (7) i (8) wynika, ze: 
dw 
do |=) cos 4w- ya cos (4w—v) + 
+ (1--28,)744 cos (Bw — v) 4-7 1! cos (4w — ү) + 
+(1+23,)7; 1! cos (80 — v) 4- (24- 2,) ть cos (4w + v) = 
(9) ==— 28, cos 40 + (а, — 28, + 38, 8,) 7 cos (440 — v) + 
(ay – 25, - 36.6.) cos (410 — v) — 2, q cos (8w — v) — 
— 26, 1 cos (Bw — v) +36, y cos (4w + v). 


Z (9) z łatwością otrzymujemy związki następujące !): 
5 


928 
(10) пу 1=36—16+66; 
2 
ts = — 3/5 в + 38185 14 7— TS ; 
1 
28, 38, 


ls 13- 23, = 24, ° 


1) Brendel I, pg. 142; Ludendorff, 1. с. $ 4. 
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Gdybyśmy w wyrażeniu 1, ograniczyli się do takiego samego 
stopnia dokładności, z jaką wyznaczylismy 2,, to otrzymalibyś- 
my wzory następujące : 


(11) u=—2Bh; 1.38, —*/28Ba; 
neh; == — 261; 15=-—28s; 
1 = ?/5Й\. 


Na str. 57 mieliśmy wzór dla Q, w którym występowały 
spółczynniki q,, q, ...; napiszmy ten wzór, ograniczając się do 
wyrazów pierwszego rzędu tylko, mieć wtedy będziemy: 


(12) ()=q, sin 4w-1- 9,7 sin vgn sin (4w — v) + 
+ qua sin (80 —v)--Qs 17 sin v, +q, sin (40 — v) + 
+q, sin (8w — v,). 


Ale na zasadzie związków, wypisanych na str. 58, 59, widzi- 
my, Ze wszystkie spólczynniki q sa takze funkcyami spólczynni- 
ków ç i 8. 

Zamieńmy w odnośnych spółczynnikach wszystkie y i B, po- 
sługując się związkami (10) i (11) i uporządkujmy rezultat tak, aby 
spółczynnik В występował jako czynnik, wtedy otrzymamy związki 
następujące: 


qa O FO BE 
a = 0 +92") Bi + Qu + 929 Ba 
аз==@'% + 90 24-93% Bs + Ф309 6, 
Q= Q) + Qu в 
qs — 59 + 9,0 B; + q; Bs 
6 = Q4? -- 060 B, + q 2 B; 
q. — 4,9 4-99 Bi 027 63 


Analogicznie mieliśmy także w owem wyrażeniu dla (Y część 
wyrazów zwykłych w postaci następującej: 


@„== g, sin 8w+ 9.1 sin (4w + v) + g, y sin (12w — v) + 
+g, 7 sin (4w 4- v,) 4- g, т’ sin (12w — v4) 


http://rcin.org.pl 


68 J. KRASSOWSKI. 
bedziemy wiec mogli napisaé dla spólczynników g wyrazenia na- 
stepujace: 
Wn ha 
Y2=9 929-930) Ë. +99, +9:% Ë, 
(14) Js = GO وو‎ 8, +g 8, + وو‎ В, 
U= 949 +g b +g 28, +90 Ba 
Js =9; +- g, 8, + 9,9 Bs +90 Bs 


Z (13 i (14) oraz ze zwiazków na str. 58 í59 widzimy, ze 
spólczynniki q” są funkcyami wielkości Аі р, a więc są to wiel- 
košci liczbowe, które bedziemy mogli bez trudnosci obliczyé. 


Zatem mieé bedziemy: 
(15) Q= А; Q0 —4 A — ВВ Asao 


=8 4oo HA 00 Hi Arno + e (оу ,)— 
— 4h? 4,00 — 84,59 62 4,300 
1°=0; ЧА "+ A, a= — (3А, + SABIA 44,, io 
OA: 4.9 — A — Bu A, 
qi = Аз; uoi тти EE 4,1 
939 =} As 9 =} тра Su 4,00 
949 = A оз 09 =} а — 2% (4451—3400); 


q) = + het + 44 4,50: 


9.® = a5 4 ep "e au PT 16y.455, ; 
4: © = — +4 400 + 4 Ayo; 
ATI: (Da 4140 ATH 4 34 —9 
Qe 8103 48 — Aio + Ao ) MBA Aa 
qi? =} E ii 


9:9 = Аза; @®==КА А01) — (BA A); 


(3) — . 
9:9 = ¿Aro 
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oraz: 
(16) g/9—4,,,; HPA + 41200) 44, ЗА, оо) 
Же: 9 = y 4%." —8y 4}, о; 9) = PAZ; 
929 = 4 Aron 128 42005 
9, — А}, + 4 (An + A) — Sa At; A 
g, 9 = A; — 1244 
BA 000 Ао" + Аш) FM (Asie 


+3 4600 — 84511): 
9 =}. e gi e — (Ao FE 4,99; 


09 = Asus Som + “у + 
T 86 (459, — MU 


12.0.0 ? 


g, 9 =4(4 ni + 84 4,50): 


Zupelnie analogicznie postepujac bedziemy mogli napisac 
funkcye P w postaci nastepujacej: 


P= p, + p, cos 4w-|- p, 1 cos v + рат cos (4w — v) + 
+ ps n cos (8w + v) -- py 19! cos v, + p," cos (4w — vy) + 
+ p." cos (Bw — v) 
gdzie spólczynniki p maja wartosci nastepujace: 
P= py POR 
руриро 
Pa=Pp, OP: By P> Ra Pa BA 
PZP: h HP: NH", 
p =p) HN RD, 
DE, 
D,—p. E le CEE f, 
Par PR; 


(17) 


em — oD СС 
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(48) p,” = Bos; BE — 4 Во; PO =4 Bo; 
pO =B o P Bico ¿Bios — Sy Bao: 
SPE" BM В) (m5 389. yk 
+242 (31 Biroo —4B,00) 
p= БВ; P=} Be + Bas) Bas 
+485 — ЗВ...) 
Da? = 4 Вр; Pa” =} Beo, ВВ, о: 
в": p; =4 (Bis + lie - Au (B + 
+ Bons РУ Bio р,” = Ву Su Boop: 
p = Bg: рио = HB," + Bus )—2W4B,, 


ЗВ, 0); p = Воо 4р В 


4.0.0 ? 


810 


pi e Boss io Bin + Вы) 8w В m 


p= + BB P5 = В, — Boni 

p, = Bo; DO EN HB AGBS + B; Ж, 

ре = 4 Bs PO = + Bis; 

p =B РО (Bao + Bas) — 4һ(3 Bot 
HER) p= 4 B; RB: 


W ten sposób mamy wiec juz w granicach przyblizenia, któ- 
re zamierzaliśmy osiągnąć, wyznaczone wszystkie spólczynniki 
rozwinięć P i Q, jako funkcye pewnych arytmetycznych wielkości, 
które, jak to później zobaczymy, możemy zawsze bez wielkich 
trudności obliczyć. 

Przejdziemy teraz do właściwego całkowania naszych rów- 
nań różniczkowych zasadniczych. Całkowanie tych równań wy- 
konamy zgodnie z metodą wskazaną przez Brendel’a'), którą 


') Brendel. Th. d. kl. Pl. I kap. VII; II kap. Ш. 
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z wielkim powodzeniem zastosowali Ludendoríf, Buchholz 
i Kramer. W niniejszej pracy ograniczymy sie, przy calkowaniu 
naszych równań, do wielkości zerowego stopnia aż do rzędu trze- 
ciego włącznie. 


Zacznijmy od całkowania równania dającego nam e. gdy sie 


dv’ 
ograniczymy do wyrazöw rzedu O, wtedy ze wzoru (2) otrzymamy 
ds 
dv x coc Qo , 


ale na zasadzie wzoru (12) jest 


@ = 49, sin 4w, 


a zatem 

as H 
(19) ро ыр 49, Sin 410 
a ponieważ 


w—(1-—g)v — B—yuW, 


zawiera w sobie wyrazy okresowe rzedu wy2szego, niz zerowy 
(a także i wyraz wiekowy, jak wiemy), to po całkowaniu w S wy- 
stąpią dzięki tej zmienności argumentu w, wyrazy rzędu conaj- 
mniej pierwszego, ponieważ w S, niema wyrazów zerowego rzędu 
kształtu Ai C. 

Całkując (19) otrzymamy: 


(20) S= a, + a, cos 4 W + wyrazy zwykłe. 


Zastanówmy się bliżej nad otrzymanym rezultatem (20). 
Przedewszystkiem wyrazy zwykłe, które oznaczymy Soo otrzy- 
mamy z następującego wzoru: 


(= — ХА, оо Sin nw 


gdzie n=8, 12, ...; po scałkowaniu wyrazów zwykłych mieć bę- 
dziemy: 
(21) Sn 00 = — Ао — cos nw; 

a 0 T9) 


W wyrazeniu (20) wyraz a, jest to pewna stata!), która póz- 


) Brendel. A. N., 3346. 
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niej określimy. Aby wyznaczyć spólczynnik a, zwróćmy się do 
wzoru dającego nam Q w funkcyi spólczynników À oraz argumen- 
tu w (pg. 55) uwzględniając związki (3) pg. 37 (10) pg. 46 oraz pg. 
58 et sqq. przez różniczkowanie otrzymamy: 


4(1— p) Apoo H2 (1—1) А50, + 168 (1 — p.) 4,07 + 
+1—u) 45, + (1— Mol B — p 4 (1-м) Aja, — 
— 12 (1 — gy) Ае! Pi 714-821 —24(1— pa) А, oo H61 ва) А 2.0087 


po scalkowaniu otrzymamy: а, = 


A... + A0 + wnt] 1 Mawo Ku 
А1) 801-6) dm iT |16 (1-9) 1 
1 Азро | i "m 3Apo | 
O Tie: MAE ао 
p | 31 4,200 A, oo | E 
2 | 4 (—ь) i-a) 
ale 14-2 
1— p = + 1 
oraz Ecce 171 
u "dh 


Uwzględnijmy te związki i uporządkujmy (a) według potęg 
rosnących 8, wtedy otrzymamy: 


(22) a=] +a" &-r o" 8*4... 
gdzie dla krótko$ci napisalismy: 
Ás.oo 
терү 


9" = A 0.0 8.4, оо 
—30--8) (+8) 
(22%) AS 4% 4* 0 Эр. А". 0 


4.0.0 4.0.0 12.0.0 4.0.0 
PETER) ani E + e~ 
6p. 4500 16g? Arno 8824,00 , Ai 00 


(1-5) (+28): (143) (232) 
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Wyrażenie (20) po uwzględnieniu (22) przejdzie w nastę- 
pujące: 


(23) S,=a (o! 4-4 ' f, F9" 8°) cos 4w-1- 8,00. 


W wyrazeniu (23) mozemy juz wszystkie wyrazy obliczyé, 
z wyjątkiem tylko stałej całkowania a,, którą na końcu wyznaczy- 
my po scałkowaniu równań, odnoszących się do pozostałych spół- 
rzędnych Gyldćn'owskich. 

Przejdźmy do scałkowania równania różniczkowego, odno- 
szącego się do spółrzędnej R. Ponieważ obecnie chcemy pomi- 
nąć wielkości stopnia pierwszego, więc przedewszystkiem musimy 
zauważyć, że (р), = 0, gdyż w (р), jak wiemy), występują tylko 
wyrązy krótkookresowe kształtu B, których najniższych stopień jest 
pierwszy. Mieć więc będziemy ze wzoru (1) str. 63?): 

d*R dR 
(24) аз +В=28, — Po e] 
gdyż 
p=(p) +R. 

Analogicznie do tego, jak poprzednio pisaliśmy, wielkość Qo, 
możemy napisać : 

P,=p,+-Pp, cos 4w. 

S, znamy już ze wzoru (23), zatem w (24) pozostaje nam obli- 
czyć tylko wielkość 

ав : 
Qo (5 | ) 
otóż wzór (46) str. 34 daje nam: 


(25) Q= A, sin nw + В ХА, sinnw=pKtnA, po cos mv, 


za$ z (12) str. 48 znajdziemy : 
aR К йш 
du ) =—481 sin 4w E^ 


czyli uwzgledniajac (10) str. 46 znajdziemy : 


(26) 


1) Brendel, I. c., pg. 29. 
?) Brendel, І. c. pg. 90. 
3) Brendel, 1. с. pg. 123. 
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(27) (Е x = — &(1-1-à,) sin 4v. 
Pomnóżmy;(25) przez (26), otrzymamy: 
dB) _ : dh AT IE dR 
(28) Qo % |=" А моо Sin nw ( d br REA oo Sin nw | a um 


‚ (АВ 
—p Kin (To), A. COS . nw 
ale 
В, = — 8, cos 4w, 
a wiec 
(29) R, (5) — — В (1--9,) cos 4w .sin 4w— — 46,2 (14-0,) sin 8w 
0 

K, = 1, sin 4w 

mieé io bedziemy : 


(30) К, =), == — 1, 5, sin 4w.cos 4w=(cos 8w—1) aa + à, 
W (28) wyrazy krótkookresowe będą te, dla których wskaź- 


nik m przyjmować będzie wartość 4, 8 lub 12 a więc mieć bę- 
dziemy: 


(31) [o, %J | =— d E a cos 4w) 4 — 


u {А + P ENT Uh cos 4w + 


tel) An ЗА 1517, COS Aw. 
Wskutek tego (24) przejdzie w nastepujace: 


(32) Ge + Ada FID, + Be Dar + Du? + 


+ (W); В, + (Ve 1,2} cos 4w 
gdzie dla krötkosci napisalismy: 


(88) gra Baad e m 2uB 
P —— TE ih + 9 4 ooPB1— BD, 00 Ti — 


= Bi: 
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(34) (M=-= E 13, 
ике оге 
= — t Bus о 
@ = — 4 GS, — Ł (Bino + Bizoo) + 
+ Жыт ; + + наш, + Ama): 
(9), = — p (Bias + Baan ) — а кыны | 


— &(12-8)(4,,, — 34:69) 


2 
(Ds = 2р? Dios z 362 B a " Ano aja Aion 


1200 — 1-55 B 14-8, n 
Ale ze wzoru (12) str. 48 wynika): 
dR 
(35) div + R—b,—(1—(1-4-2,)*]8, cos 4w, 
czyli 
ФВ 5. -L 8.2) 8 4 
(36) dv? + R=b,— (22, + 2,?) B: cos 4w, 
gdzie dla krötkosci piszemy: ` 
b,— 2a, + po ; 
(a) 2% +8? — —1 (3), + (8), 6 + sa + (VE +0) Bir 
+ (Dhs 
zamiast 1, napiszmy: 
NE 
14-2," 


wtedy analogicznie do równania (23) str. 73 bedziemy mogli napi- 


) Brendel, l. c, pg. 120. Ludendorff, I. c, pg. 17. Kramer I, 
pg. 45. 


http://rcin.org.pl 


76 J. KRASSOWSKI. 


sac spółczynnik przy cos 4w, jako rozłożony według potęg B,, w po- 
staci następującej: 


PFP" f w AC 
przez porównanie z (a) wynika: 
(37) (20, +27)  — p' 3- p" ё. Fp" В, 
gdzie dla krótkosci napisalismy : 
2А 
1 — PRES, e 
(38) pr Boo 14 à, E 
1.0 
" 1.0 1 4 Aso, 14 0 7 Er 4 B5:0.0 
p = Bio ++ Bin. = 143,7 9 * Agoo "LE UE 
165 Agoo 
Tad 


p"=) Ha T Bios $ Bizoo А оо $4 00 — i Aż = 
= Ао, 3 Ао» —2 y Bis — бу. Biss — 4 A 
+ 124 A — 2 4,00 F бе 200 — SRB, + 
+ 728° B aoo + 488° Aygo — 1449? А... 
Uwzgledniajac (36) i (38) otrzymamy: 
J Q0 P+P'B р 6? 
(39) R= 2a44- p — D + 22 1 cos 410, 
gdzie znów dla krótko$ci napisalismy: 
4 D 1 д t 4 r, 
(40) Ро =— Воо; — 5 He Ë + $ A, 096 + 1 zt B,00 i — 
— Booofi*- 
Otrzymali$my wiec wszystkie wyrazy typowe funkcyi R--po- 
zostaje jeszcze obliczenie tych z wyrazów zwyklych, których 


w naszym przypadku pominaé nie mozemy. 
Możemy bowiem napisać: 


с +R=28, — Р, = Y (28,00 cos nw — В, оо cos nw), 


albo, gdy oznaczymy 28,00 — Воо przez Ч оо, to mieć będziemy: 
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2 
(41) e + В = W, оо COS nw. 


СаїКа równañ (41) bedzie: 

(42) R, = g, sin v — g, cos v. 1) 

gdzie g, i g, Sa pewne funkcye, które należy wyznaczyć. 
Z łatwością otrzymamy związki następujące : 


p. = } ХК, оо cos (nw + v) + $ EW, оо cos (nw — v 


Sh =}, оо sin (mw + v) — XW,oosin (nw— v). 
Scalkujmy te równania według metody podanej przez Gyl- 
dén'a Brendel'a 2). 
Jest: 


| cos nw + в) = sin (nw + v) + 


1 
SES] 
n y ` dwa 

+ 
forti = cos (nw + v)dv 


nu + v) dv == cos (nw + v) + 


1 

п(1— у) +1 

ne [йй a 
+ эй =й ‚| FIT sin (nt 3-v) dv 
a ponieważ Wa zawiera tylko wyrazy długookresowe kształtu А 
lub C, a wyrazy kształtu C mogą występować dopiero gdy uwzglę- 
dnimy wielkości pierwszego stopnia®), dla tego też w naszych wy- 
razeniach w całkach zależnych od Wa nie będzie żadnych wyrazów 
stopnia zerowego. Wstawiając ostatnie dwie równości w dwie po- 
przednie, znajdziemy: 


1) Kramer, 1. c. L, pg. 48. 
2) Brendel. Om Användningen $. 11, 
3) Brendel. Th. d. kl. Pl. I, pg. 120. 
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| A sin (nw -- v) + 
E + Реа |) Sin (mo — v) 
Ja = —4X d n cos (nw + v) + 
+12, а Yas cos (nw— v), 


gdy (43) wstawimy w (42), otrzymamy: 


R, = X R,oo cos nw, 
gdzie 
R TT ЖЩ 
EAS 
ale 


2 1223552 
Ҹ,оо== 28,00 — Bo i (1—p1)= ( de 1) 
wiec ostatecznie otrzymamy: 


28,00 — Doo 


(44) Run | g +š) y 
= = SST 


dla wartości п = 4 gdzie 5,0, jest określone przez (21). 


Pozostaje nam jeszcze całkowanie równania dającego zwią- 
zek pomiędzy W i v Zajmijmy się obecnie tem równaniem. 
Mieliśmy na str. 20 równ. VI, w którem obecnie uważajmy tylko 
wyrazy stopnia zerowego, mieć tedy będziemy: 


(45) “> =S — 2,R, — 28,8, + ЗЕ, + 38,2, — 48}. 


Wielkości S, i R, są nam już znane, jakoż uwzględniając 
związki na str. 48 mieć będziemy: 


So = a, + a, cos 4w; R,=b,-+-B, cos 4w, 


z których odrazu znajdziemy : 
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[S,R,=4,b,+a,b, cos 4w-]-a,8, cos 4w + а, В, cos? 4ш; 
R? =b? + 26,8, cos 4w + $,? cos? 4w; 
(46); R: =b,* + 3b,*8, cos 4w + 35,8,? cos? 4w--6,* cos? 4w; 
S,R2= a,b, + 2a, 6, В, cos 4w + a, 0? cos? 4w + 
| + a, b, гов 4w + 25, а, f, cos? 4w + 8,2 а, cost 4w. 


Zbadajmy wyrażenia (46). Wielkości: a,, bo, а,, b. sa 
istotnie rzędu masy zakłócającej m’, a więc a2, 0,2, a,b,, a b?, 
a,b, są istotnie rzędu т’? — należy je pominąć”); %, jest rzędu = 
gdzie à oznacza małą wielkość; zatem mieć będziemy: 

R? —2b,8, cos 4w-+-$,? cos? 4w 
ale poniewaz 


cos? 4w = cos 4w . cos 4w =4 cos Św + 1, 
to 


(47) В? = 2b, B, cos 4w 4-4 8,2 + B,* cos 8w 
SR, = аъ В, cos 4w + } а, 8, + 4 a, b, cos 8w 
SR? = $ a,f,? + aß? cos 8w +3 a,8,* cos 4w +} а, 8,2 cos 12w 
В = $ babi? + $5,8,* cos 8v + 1 6,5 cos 4w + 1 8? cos 12w. 
Gdy uwzględnimy powyższe związki, otrzymamy ostateczne 


równ. różniczkowe dające S. 


an © — e, +1 + (a, — 22,8 — 28 + 6bB + 


+} a,b? — 38,3) cos 4w — (a, b, — $ 8 
+ 4 a, fı? — 65, 812) cos 8w + (3 a, 5,2 — 8,2) cos 12 w, 


gdzie dla krótkości oznaczyli$my wielkości następujące : 


1) Brendel w przypadku Hekuby i Hildy pomijał wszystkie wyrazy 
rzędu wyższego, niż drugi, a z wyrazów drugiego rzędu odrzucał te, które są 
rzędu istotnie pierwszego lub wyższego. Ponieważ w moich rachunkach posu- 
wam dokładność o jeden krok dalej, przeto pomijam dopiero wyrazy rzędu 
wyższego, niż m’, oraz te, które są rzędu istotnie drugiego lub wyższego. 
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т= $ > 

Со = ho — 418, + š @ * — 60081? — 2b. 
W tej wielkosci с, wystepuje јако wspólny czynnik masa m', pod- 
czas gdy w wyrazie 3 6,2 występuje таза т’ podzielona przez ma- 
la wielkosé 2,. Wykonajmy calkowanie równania (47). Mielismy 
przedewszystkiem na str. 32 równanie (42): 


W= Tv + Wi; = K, 
różniczkujmy ten związek względem v, ROM 
aW — ай’. 
Ра, 


y oznacza część wiekową, a ponieważ długookresowe wyrazy 0 
stopnia nie istnieją?) w funkcyi W; ani też tembardziej w К, prze- 
to będziemy mogli napisać: 


“= 
(48) i42 x 
Porównajmy (47) i (48). Otóż widzimy, że 
(49) y= Co + Y , 


z drugiej zaš strony, uwzgledniajac zwiazek (13) str. 49, oraz (10) 
str. 46, przez rózniczkowanie znajdziemy: 
(49) = = (148) cos 40 + 2 (1 + 8,) K', cos 8w + 
+ 3 (1 + à) К", cos 12w, 
gdzie spólczynniki K', i K', oznaczaja wyrazy zwykle. 
Podstawmy w (48) związek (49) znajdziemy : 


ай 


(50) р 


= 1-1, (1 + 4) cos 4w + 2 (1 + 2,) К’, cos 8w + 
+3 (1 + 2,) К’, cos 12w. 


Porównajmy wzór (50) ze wzorem (47) to otrzymamy : 


1) Patrz wzór 14, str. 49; Brendel. Th. d. kl. PL I, p. 68. 
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1= Co +1 
کم‎ 24,8, — 28, + 60, а, 1° — 36" 
1+3, 
@) — _ 3 — 2a, B, —3a &* + 120,1? 
я 4 (1 + 3,) 
gi. За, 6,2 — 45? 
| "432 4-8) 


Do wyznaczenia stalych calkowania jeszcze powrócimy, tym- 
czasem ograniczymy sig do uwagi, że ponieważ uważana przez 
nas planetoida Thule należy do krytycznych planetoid pierwszej 
klasy, według podziału Brendel'a'), to dla takiej planetoidy 
mieć będziemy: 

€, =0, ale 12-0, 
a więc w naszym przypadku =r °) i równ. (48) przejdzie ро 
scałkowaniu w następujące: 
(52) W 170 + К. 
Przypuśćmy, że stałą a, juz wyznaczyli$my, wtedy uwzgled- 


niając związek 5, =2a,+-p'" na str. 70 i zw. (40) str. 76 będzie- 
my mogli napisać po uporządkowaniu rezultatu według potęg f,: 


b,—=24,+0Q0+ 4181 + q, 8, 


gdzie 
q = — Bo00; 
(53) y E 4 
nee — 4 Boot aM, g Аш РТ Ё. | Bus 
q | =s —3 n 


Zatem, gdy znamy stałą a,, to na zasadzie równ. (51) i (53) 
możemy z łatwością obliczyć już poszczególne części równania (52) 
i tym samym otrzymać funkcyę W, ograniczając się ciągle do wiel- 
kości stopnia zerowego. Postępując zupełnie analogicznie do te- 
go, cośmy zrobili gdyśmy szukali wyrazów zwykłych w funkcyi R 
(str. 73), znajdziemy z łatwością podobne wyrazy i dla W, miano- 
wicie będą one: 


1) Brendel, 1. с., pg. 92, 100. 
2) Brendel,l.c; Kramer I. $ 1. 
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Sn.0.0 PEE 2 Roo 
n 
(1+8) 


Przejdzmy obecnie do wyznaczenia stalych calkowania. Po- 
nieważ obecnie ograniczymy się tylko do badania ruchu planetoi- 
dy w samej płaszczyznie jej drogi, więc mieć będziemy tylko 4 sta- 
łe, które powinny występować w naszych równaniach po ich scał- 
kowaniu. Te stałe całkowania są, jak to wiemy, x moduł mimo- 
środowy, I' długość perihelium, a połowa osi drogi planetoidy 
oraz tak zwana stała ruchu n. Jak widzieliśmy po wykonaniu cał- 
kowania wyrażeń różniczkowych spółrzędnych Gyldén'owskich 
występuje jeszcze jedna nowa stała a,. którą, dzięki temu, że jest 
nadliczbową, możemy uważać jako stałą dowolną; tę stałą a, 
określmy za pomocą warunku, że b,, część stała, w wyrażeniu na R 
ma być rzędu istotnie pierwszego— wtedy także i b, będzie takiego 
samego rzędu (str. 70). Ta stała a, będzie tak dobrana, aby część 
stała c, w wyrażeniu (47) przyjmowała wartość zero. Ponieważ 
W—=W, +-K, а К składa się z funkcyi R i S, to stała część w W 
składać się będzie ze stałych części funkcyi R oraz S; wskutek te- 
go mieć będziemy 
(55) €,4-1— a, —29,. 


Ale Brendel!) wykazał, że o ile mamy do czynienia z plane- 
toidą charakterystyczną. to musi zachodzić związek następujący: 


a, — 2, =0 = e 1; 
dla krytycznych planetoid, a taką jest Thule, y nie może być ze- 
гет 2); a równość 


(54) W.00= dla n=4. 


a, —2b0.=C, +7 
pozostaje ważną i dla planetoid krytycznych, będziemy więc mieć 
związki następujące: 
(56) €,—0; 1=F0. 
Z pierwszego związku (53) widzimy, że 5, da się wyrazić ja- 
ko funkcya liczbowa wielkości a, oraz potęg В, , a więc w przypad- 


ku ruchu planetoidy typu Thule, należącej do klasy pierwszej, 
równania (53) i (56) najzupełniej określą nam wartość stałej a,. 


1) Brendel, |. c., rozdz. V sparsim. 
*) Brendel, Astr. Nachr. 3346; id. Th. d. kl. Plan. 1 rozdz. VII spars. 


http://rcin.org.pl 


О RUCHU PLANETOID. 83 


Wielkość W (str. 37) możemy napisać w postaci następu- 
Jace]: 
(57) W = (e,-- 1) v + wyrazy okresowe, 
czyli 

W =(C, +1 + fo) 2 + wyrazy okresowe, 
gdzie jak zwykle 
Tos TH. 

ale т, jest przynajmniej stopnia drugiego, więc obecnie możemy 
je opuścić, zatem 
(58) 1= +. 

Mielišmy zwiazek nastepujacy: 


| t = 0 — р К 
і 
w, = (1—8) 0— В – рМ, 
а міес 
w=w; +y К — (1 —)v — B"uW-uK, 
ale 
W = Ма ть 
w= (1— во — вто + wyrazy okresowe = 
= (1 — №.) о + wyrazy okresowe, 
gdzie i 
up) 
zatem czesé wiekowa W jest nastepujaca: 
(59) pars sec w=(1 + p) 0. 


Przypomnijmy sobie określenie wielkości p. i 1. 
Mieliśmy związek następujący: 
т п 
(60) $t оа? q = 
We wzorach (60) jak wiemy z poprzedniego, przez n ozna- 
czamy stałą całkowania; wielkość, którą Brendel nazwał „stałą 
ruchu“; n, zaś oznacza rzeczywisty ruch średni dzienny. 
Wprowadźmy dalej wielkości à i à, zdefiniowane w sposób 
następujący: 
_ 8—8. 3— à. 
= 4 я E 4 , 
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Z wymienionych przez nas małych wielkości à i à,, występu - 
je w naszych wyrażeniach, jako dzielnik tylko 2,. 

Brendel wykazał, że à, nie może być zerem w przypadku 
ścisłej wymierności, jest ono rzędu m’ o wykładniku = przeto za- 
wsze ono jest znacznie większe od masy zakłócającej. Co zaś do 
wielkości 2, to ona jest rzędu m' z wykładnikiem < 1 a w przy- 
padku wymierności staje się nawet równą zeru. Otóż zasadnicza 
myśl metody podanej przez Brendel'a polega właśnie na tem, 
że w wyrażeniach perturbacyi występują, jako małe dzielniki nie 
à lecz 2,. Brendel wykazał, że stosunek ścisłej wymierności mo- 
że zachodzić tylko pomiędzy wielkościami m i n' (a nie pomiędzy 
n, i n) wtedy 8—0. 

Na zasadzie wzorów poprzedzających mieć będziemy !): 


=p» (14+1)=p(1++00) 
3—8, —(3 — 2) (1 +1+00) 
3, =ô Co + 914-6 — 3с,— 33; 


ale 


poniewaz тату do czynienia z planetoida krytyczna, dla której 
0,550! Oraz. =, 
0, —21-4-2— 31, 
4p —3— 0, 


to mieé bedziemy: 


ale poniewaz 


to ostatecznie otrzymamy : 
(61) 0, =85—47p.. 
Zupełnie analogicznie możemy napisać, że 


gdy uwzględnimy (59), to п, możemy uważać jak ruch średni pla- 
netoidy w długości. Postępując zupełnie analogicznie jak poprze- 
dnio robiliśmy, otrzymamy: 


(62) ò, =ò — 4p 1—9,—4qy1,, 
Gdy ô, staje się bardzo małem lub nawet zerem, wtedy za- 


5 Brendel, A. N 3346. 
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chodzi przypadek libracyi, który w teoryi О yldén'a jest łatwy do 
uwzględnienia. n, właśnie określa nam ten ruch średni dzienny, 
który powoduje libracyę; musimy też znać dokładnie tę wielkość m, , 
aby módz wyznaczyć czas obiegu planetoidy. Należy zauwa- 
żyć, że р, i wą różnią się tylko o stałą yə. Wielkość y, jest dru- 
giego stopnia, a więc obecnie ją opuszczamy, podczas gdy += 
jest stale dodatnia i w przypadku planet charakterystycznych 
pierwszej klasy, do jakich zreszta nalezy Thule, jest juz wielko- 
$cig stopnia zerowego, a wiec w obecnych warunkach musi byé 
ona uwzgledniona. 


У. Zastosowanie liczbowe otrzymanych wzorów 
do ruchu planetoidy (279) Thule. 


Planetoide Thule (279) jako gwiazde 13.5 odkryt w Wiedniu 
J. Palisa dnia 25 pazdziernika 1888 r. Po kilkunastu dniach 
obserwacyi, ten astronom odrazu sie zoryentowal po wstepnym ra- 
chunku, Ze droga tej planetoidy lezy w bardzo znacznej odleglo- 
ści od słońca, mianowicie prawie 4 razy dalej niż średnia odległość 
ziemi od słońca. Przypuszczenia Palis'y były potwierdzone 
przez rachunki Langego z Berlina i Halma w Kiel. Nieby- 
wałe dotąd wśród małych płanet rozmiary drogi nowo odkrytej 
planetoidy wskazywały, że ta planetoida, znajdująca się na ze- 
wnętrznej części pierścienia planetoidalnego będzie bardzo cieka- 
wym objektem dla badań teoretycznych. Można było bowiem 
przewidywać, że planetoida Thule z powodu swego znacznego od- 
dalenia od Słońca oraz z powodu znacznego zbliżenia się do Jowi- 
sza — może podlegać okresowo nadzwyczaj silnym zakłóceniom 
spowodowanym przez przyciąganie Jowisza. Z tych wielkich per- 
turbacyi ruchu Thule, znając dokładnie jej pierwiastki drogi, 
z czasem można będzie obliczyć masę układu Jowisza, analogicznie 
jak to zrobił Newcomb, badając ruchy (33) Polyhymnii '), ale 
uczynić to będzie można, jak to słusznie zauważył Bidschof?), 


!' S. Newcomb. Astronomical Papers prepared for the use of the 
American Ephemeris vol. V, part V. 

2) Sitzungsberichte der Ak. der Wissenschaften, Wien Bd. 100: 
pg. 937. 
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dopiero wtedy, gdy mieé bedziemy przynajmniej 25 lat obserwa- 
cyi tej planety. Pomimo stabego blasku (1.35—15.5) jakim $wie- 
ci Thule, wskutek czego może być ona obserwowana jedynie przy 
pomocy wielkich narzędzi, obserwowano tę planetoidę podczas 
wszystkich opozycyi, które miały miejsce od chwili odkrycia, t. j. 
25/X 1888 г. aż do 10/11 1891 r. Te wszystkie obserwacye zostały 
bardzo szczegółowo opracowane przez Bidschofa!), który na 
ich zasadzie podał wartości pierwiastków drogi tej planetoidy, 
oraz efemerydy dla czasów bliskich opozycyi od czasu odkrycia do 
r. 1894. W pracy swej wspomina Bidschof, że 3 razy w ciągu 
stulecia następuje znaczne zbliżenie pomiędzy Thule i Jowiszem, 
w ruchu planetoidy występować muszą bardzo znaczne pertur- 
bacye. | 

Na zasadzie wszystkich obserwacyi, które mu byty dostepne, 
wyprowadzit Bidschof nastepujace ostateczne elementy: 


Epoka: 1881 Luty 20.0 cz. $r. Berl. 
L 104 21'30'".8 


M — 155 36 48.8 

т — 308 44 42.0 

8 = 75 26 12.1 

© — 933 18 299 równ. 1890.0 
== 2 29 34.2 | 

ф = 4 43 142 


log.a — 0.629 6674 
№ = 403”. 1860, 


oraz wyznaczył wielkość Thule dla średniej opozycyi, 13776. We- 
dług obliczeń Bidschof'a, najkorzystniejsze opozycye bywają 
w lipcu i sierpniu; wtedy planetoida osięga 13" 3, ale zato, dla 
obserwatoryów leżących w większych szerokościach północnych, 
jest ona wtedy dość niekorzystnie położona na niebie, gdyż jej 
deklinacya wówczas waha się pomiędzy — 24°. 5 a — 15°. 5. Ele- 
menty, które wyżej przytoczyłem zostały wyprowadzone na pod: 
stawie 8 miejsc normalnych. Te miejsca normalne są dość do- 


1) Sitzungsberichte d. Ak. d. Wiss. Wien Bd. 100, pg. 937. 
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brze przedstawione przez ostateczne elementy; najwieksze odchy 
lenie wynosi dla 
cos 6d a — +31. 1; 


do -—-— 3". 0. 


zaś dla 


W т. 1908 A. Wedemeyer!) ogłosił nowe rachunki, odno- 
szace sie do planetoidy Thule. Powodem podjecia na nowo tych 
rachunków bylo z jednej strony to, iz od chwili odkrycia do roku 
1906 Thule dwa razy obiegia swoja droge, z drugiej za$ strony, 
dzieki pracy Bidschof'a, który wprawdzie uwzglednit w swych 
rachunkach obserwacye do 1891 r. lecz podat efemerydy az do r. 
1894 i dalej, zdolano podczas wielu opozycyi zaobserwowaé (279). 
Ze wzgledu na ten nagromadzony materyal, A. Wedemeyer 
zdecydowal sie, wychodzac z elementów Bidschofa, opracowaé 
wszystkie obserwacye, jakie tylko byly zrobione od r. 1888 do r. 
1906 (ogólem 53), uwzglednié dokladne pozycye gwiazd porów- 
nawczych ina podstawie tego calego materyalu ugrupowanego 
w 14 miejsc normalnych wyprowadzié nowe elementy. 

W swoich rachunkach Wedemeyer postugiwal sie bardzo 
dogodna metoda Oppolzer'a, obliczania perturbacyi pierwiast- 
ków drogi.  Ostateczne elementy, otrzymane po wyrównaniu 
miejsc normalnych, byty nastepujace ?): 


Epoka i oskulacya 1907 grudzieñ 6.5 cz. $. B 
M — 210? 15' 47". 26 
w — 234 28 20.52 | 
Q = 75 35 57.14 1910.0 
= 2 2229.98 | 
= 4 37 35.94 
w = 404" 29204. 
Obserwacyi pózniejszych nie uwzglednionych przez Wede- 
meyera wykonano niewiele, mianowicie w 1906 r. dn. 17/X 


w Wiedniu i w 1911 r. w dniach 20 i 30 kwietnia fotograficzna 
w Heidelbergu?). Planetoida coraz bardziej przyblizala sie do Jo- 


1) Archiv. der deutschen Seewarte XXXI Jahrgang. 
23) Elementy Wedemeyer'a służą za podstawę niniejszych rachunków. 
3) A. N. 4591, 4563. 
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wisza przez co perturbacye jej biegu stawaly sie coraz znaczniej- 
sze i zgodność pomiędzy elementami Wedemeyera a rzeczy- 
wistością stawała się gorszą. Krytyczne zbliżenie Thule i Jowisza 
następowało w r. 1912. W lipcu 1913 r. Viljev!) ogłosił nowe 
oskulacyjne elementy Thule. Te jednak elementy nie mogą być 
uważane za bardzo dokładne, gdyż obliczone zostały wprost z ele- 
mentów Wedemeyer'a przez dodanie tylko perturbacyi od osku- 
lacyi 1907 r., które obliczył przy pomocy jednej z metod perturba- 
cyi szczególnych bez uwzględnienia obserwacyi z lat 1906 i 1911. 
Viljev podał efemerydy na podstawie swoich elementów. Dla 
chwili opozycyi miejsce wskazane przez Viljeva znacznie się 
różniło od tego, jakie było podane przez B. J. na podstawie ele- 
mentów Wedemeyera, planetoidy jednak nie zaobserwowano. 
Elementy podane przez Viljev'a są następujące: 


Epoka i oskulacya 1913 czerwiec 17.5 cz. śr. Berl. 
M = 358" 35' 20". 7 


w = 220 43 389 | 
RD #6 lomo 


fm 2/31 126 | 
Ф = 3 39 490 
p = 397". 60002 


log а = 0.6337068. 


Aby więc zbadać ruch tej ciekawej planetoidy, podlegającej 
tak znacznem zaburzeniom, wydaje mi się, że jedyna droga jaka 
nam pozostaje, jest podać, przynajmniej w przybliżeniu, wyrażenia 
analityczne perturbacyi ogólnych aby umożliwić na długi czas 
obliczanie efemeryd, pozwalających obserwatorom odnalezienie 
tego ciała niebieskiego. Metoda Gylden'a-Brendel'a wydaje 
mi się tutaj najodpowiedniejszą ze względu na to, że odrazu, przy 
pierwszem obliczeniu perturbacyi, pozwala uwzględnić wpływ 
wyrazów zależnych od wyższych potęg masy Jowisza. Gdy nowy 
materyał obserwacyjny zostanie zdobyty, to wtedy okaże się ko- 
nieczność obliczenia nowych pierwiastków drogi Thule, a rozpo- 
rządzając już całym materyałem obserwacyjnym, obejmującym 


1) A. N. 4661. 
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znaczny przeciąg czasu, będzie można dopiero wówczas podjąć 
zagadnienie wyznaczenia; na zasadzie perturbacyi Thule, masy 
układu Jowisza. 

Do naszych rachunków użyliśmy elementów Jowisza według 
„Connaissance des Temps pour 1914“. Masę Jowisza m’ wzięłem 
według Newcomb'a, wynosi ona: 


1 


m = 1047.355 


(lg m' — 6.9799061). 


Poniewaz w niniejszej pracy ograniczymy sig do wyszukania 
przerwy (luki) dla planetoid typu 2: pozostawiajac do nastepnej 


części wyrażenia liczbowe spółczynników dających rozwinięcie 
perturbacyi Thule, dla tego też przytoczymy tutaj tylko rachun- 
ki, które odnoszą się bezpośrednio do tego zagadnienia. W roz- 
dziale drugim podaliśmy, zgodnie z teoryą, którą do naszego 
przypadku stosujemy, wyrażenia analityczne rozwinięcia funkcyi 
perturbacyjnej, w których występowały pewne arytmetyczne wiel- 
kości, dające się wyrazić przez stosunek a połowy osi drogi Thule 
do połowy osi drogi Jowisza. Te spółczynniki liczbowe były to 
wielkości ß i y zaopatrzone odpowiedniemi wskaźnikami. 

Przystąpmy do ich obliczenia. Z elementów Jowisza i Thu- 
le znajdziemy wielkość a: 


lg a=9.912 6572. 


Ta wielkość malo się różni od jedności i widzimy, że spól- 
czynniki @ i y, zależne od tej wielkości, nie będą się znajdować 
ani w tablicach Masal'a'), w których mamy tabułowane wartości 
B i y dla argumentu lg a od 9.30 az do 9.60, ani w tablicach G yl- 
dé n'a?), w których jest zawarty lg « od 9.620 do 9.800. 

Musimy więc bezpośrednio z elementów obliczyć wszystkie В. 
Na str. 24 mieliśmy wzór (8): 


z 
1) Masal H: Tables de l'intégrale paa . Astronomi- 
o fl—a? sin? p} 
ska iakttagelser och undersókningar ánstalda pá Stockholms observatorium, Fje- 
rde Bandet, Stockholm 1891, pg. 51. 


2) Gyldén Н.: Astr. Gesell. Publ. 21. Hülfstafeln zur Berechnung der 
Hauptungleichheiten der Kl. Planeten. 
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r 


1 i 
I g = 3! sin?" ody 


к. [I—a? sin? çË 
Znając a możemy przystąpić do obliczenia wszystkich f . 
Gylden!), Harzer?) obliczali te wielkości w ten sposób, 
że przedewszystkiem znajdowali wszystkie @ dla n=1. Przy po- 
mocy pewnych wielkości pomocniczych fa, À, i #,, określonych 
w sposób następujący: 


fw ee MENS ne, Le 
" Q@n+2)@2n+4 (Ha)? " 2n+2 FEL 
I B 
Ti "n Я 
Wielkość >, można rozwinąć w шашек ciągły : 
1 
РА 
1 — fayi 
I—fnte 
1—... futm—1 


1—f n Ya m1 | 
Użycie tego wzoru jest następujące: najpierw oblicza sie >, 
z największem m zachodzącem w naszych wzorach, wprost ze 
wzoru: 
=: b 
pr jx 7 
l— fer: a 
1— fcm 19 m 1 
a potem przez koleine stosowanie wzoru redukcyjnego 
1 
1 — fn LASS! 
otrzymuje sie wszystkie poprzedzajace y. 
Wprowadźmy jeszcze wielkość pomocniczą p,—,9,, wte- 
dy mieć będziemy wzór następujący: 


Bl, = В Do D1 Р»... Du, 


Un = 


1) Gyldén Н. Orbites absolues I, livre Ш. 
*) Harzer P. Untersuchungen etc. p. 26. 
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który może posłużyć do tego, aby po znalezieniu f*, obliczyć 
wszystkie f!,. Potem wszystkie inne 8, otrzymamy na zasadzie 
następującego wzoru rekurencyjnego: 
п) KU = + 

Trudność w obliczeniu wielkości В’, polega na tem, ze uła- 
mek ciągły, wyrażający 9,, jest w niektórych przypadkach mało 
zbieżny. Harzer zauważył, że dla planetoid zwykłych wystar- 
czy przyjąć ostatnie п = 12; ale w naszym przypadku stosunek a 
jest mało różny od 1, wielkości f, maleją bardzo powoli, nawet 
dla n — 30 jeszcze nie możemy osiągnąć tego, aby dwa następują- 
ce po sobie przybliżenia w rozwinięciu ułamka?), różniły się od 
siebie o wielkości zaniedbywalne, musimy więc, aby módz obli- 
czyć ß!„, postąpić inaczej, bez obliczenia wielkości pomocniczych 
A, i fa. W tym celu przedewszystkiem obliczymy spółczynnik f!;; 
możemy to zrobić najlepiej przy pomocy metody średniej arytme- 
tyczno - geometrycznej. Na tę metodę poraz pierwszy wskazał 
C. F. Gauss*), jako na praktyczny sposób obliczania całki elip- 
tycznej zupełnej, mamy bowiem: 


кет f nS 
° 84 Yı-a?sin?g ' 
0 
jako pierwsze przyblizenie kladziemy 
E S E y—V1—ao? 
obliczmy z tych dwu wartości średnią arytmetyczną i geome- 
tryczną: 
а = zał y=Vry, 
następnie z z' i y! znowu średnią arytmetyczną i geometryczną: 
r , mes 
E TU. у" = Vary (Жа; 
az przy pewnej wartosci m znajdziemy: 
gm — у“ 


1) Gyldén H. Uudersókningar af theorien för himlakropparnas rörel- 
ser. Il, str. 23, (Bihang till sv. Vet.-Ak. Handl. Bd. 6). 

+, Hansen zauważył, że dla п dużego Ya = 1 + a2. 

$) Gauss C. Е. Werke Ш, pg. 352. 
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wtenczas: 
1 1 


Po ze) == ym ^ 
Ta metoda obliczenia В', jest bardzo praktyczna i nadzwy- 
czaj predko prowadzi do celu. 
W naszym przypadku mamy: 


lg z—0; 1#И1—а? = 9.7600268 
2) —9.8882366 y!" — 9.8882666 
lg B! = 0.1117634. 


Najwyższa wartość wskaźnika n jaka nam jest potrzebna jest 
== 12, obliczyć więc musimy &,! dla n—1 ..., 12. 

Te spółczynniki otrzymamy najprędzej przy pomocy metody 
podanej przez Ludendoríf'a!). Tablice pomocniczą podaną 
przez Ludendorffa rozszerzyłem tak co do rozciągłości, jak 
i stopnia dokładności i przy pomocy jego wzorów znalazłem sto- 
sunkowo łatwo pozostałe fl, (n—1,2...12). Dla kontroli drugi 
raz te same spółczynniki f!, (n=0, 1, 2... 12) zostały obliczone 
przy pomocy kwadratur mechanicznych. Tę metodę do obliczenia 
spółczynników В" poraz pierwszy zastosował Harzer, później zaś 
O. Callandreau. Dokładność w rachunkach jest znaczna i kon- 
trola łatwa, gdyż każdy spółczynnik jest obliczony jako średnia 
arytmetyczna z dwóch sum mało różniących się od siebie; ze zgo- 
dności tych dwóch rezultatów widzimy, czy nie został popełniony 
jakiś gruby błąd rachunkowy. W naszym przypadku kąt 90° dzie- 
lilem na przedziały po 5° i dopiero wykonałem całkowanie. Po- 
nieważ cały rachunek wykonałem według wzorów podanych przez 
Harzer'a i O. Callandreau, dlatego też w bliższe szczegóły, 
odnoszące się do tej metody, nie wdaję się, odsyłam do rozpraw 
skąd one zostały wzięte 2). 

Wartości 81, obliczone za pomocą kwadratur mechanicznych 
wypadły zupełnie zgodne z wartościami obliczonemi przy pomocy 
powyższej tablicy. Aby otrzymać pozostałe spółczynniki $',,, 6%, 
Bias Piro, Po, które mają posłużyć jako punkt wyjścia do użycia 
wzoru II, znów zastosowałem metodę kwadratur mechanicznych, 


a więc 


1) Ludendorff. А. N. Bd. 160. 
*) Callandreau. Journal de l'Ecole Polytéchnique, Il série, t. VII; 
Harzer, 1. с.; pg. 26. 


http://rcin.org.pl 


О RUCHU PLANETOID. 93 


oraz kontrolę przy pomocy tablicy Ludendorffa; zacząłem od 
obliczenia powyższych spolczynników dlatego, ponieważ wedle 
uwagi Harzer'a') najkorzystniej obliczyć spółczynniki 8 naj- 
pierw dla najwyższej wartości n, gdyż w ten sposób najlepiej 
zmniejszymy wpływ błędów, pochodzących od zaokrągleń. Gdy 
już spółczynniki 8 zostały znalezione, zastosowałem wzór Il 
i w ten sposób otrzymałem wszystkie pozostałe spółczynniki. Aby 
skontrolować cały rachunek obliczyłem kilka dowolnie wziętych 
spółczynników metodą kwadratur. Rezultaty otrzymane różniły się 
tylko o pare jednostek w 7-ej cyfrze logarytmów, co dla naszego 
przypadku, gdzie się ograniczamy do rachunku 6-cio cyfrowego, 
jest najzupełniej obojętne. 

W ten sposób znalezione logarytmy spółczynników £ są na- 
stępujące: 

Logarytmy spółczynników 6, 

n/s 1 3 5 7 9 
0.111763 0.388660 0.719219 1.092421 1.495190 
9.887778 0.236762 0.620523 1.027947 1.451216 
9.761532 0.153599 0.563644 0.986993 1.420154 
9.693387 0.102476 0.524368 0.955978 1.394980 
9.642434 0.062590 0.492468 0.929894 1.373181 
9.601597 0.029631 0.465363 0.907191 1.353821 
9.567459 0.001428 0.441683 0.887002 1.336327 
9.538099 9.976724 0.420599 0.868721 1.320326 
9.512324 9.954718 0.401561 0.852056 1.305577 
9.489344 9.934848 0.384186 0.836699 1.291872 
10 9.468604 9.916727 0.368193 0.822446 1.279062 
11 9.449703 9.900062 0.353368 0.809138 1.267029 
12 9.432337 9.884629 0.339546 0.796650 1.255680 


© o NO сл > G N — O 


Gdy już znaleźliśmy wszystkie wielkości I, to na zasadzie 
wzoru następującego 


AE DE Ob (26 — 1) „агт gasta 
A 99 Te 


otrzymamy spółczynniki 7... 


Ш Yes 


1) ‚Harzer P, I. c. pg. 26. 
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Te spólczynniki sa nastepujace: 
Logarytmy spólczynników т». - 


n/s 0 1 2 3 4 
0 0.024421 9.673703 9.700962 9.839429 0.023954 
1 9.713092 9583198 9.574342 9.726002 9.917251 
2 9.499504 9.364732 9.455100 9.615956 9.812414 
3 9.344016 9.237503 9.340652 9.508424 9.709071 
4 9.205720 9.117201 9.229629 9.402800 9.606979 
5 9.077540 9.001658 9.121202 9.298793 9.505931 
6 8.956059 8.889613 9.014822 9.196093 9.405778 
7 8.839356 8.780260 8.910104 9.094497 9.306403 
8 8.726239 8.673047 8.806768 8.993846 9.207710 
9 8.615916 8.567583 8.704600 8.894015 

10 8.507834 8.463575 8.603436 8.794901 

11 8.401590 8.360800 8.503145 8.696416 

12 8.296881 8.259083 8.403622 8.598585 


lg тоо = 9.859395. 


Znaleźliśmy w ten sposób spółczynniki &,9 i gna, możemy 
obecnie przystąpić do obliczenia wielkości A i B opatrzonych od- 
powiedniemi wskaźnikami u góry i u dołu. Te wielkości podsta- 
wowe 4 i B mogą być obliczone dwoma sposobami. Jeden stano- 
wi zastosowanie wzorów (20) i (24) (str. 26, 27) przy pomocy 
wielkości pomocniczych Q», m, p które sa funkcyami 4, a, oraz wiel- 
kości P, m, pi Qu mp, Znów będących funkcyami 9, „ ,. Podo- 
bne postępowanie Brendel zaleca; jest ono istotnie bardzo prak- 
tyczne, gdy chcemy odrazu obliczyć wszystkie potrzebne spół- 
czynniki w rozwinięciu funkcyi perturbacyjnej dla danćj planetoi- 
dy. W naszym jednak przypadku podobny sposób postępowania: 
nie zalecałby się, gdyż ograniczając się do wyrażeń zerowego 
stopnia, daleko prościej odrazu obliczyć spółczynniki A i B wprost 
we funkcyi us. Wzory stosowne dla takiego rachunku zostały po- 
dane przez Masal'a!). Są one wygodne w użyciu o ile nie posu- 


1) H. Masal: Formeln und Tafeln zur Berechnung der Absoluten Stör- 
gen der Planeten, Kongl. Svenska Vetensk. Ak. Handlingar. Bd. 23. 
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wamy się w rozwinieciach do wyższych stopni, gdyż wtedy stają 
się one skomplikowane. W granicach dokładności przyjętych przez 
nas, wzory Masal'a są następujące: 


Алоо = — 2N Yno 
Boo = + Mino + 415a 
AG = Tao [+ 2n) — 2n? q] + 227. 
BR = Yno [— (m+n) J-2n* €] + inil (4n + 6)+-4n. q] — 
— Bus 
AGA = To [ (n + 2n) + 2n? ç] + 2n ть 
BC = Tao [— (n? 4-n) — nę] ra [— (4n +6) — Ang] — 
—8 Tas 
Año = Tno(2n74-4m) 4-42 ns 
Bi, = — Tno 20 (n+ 1 — qua (80 + 12) — 161, 
Arto = Tuo (nt + 5n? -- 6n) — тл (4n* + 14n) — 8n» 
BA, = Hron (n+ 1) (n +2) 2-104 6 (n 4- 2)? + 
+ 142. 24(n--3)--481,5. 
gdzie ç = =” . Na zasadzie tych wzorów otrzymali$my nastepu- 


jące wartości logarytmów spółczynników 4 i B. 


n 0 4 8 12 
Anoo 0.108810, 9.930359, 9.677 092, 
Воо 0.275763 0.257352 0.017139 9.738 837 
4, 0.991403 1.001115 0.876663 
pi^. 1.136647, — 1.173232, 1.088907, 0.931817, 
ba^. 1.682384, 1.834337, 1.818637, 
ps 1.906612 1.946349 1.945695 1.972880. 
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Znaleźliśmy więc ostatecznie potrzebne spółczynniki rozwi- 
nięcia funkcyi perturbacyjnej; możemy obecnie przystąpić do wy- 
rachowania przerwy dla planetoidy typu Thule w pierścieniu pla- 
netoidalnym. 

Mieliśmy równanie (37): 


(1) (28, 4- 28,7) в, — p' Fp 6, + p" R’, 


które zawiera tylko jedna niewiadoma B, (8, wyraza sie przez ß,) 
przy pomocy zwiazku (61) str. 84 i jest wzgledem tej niewiadomej 
stopnia trzeciego; (jak zobaczymy spólczynnik przy B, w tym rów- 
naniu zawiara 2). 

Brendel okazat, ze w celu znalezienia obu granic przerwy 
nalezy wyznaczyé ten spólczynnik, przez nadanie odpowiedniej 
wartości na 2, w taki sposób, ażeby to równanie dawało na 6, tylko 
dwie, od siebie różne wartości, wówczas dwie wartości 2,, odpo- 
wiadające tym dwum różnym wartościom f,, określają właśnie gra- 
nice przerwy. 

Brendel!) wykazał, że В, w przypadku planetoidy charak- 
terystycznej, nigdy nie może się stać większem co do rzędu wiel- 
kości od pierwiastka trzeciego stopnia z masy zakłócającej; gdy 
mamy do czynienia z planetoidą tego rodzaju, że stosunek ruchu 
średniego dziennego planety zakłócającej do takiegoż ruchu dzien- 
nego planetoidy wyraża się ściśle przez prosty ułamek, wtedy 


$—0, zaś 8, =--4вт1; 
ale jak to już poprzednio widzieliśmy 


3 
FRZ ви". 
Poniewaz В, w przypadku planetoidy charakterystycznej со 


RR 
najwyżej może być rzędu Vm’, przeto najmniejsza wartość (à,) 


jest rzędu mniejszego (co do rzędu wielkości) niż m'*. Jeżeli mia- 
nowicie ta wartość graniczna byłaby przekroczoną, to wtedy nasze 
rozwinięcia stalyby się rozbieżne 2), ruch planetoidy, któraby mia- 
ła takie 0, byłby niestały. Dla tych krytycznych wartości à, 


1) Brendel. А. М. 3346. 
2 Brendel, l. c; Brendel. Theorie d. kl. Pl. I, Kap. VII; 
Gyldćn Acta Math: t. IX. 
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a więc i B, występuje (w znaczeniu Brendlowskim) przerwa (la- 
kuna) w pierścieniu planetoidalnym; przerwa ta przychodzi do 
skutku w ten sposób, że à, zrywa ciągłość; wartości jego po obu 


stronach przerwy są przeciwnego znaku. Ponieważ 2,, jak to wy- 
, 
nika z poprzedniego, nigdy nie może być zerem, przeto y, = 5 


nigdy nie może się równać jakiemuś prostemu ułamkowi. Nato- 
miast à może się stać równem zeru: ma to miejsce wtedy, gdy y 
jest równe prostemu ułamkowi. dla tego też we wszystkich na- 
szych dotychczasowych wzorach, gdzie zachodzi mały mianownik 
2,, nie wolno zastąpić go przez à, jakkolwiek à może być bliz- 
kie à,. 

Stosownie do tego, cośmy wyżej powiedzieli, zajmijmy się 
przedewszystkiem wyznaczeniem wartości В,. W tym celu weźmy 
na uwagę przytoczone przed chwilą równanie (1). Ponieważ 2, 
nie jest rzędu zero, to 2,78, jest rzędu wyższego, niż pozostałe wy- 
razy w tym wzorze, możemy je więc opuścić. Zamiast 8, pod- 
stawmy 

ò — 4p y= ô— buh’, 


wtedy otrzymamy: 
268, — 12u p =p +p" 8, 3-2" B, 
co możemy napisać: 


; p р"— 28 q n 2 = 
czyli 
(3) Bi” +28, +78, =0, 
gdzie dla krótkości położyliśmy: 
TS p" А ^T p'—25 "i m. 
(5 Cams om tt» 


Z (4) widzimy, ze tylko spölczynniki ö oraz k sa niezalazne 
od à. Równanie (2) musimy tak przekształcić, aby sig pozbyć wy- 
razu 8,2; w tym celu połóżmy: 

(a) = ¬ 3 
wtedy znajdziemy: 


О ruchu planetoid. 


m 
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(5) b `+ db, +e=0, 
gdzie oznaczylišmy: 
o. $ 
(6) d=j— z 
oraz Mh 
miss. 


Równanie (4) bedzie posiadalo pierwiastek podwójny wtedy, 
gdy wyróznik równania (5) jest równy 0. Po obu stronach przer- 
wy, b, ma wartości przeciwnego znaku; ponieważ В, —b, jest małe 
w porównaniu z b,, przeto także i wartości В, będą przeciwnego 
znaku. 

Wyróżnik ten jest 


4 e 
Карат 9: 


z teoryi równań trzeciego stopnia wiemy, że gdy D > 0, wtedy 
mieć będziemy tylko jeden rzeczywisty pierwiastek rówńania (5), 
a więc dla B, otrzymamy tylko jedną wartość; gdy D —0, wtedy 
otrzymujemy trzy rozwiązania, z których dwa są równe i wartość f, 
znaleziona z tego warunku daje nam właśnie krytyczną wartość 2; 
wreśzcie gdy D < 0 mieć będziemy wszystkie trzy pierwiastki rze- 
czywiste, ale tylko jeden, jak to możemy się przekonać przez pod- 
stawienie, spełniać będzie nasze równania. 

Aby znaleźć ten krytyczny warunek, z otrzymanych poprze- 
dnio rozwinięć spółczynników spółrzędnych Gyldćnowskich 
obliczmy stałe wielkości À oraz k; w naszym przypadku znale- 
zlismy: 

lgi=8.442234_0; Igk=6.672816_1 


te wielkości są stałe; wielkość d zależy od j, a więc od 2, po sze- : 
regu prób dokonanych co do wartości j znajdujemy, że warunko- 
wi D=0 czyni zadość ta wartość j, która daje ze wzoru (5) 


1gd—8.19508. 9. i Ige—7.773160 10, 


skąd znajdujemy: 
lg 5,'— lg 5," = 8.824044 
zaš 
Ig б,''——= 9.125074, 
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na zasadzie wzoru (a) z łatwością otrzymamy: 
Ig(Bh = 8.759360; Ig(8,)u= 9.154128, 
te wartości B, за krytyczne. Znając wartość krytyczną j z latwo- 
ścią znajdziemy tej wartości odpowiadającą wartość 2. W naszym 
przypadku znaleźliśmy: 
lg 2— 8.559884 ; 


(Gr i (hu 
otrzymamy dwie wartosci dla y na zasadzie nastepujacego wzoru, 
dokładniejszego niż ten, który poprzednio stosowaliśmy !): 


"IKE ZY 
0 кы NÉS 


znajac 


Wartosci te w przypadku naszym sa nastepujace: 
lg yı = 7.696539; lg yu = 8.495577. 


Teraz ze wzoru (2), rozdzialu poprzedzajacego, znajdziemy 
dwie warto$ci ö,, odpowiadajace tem wartosciom y, mianowicie : 


1g(2,))— 8,334772;  1g(2,)ii = 8.750883, 
dalej ze wzoru 
LEM o 
MATE 
znajdziemy obie warto$ci prawdziwego ruchu Sredniego, które 
ograniczają szukaną przerwę dla typu */,. Te wartości są: 


(n) = 401". 732; (n4) = 391". 485 


wreszcie stałą n, dla wartości krytycznych 2,, znajdziemy ze wzoru 
następującego: 
_ 4n 
„gg 
skad: 
n — 403". 723. 

Rezultaty otrzymane możemy więc wysłowić w sposób na- 
stępujący: 

Dla planetoid typu 3/, wartość B, musi spełniać 
warunek następujący: 


1, Brendel А. N. 3346 wzór № 38. 
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lg 8, <8.759360 lub lg 8, =9 154128, 
lg 8, <8.334772 lub lg 3, >8 750883, 


zatem prawdziwy ruch dzienny n, dla typu Thule 
musi byé 


zaś 


n, > 401". 732 lub m, < 391". 485. 


Pomiedzy temi warto$ciami m, znajduje sie przerwa w zna- 
czeniu Brendelowskim w pierscieniu planetoidalnym; planetoida, 
która posiadalaby ruch Sredni dzienny prawdziwy (w definicyi 
Brendelowskiej), zawarty pomiędzy powyższemi granicznemi war- 
tosciami n, poruszalaby sie po drodze niestalej. 

[Musimy tu zwrócić uwagę, że ze względu na definicyę m, ta 
wielkość nie jest bezpośrednio porównywalną z wartością ruchu 
średniego dziennego oskulacyjnego, jaki znamy z elementów pla- 
netoidy]. 

W ten sposób otrzymaliśmy wartości krytyczne prawdziwe- 
go ruchu średniego dziennego, warunkującego przerwę dla typu */,; 
taka przerwa dotąd była wyznaczoną dla Hekuby i Hildy (zrobił 
to Brendel!) Prace Schwarzschild'a dla Hekuby*) i Dzie- 
wulskiego*) dla Hildy, wykonane na zupełnie odmiennej dro- 
dze, mianowicie metodą Poincarć'go wynajdywania rozwiązań 
okresowych, potwierdzają rezultaty otrzymane przez Brendel'a 
jego własną metodą, co daje nam najlepsze świadectwo o war- 
tości tej metody, stosunkowo łatwej w zastosowaniach prak- 
tycznych. 


' Brendel. A № 3396. Buchholz H. Untersuchungen etc. 
1) Schwarzschild. A. N. Bd. 160. 
3) Dziewulski. A. N. Bd. 181. 
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RÉSUMÉ. 


JAN KRASSOWSKI. 


Sur le mouvement des petites planètes du type de Thule (2). 


Détermination des termes élémentaires et de la lacune 
pour le type Thule. 


Le nombre croissant des petites planètes a suscité l’appari- 
tion des méthodes spéciales en mécanique céleste qui permettent 
de calculer plus ou moins rapidement les perturbations. 

Gyldén et ses nombreux élèves ont divisé les petites pla- 
nètes en groupes, appelés ,types“ où la relation de commensura- 
bilité entre le mouvement moyen diurne de la petite planète et ce- 
lui de Jupiter joue un rôle très important. On sait que si ce rap- 
port peut s'exprimer approximativement par un rapport simple 
des deux nombres entiers, les perturbations, causées par la planète 
principale, deviennent de plus en plus grandes à mesure que le 
rapport de ces deux moyens mouvements diurnes s'approche de 
la commensurabilité exacte. Jupiter produit dans l'essaim de 
petites planétes en général des grandes perturbations gráce à sa 
masse énorme. Les perturbations Joviennes deviennent déjà trés 


sensibles pour les planétoides du grouppe 3 (moyen mouvement 


diurne ca 900”) et du grouppe Y (moyen mouvement diurne ca 
600"). 

Si la commensurabilité devient trés approchée, les perturba- 
tions croissent trés vite et les méthodes ordinaires du calcul des 
perturbations deviennent insuffisantes ou trés difficiles et pénibles 
dans les applications, surtout s'il s'agit des planétoides-types. 
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C'est à Gyldén que nous devons une des méthodes du 
calcul approché des pérturbations. Cette méthode, remaniée et 
perfectionnée par ses éminents éléves (M. M. Backlund, Bren- 
del, Harzer, Wolf etc.), peut étre méme appliquée dans les 
cas les plus difficiles où la commensurabilité des moyens mouve- 
ments de la planété principale et de la petite planéte est à peu 
prés exacte. 

Cette méthode a fait ses preuves; elle a été déjà plusieurs 
fois appliquée, avec bon résultat, aux différents cas difficiles du 
calcul des perturbations des petites planètes appartenant aux grou- 
pes d'Hécube, de Hestia ou Hilda. 

M. M. Brendel et Kramer nous ont donné des tables 
étendues, qui nous facilitent notablement le calcul des perturba- 
tions dans le cas des types Hestia et Hécube et permettent rapide- 
ment de déterminer leurs positions géocentriques. Les calculs de 
M. Harzer ont été sur une autre voie brillamment controlés par 


M. Simonin pour Hécube. Le cas difficile du type de Hilda 5 | 


a été commencé d'étre étudié par M. Buchholz, qui a déterminé 
les termes caractéristiques, sans donner l'expréssion méme des per- 
turbations. 

La petite planète Thule (279) appartient au type */, des péti- 
tes planètes et jusqu'à présent c'est l'unique réprésentant de ce 
type. Thule peut s'approcher de l'orbite de Jupiter beaucoup 
plus que Hilda. M. Bidschof trouve que des grandes proximi- 
tés entre Jupiter et Thule arrivent jusqu'à trois fois dant l'éspace 
d'un siécle. Pendant un pareil rapprochement le mouvement de 
Thule est fortement perturbé par la puissante action de Jupiter, 
qui produit des variations trés notables dans son orbite. Ainsi 
pendant la proximité qui a eu lieu en 1912, les calculs de M. Vil- 
jev nous montrent que le moyen mouvement de Thule a subit de 
trés fortes variations. 

M. Kobb, en se basant sur les travaux de M. Darwin sur 
les orbits périodiques, a cru pouvoir admettre que l'orbite de Thu- 
le peut devenir instable et s'entrecouper méme avec celle de Ju- 
piter. 

Jusqu'à présent le mouvement de cette petite planéte, appar- 
tenant à un type planétoidal si intéréssant, n'a pas été étudié 
d'une maniére générale. A ce que je sais, M. M. Bidschof et 
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Wedemeyer ont seulement calculé et corrigé les éléments de 
Thule. M. Viljev de son côté n'a fait que le calcul des pertur- 
bations de son mouvement, limité à un court espace de temps en 
se basant sur les éléments de M. Wedemeyer. 

Il m'a semblé intéressant d'essayer d'appliquer à Thule les 
méthodes de Gyldén, modifiées et perfectionnées avec tant de 
succés par M. Brendel. L'orbite de départ de nos calculs a été 
celle de M. Wedemeyer. 

Dans le présent mémoire je donne un aperqu succint de la 
méthode de Gyld én, modifiée par M. Brendel, je détermine pour 
le type Thule les termes caractéristiques ainsi que la lacune (dans 
le sens de M. Brendel) dans lanneau planétoidal, qui jusqu'à 
présent n'a pas été determinée pour ce type. 

Pour ne pas étre trop long je me permets de renvoyer le lec- 
teur aux formules contenues dans le texte polonais, aux quelles, 
du reste, je fais les renvois donnant les numéros des formules et 
les pages correspondantes du texte. 


1. Dans le premier chapitre je me suis attaché d'exposer, 
aussi brièvement qu'il m'a été possible, les idées générales sur les- 
quelles repose la théorie de Gyldén ainsi que les équations les 
plus importantes de la théorie qui nous intéresse, 

On sait que c'est la méthode de Hansen qui a sérvi à G yl- 
dén du point de départ. Le équations fondamentales sont: 1 pg. 
4; II pg. 8; et III pg. 9, texte pol. Dans ces équations v désigne 
la longitude vraie de la planéte et r son rayon vecteur. 

Les équations de Hansen en coordonnées polaires, G yl- 
dén transforme, en introduisant des nouvelles variables p, n S 
et v, qu'il appele ,longitude vraie“ et qui est une variable in- 
dépendante. Les équations nouvelles ainsi obtenues sont: IV, 
ГУ? pg. 11, V, V*is pg. 12. Nous introduisons maintenant les no- 
tations de M. Brendel!) et nous rappelons la définition des ter- 
mes élémentaires et caractéristiques à longue et courte 
période. Parmi ces termes nous distinguons, suivant M. Harzer?) 
ceux qui appartiennent au genre А. B, C, D et ceux que M. Bren- 


1) Brendel. Theorie d. kl. Planeten I. 
2) Harzer P. Untersuchungen etc: pg. 3. 


http://rcin.org.pl 


106 J.' KRASSOWSKI, 


del appelle ,rein erster Ordnung*. Avec M. Brendel nous 
désignerons par les lettres latines les termes du premier ordre, et 
par les lettres grecques ceux qui contiennent, comme diviseur, la 
petite quantité 2. 

Dans le présent mémoire nous nous bornons à ne considérer 
que le mouvement de la petite planéte qui s'effectue dans le plan 
méme de l'orbite: c'est pour celà que nous laissons de cóté les per- 
turbations de l'inclinaison 

En nous servant des rélations connues, données par M. Bren- 
del, entre les quantités Gyldéniennes 1, П, А nous obtenons: 


= am 28 
(1) "ie l+n cos У 


où (r) c'est le rayon vecteur „absolu*. 
Aprés quelques transformations on parvient à l'équation sui- 
vante: 


(1 — 3»? 
(2) Toute res cos nv. 

Cette équation obtenue par Gyldén est trés importante, 
car le second membre ne contient aucun élément à longue période, 
tandis que le premier membre de cette équation est formé par une 
fonction pareille. 

Pour achever toutes les transformations nécessaires on intro- 
duit encore le „temps reduit“ (t) qui est donné par l'équation 
suivante: 


(3) о A —v-|-XB,sin nv où 
L=nt+4-A 


et la quantité E est définie par: 


dE _ y . B, cos nl 


dE y dB, sin nll 
dv dv 


sin nv — Z — cos nv 
dv y 


alors (12) devient: 
at) | ü—3»9 ‚ав 


(4) dv (I+ncos v) ' dv ` 


Si W nous donne la différence entre l'angle déterminé par le 
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temps vrai nt et celui qui correspond au „temps réduit“ 
m(t), nous aurons encore la relation suivante: 


(5) W —nt—n(t) 


qui enfin nous donnera l'équation différentielle VI pg. 20 où ent- 
rent aussi les quantités R et S. 

Pour les applications de la méthode les relations IV, IVP, V, 
Vois et VI sont les plus importantes. Nous allons maintenant tâ- 
cher d'intégrer ces équations. Ces équations contiennent les quanti- 
tés P et Q qui sont des dérivées partielles par rapport à v et à r de 
la fonction perturbatrice Q que nous devons développer en série 
suivant les procédés, donnés par Gyldén. 


2. On sait que Gyldén pose: 


(6) + = R, + 2R, cos H+-2R, cos 2 H+-... 
où 

2.0! sin" ed | | 
(7) Ем — r! ść | à sf 22 sin? g + а?) sin? y 


en substituant dans cette expression les coefficients numériques 
: a 
connus Yns, ne dépendant que de а (on se rappele que e= le 


rapport des demi-axes de la planëte troublée et troublante) on 
trouve: 


, " 
(8) R= (1 — 3) nme. 


Ce développement converge encore dans le cas de Thule. 
En introduisant dans (6) l'expression trouvée, nous obtenons 


le développement en série seulement de la partie A de la fonction 


perturbatrice, qui, comme on le sait d'ailleurs, étant complète, 
s'exprime au moyen de la relation suivante: 
a — 


(9) a (Q) — m’ [ 


À 2 a? (1 Lx cos H] 


Si nous développons cette expression en série, en introdui- 
sant les coefficients numériques connus, nous obtenons: 
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(10) a (9) —2m' b © (n, s, 8), v pen" q” cos nH 


oü n, 85,8, Y, V 
(11) Qa =m > Q(n.s.s') v.v pe py» 1/27. 


Ces coefficients sont des fonctions numériques de a et peu- 
vent s'exprimer au moyen des 7„,. !) Cette relation est donnée pg. 
25 texte pol. (20). 

Les quantités P et Q peuvent s'exprimer, aprés quelques 
transformations, au moyen des relations suivantes: 


+ 
P=22P mas, 1 pp" 12° 2" cos n(o—v' 
(12) ( pp" т ( ) 


9=2 Оль) w pip” qq? sin n(v — v) 


(où * signifie que pour n=0 le facteur 2 doit être omis). 

Les coefficients de ces expressions sont données par les for- 
mules pg. 27 — 30 texte pol., qui n'ont pas besoin d'aucune expli- 
cation. 

Dans les expressions (12) nous écrivons p en séparant la par- 
tie élémentaire du genre B et les termes caractéristiques: 


(13) p=(p) +R 


où (p) désigne les termes élémentaires du genre Bet P contient 
les termes ordinaires et caractéristiques. En ne considérant dans 
р que les termes caractéristiques seuls et les substituant dans P 
et Q, nous obtenons l'expression (39) pg. 31 texte pol. pour P et 
une expression analogue pour Q. 

Pour avoir les expressions définitives des quantités P et Q 
il est nécessaire d'effectuer encore une transformation, parce que 
dans ces quantités entrent les S, R, W, qui de leur côté contiennent 
les termes des genres différents (ordinaires, à courte et longue pé- 
riodes, élémentaires etc.). En effectuant des transformations ana- 
logues aux précédentes (pg. 32 texte pol.) on obtient les expres- 
sions définitives pour P et Q (pg. 32, 33 texte pol.) 2). 


') Brendel, Th. d. kl. Plan. I, pg. 52; Harzer, Untersuch., pg. 49; 
Kramer I, Th. d. kl. Pl. pg. 9. 
3) Brendel, c. pg. 70; Harzer: I. с. pg. 30. 
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3. Nous allons déterminer maintenant les termes caractéri- 
stiques et élémentaires dans les cas de Thule (279). 

Aux pg. 12 et 14 texte pol. nous avons donné les équations 
différentielles qui déterminent les quantités 2 et p; l'équation dif- 
férentielle en W leur est analogue: 


(14) W=c,-+ y sin (,v— B,). 


Dans toutes ces équations а„ Du, c, sont des petites quan- 
tités, À, est un petit diviseur. Quelquefois cependant, après avoir 
effectué l'intégration, il arrive que malgré que les quantités а», On, 
c, sont petites, les termes renfermant le diviseur X, peuvent être 
ou peu différents de zéro ou trés grands!), car, aprés cette opéra- 
tion, les termes considérés peuvent avoir comme facteur x ou 
Is Si À, est peu différent de zéro, l'argument de la fonction 
cosinus est un trés petit angle, ce qui nous donne pour le terme 
considéré une période trés longue. 

Si À, diffère peu de 1 alors la période du terme, renfermant 
un pareil À est peu différente de la période de la variation de v, 
par conséquent elle diffère peu de la durée de révolution de la pla- 
néte considérée. Des pareilles expressions nous donneront les ter- 
mes à courte période. 

Par l'intégration, dans les expressions se rapportant à S, se- 
ront agrandis seulement les termes à longue période, tandis que 
ceux à courte période vont démeurer sans changement; inversé- 
ment dans les expressions de p aprés l'intégration augmenteront 
seulement les termes à courte période et ceux qui auront une pé- 
riode longue ne seront pas changés. 

M. Brendel a remarqué, que dans le cas d'une commen- 
surabilité trés approchée entre les moyens mouvements de la pla- 
nète troublée et troublante, il arrive que les termes du second ordre 
sont encore trés importants. L'idée de Gyldén répose sur ce 
principe que dans les expressions approchées, obtenues pour P 
et Q, nous ordonnons les termes non suivant leur ordre ou degré 


1) Brendel, Om amvändningen etc., pg. 6. Ludendorff, Die Jupi- 
terstórungen d. kl. Pl. von Hecuba-Typus, pg. 4. 
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mais selon leur valeur véritable, qu'ils obtiennent aprés l'intégra- 
tion. C'est pour celà que nous devons d'abord pour Thule déter- 
miner dans les développements de S, p, R, W les termes à longue 
et courte période. 

Nous commençerons par séparer avec M. Brendel dans 
S, p, W, R les parties à longue et à courte période ainsi que les 
termes ordinaires, en souscrivant sous chaque lettre un des indices 
а, К, =. Nous aurons donc ainsi 


S—S8,4-$--S, ВВ, R + B, 
1=(0+R W-—tv--Wa4-K 
K —W,--W 


où 


et qv c'est le terme séculaire. 

Dans notre cas spécial le rapport du moyen nouvement 
diurne n’ de Jupiter au méme mouvement n, de la planéte Thule 
(279) est: 


"T ==) _ 3—8, 


(2, représente une petite quantité) d'où : 


Généralement l'argument des fonctions trigonométriques qui 
entrent dans nos développements a la forme suivante: 


nw-mn(l—y)v—nB-—nyW,. 


Donc, dans les termes à courte période, le coefficient de v sera 
peu différent de 1 et dans ceux à longue période ce méme argu- 
ment différera peu de 0. Maintenant, dans les expressions préce- 
dentes pour P et Q, il faut déterminer les termes de ces genres dif- 
férents; les expressions de ces deux quantités, données aux pg.42— 
43 texte pol. nous donnent tous les termes particuliers qui cara- 
ctérisent le type planétoidal considéré. Nous devons encore trans- 
former ces expressions pour pouvoir calculer les différentes par- 
ties mentionnées plus haut, des expressions de S, R et W. 
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Comme on a remarqué précédemment, ces quantités sont 
formées avec des termes à longues et à courtes périodes, qui d'ail- 
leurs dépendent des mémes arguments que dans les expressions 
pour P et Q. 

En introduisant la classification de M. Harzer des diffé- 
rents termes particuliers et faisant des transformations nécessaires, 
nous obtenons les expressions (10 pg. 40 texte pol.) où les diffé- 
rents termes de chaque catégorie de M. Harzer sont spécifiés 
par la majuscule correspondante. 

Dans (10) pag. 46 texte pol. cinq expréssions ne contiennent 
point la quantité 2,, qui, comme on le sait, est seule liée avec le 
mouvement d'une petite planéte déterminée, ces expressions se- 
ront donc appliquables au mouvement d'une petite planéte quel- 
conque!), par contre les termes contenus dans (10) l. c. doivent 
étre calculés à nouveau dans chaque cas traité spécialement. 
C'est pour celà que Gyldén a appelé ces termes caractéri- 
stiques, comme ceux qui nous déterminent dans chaque cas 
spécial les conditions du mouvement d'une certaine petite planéte. 
Dans les termes dela catégorie (C)et (D) nous trouvons dans 
l'argument la quantité W;, qui est une fonction variant peu avec 
le temps; c'est d'ailleurs la partie à longue periode de W, qui est, 
suivant la nomenclature de M. Brendel, „1а reduction du 
temps*. Cette quantité nous rend diíficile l'intégration de nos 
équations différentielles. Dans les termes de la catégorie (4) l'in- 
tégration introduit un petit diviseur à, de l'ordre de la masse per- 
turbatrice. Dans (D), (C) nous avons à en dénominateur après 
avoir integré l'expression pour S, qui comme on le sait, contient 
la masse avec l'exposant moindre que l'unité. Ces termes devien- 
nent zéro avec la masse. Nous devons aussi remarquer que l'in- 
tégration de l'expression pour S fait augmenter les termes du genre 
(4) et (C) tandis que dans p, sont augmentés seulement les termes 
de la catégorie (B) et (D). Nous pourrons donc constater que Sa 


ГА 
sera de l'ordre de m’ et Si, Ri, Ra, Wa, Wa, K,— de l'ordre A 


et que les termes ordinaires qui entrent dans S, E, K seront de 
l'ordre de la masse perturbatrice. 


') Brendel. Theorie d. kl. Plan I, pg. 106. 
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Nous trouverons donc les expressions définitives des diffé- 
rentes parties de S, R, K, qui sont données aux pages 48 et 49 du 
texte pol. Nous devons encore substituer les séries trouvées 
dans les expressions pour R et Q. Les expressions ainsi obtenues 
sont données aux pages 50—60 du texte pol. G, désigne la somme 
des termes ordinaires; Pa, qu et g, sont des coefficients numćri- 
ques qu'on peut calculer sans aucune difficulté. 

4. Nous allons effectuer maintenant l'intégration des équa- 
tions différentielles données, (IV pg. 11, V pg. 12, VI pg. 20 texte 
pol). Plus haut nous avons déjà séparé les différents termes sin- 
guliers qui entrent dans les expressions des mouvements des peti- 
tes planètes. M. M. Brendel et Kramer ont divisé les petites 
planétes suivant leur moyen mouvement diurne en certaines clas- 
ses. Suivant cette classification la planéte Thule appartient à la 
première classe). Le mouvement des petites planètes apparte- 
nant à cette classe est caractérisé par des termes de degré 0, qui 
entrent dans R et K, qui sont grands par rapport à la masse de 
Jupiter. Dans cette classe on aura aussi dans le développement 
de R un terme important qui, comme l'a démontré M. Brendel, 
appartiendra à la catégorie D de M. Harzer c. à d. caractéri- 
stique à courte période. 

Comme Thule peut se trouver en grande proximité avec Ju- 
piter, les séries développées suivant le rapport des axes ou des 
rayons vecteurs, deviennent lentement convergentes, il faut donc 
dans les limites d'exactitude?) que nous nous sommes assignées, 
dans les développements des coordonnées et dans les équat. IV, 
V, VI, omettre la troisième puissance de la masse de Jupiter, ainsi 
que dans, les expressions du degré 0, les quantités qui seront du 
troisième ordre par rapport à la masse perturbatrice, divisée par 
une petite quantité <, parce qu'elles ont encore en facteur une 
petite quantité q ou n’ ou les carrés de ces quantités. Comme 
nous poussons nos développements jusqu'à la troisième puissance 
de la masse perturbatrice, nous devons à l'exemple de M. Kra- 
mer (dans le cas de Hécube) compléter les développements de 
M. Brendel. Nous obtenons ainsi les équations (1) pg. 63, (2), 


1) Brendel, Th. d. kl. Pl. 1 pg. 92; Kramer I, pg. 20, 144. 

2) L'exactitude que nous nous sommes assignés dans la représenta- 
tion du mouvement de Thule est celle que M. Brendel admet dans sa 
„Theor. d. kl. Plan.*. 
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(3), (4) pg. 64. (5), (6) pg. 65, (9) pg. 66. (12), (13) pg. 67 et les re- 
lations (14) —(18) pgs, 68—70. Nous avons donc trouvé, dans les 
limites d'éxactitude désirée, tous les coéfficients de P et Q. Ces 
coéfficients sont des grandeurs numériques que Гоп peut calcu- 
ler sans beaucoup de difficultés. 

Maintenant nous pouvons aborder l'intégration en se ser- 
vant de la méthode de M. Brendel, que M. M. Ludendorff, 
Buchholz et Kramer ont appliqué avec tant de succés. Nous 


+ Ф j : Р ds 
commencerons par intégrer l'équation qui nous donne do En se 


limitant aux termes du degré 0 nous aurons ——Q=-4gysintao 
et comme w=(1—y»,)v— B —y Wy, contient aussi des termes 
périodiques d'un ordre supérieur à zéro (nous savons aussi 
que le terme séculaire lui appartient) nous trouverons donc après 
l'intégration de S des termes qui seront au moins de l'ordre 1, 
parceque dans 5, il ne peut être des termes du genre А et С. 

En intégrant (1) nous trouvons: (2) S;=a+-a, cos 4t -1-S,.o0 
(Sn00 désigne les termes ordinaires), aprés quelques transforma- 
tions nous obtenons l'équation (2) pg. 73 (texte pol.). 


(3) S, — a, +(9 3- 9" Bi + 0" 2) cos 4w + Ss oo. 

Dans (3) nous pouvons calculer tous les termes, sauf а, que 
nous déterminerons ultérieurement. 

Intégrons maintenant l'équation différentielle qui détermine 
la coordonnée R. Comme nous omettons pour le moment les 
quantités du premier ordre, nous devons remarquer qu'alors (p)=0, 
parceque cette quantité contient seulement les termes à courte pé- 
riode du genre B, qui sont au moins du premier degré. L'équa- 
tion (1) pg. 63 texte pol. va nous donner !): 


ФЕ dR ls. 
(4) ая + R = 28, — P, - Qo (gp), cr = +R. 


Nous pouvons écrire aussi: 
P, =p, + p, cos 4 w. 


1) Les indices 0 indiquent que les fonctions correspondantes sont de 
l'ordre 0. 


O ruchu planetoid. 8 
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S, nous est connue par (3), il faut donc calculer Q, (z) ; 
0 


l'équation 6 pag. 34 nous donne l'équation (25) pag. 73 texte pol. 
d'où, après quelques transformations nous obtenons l'équation sui- 


vante: 
2R 
O ФЕ + R=20, Ep ON HO a + On + DUB 
+ (9), Bin + (Der? cos 4w 


où les constantes autres que a, sont données par (33) et (34) pg. 
74—75 texte pol. 


L'équation (12) pg. 48 texte pol. nous donne: 
d*R и 
(6) div + R—5,—í(1—(1-4-8)*]8, cos 4w, 
= b, — (20, + 2,2) 8, cos 4w 
où 
5, — 2a, |- py et 28, +’ =— | (0), + (9), bi + (9); 11 + 
+ (DP), (D), Bra Won”; 


Au lieu de ү, ćcrivons 


NN. 
1+ à,” 
alors après quelques transformations nous trouvons: 
Ж РТИ в" 
(7) R =. + ро — — 23, аз 1 cos 4w, 


où les coefficients sont donnés par les (29) (37) —(40) pg. 66 tex. p. 

Nous avons donc déterminé tous les termes singuliers de la 
fonction Ё; nous devons encore calculer entre les termes ordinai- 
res ceux qui seront grands. Nous pouvons écrire: 


2 

(8) 21 + В = EW, оо COS nw. 
ой Ч„оо = 2 8,00 — Воо. 

En appliquant à l'intégration de l'équation (8) la méthode de 
M. Brendel, nous obtenons: 
(9) R,—9, Sin v — gą COS 0 
où les g, et g, sont donnés par les relations (43) pg. 78 tex. pol. 

Après avoir éffectué la substitution dans (9), des valeurs trou- 
vées pour les constantes, nous arrivons à la relation (44) pg. 78 
tex. pol. qui est: 
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Er 2 8.00 EE B,oo 4 
(9) R,00 arc E Ë (1 +3.) | 
4 


où, pour n= 4, S,oo est donné par (21) pg. 71 tex. pol. 


Il nous reste encore à intégrer l'équation différentielle qui 
nous donne la relation entre W et v. En nous limitant aux termes 
de l'ordre 0 nous aurons: 


(10) d = S,—2R, — 2 R, 5, + 3R,* + 3RS, — 4R}. 
et les relations (47) pg. 79 tex. pol. En substituant ces relations 
dans l’equation (10) nous trouvons enfin: 
aw 
(11) "de = Co + Y + (а, — 2a,8, — 28, + 65,8, + 
+3 а, 8:2 — 38,3) cos 4w — (a,b, —3B, + 
+ f a, ?ر‎ — 60, 2,2) cos 8w + (1 a, f,? — B,*) cos 12w, 


1—48 et Co = ao — af, +348, — 6b,f,? -- 25,. 


c, contient m’ en facteur; šB, contient m’ divisé par la petite quan- 
tité 2,. 
Nous allons intégrer (11) en se rappelant que 


(12) = + Ws +K, 


Mais y c'est la partie séculaire, donc dans W; et a fortiori 
dans K les termes à longue période du degré 0 n'existant pas, 
nous aurons: 


où 


ат -,ак 
©) du TT dy 
En comparant (11) et (13) nous trouvons 
(14) 1=0+1, 


En ayant égard à la relation (13) pg. 49 tex. pol. et différen- 
tiant (10) pg. 49 tex. pol., nous trouverons l'équation (49) pg. 80 
tex. pol. 

En substituant dans (13) cette dernière équation, on trouve 
l'équation 50 pg. 80, d'où nous pourrons conclure les relations 
entre les coéfficients qui y entrent. Puisque Thule suivant M. 
Brendel appartient aux petites planètes critiques de la première 
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classe, il est évident que nous aurons с, =0 et y3=0. Dans notre 


cas spécial y = + !) et l'équation (13) deviendra 
(14) W=yv+K. 

Si nous connaissons déjà la constante a,, nous pouvons, еп 
se rapportant aux formules (51) et (53) pg. 81 tex. pol., calculer les 
membres singuliers de l'équation (14), en se limitant toujours aux 
termes du degré zéro. 

En éffectuant maintenant sur la fonction W les mêmes trans- 
formations que nous avons appliquées à R, nous obtiendrons sans 


peine. 
(15) W. — 262 2 Bass , (nz 1). 


n 
AA) 


Déterminons à présent les constantes d'intégration. Ces con- 
stantes se reduiront au nombre de quatre, puisque nous étudions 
pour le moment le mouvement de la petite planéte dans le plan 
méme de son orbite. Les constantes à déterminer seront: k le mo- 
dule de l'excentricité, Г — la longitude du périhélie, а — le demi- 
axe de l'orbite et n la constante du mouvement; а, joue dans nos 
formules le róle de constante surnuméraire, que nous pouvons 
pour celà traiter comme une constante arbitraire. 

La constante a, nous déterminons au moyen de b,, en sup- 
posant que dans le développement de Ё la partie constante doit 
étre effectivement du premier ordre (rein erster Ordnung). La fon- 
ction W = W;+K et K est formée au moyen des fonctions А et S, 
pour celà la partie constante de W se composera des parties con- 
stantes de R et de S. Nous aurons donc: 


(16) Co t — a — 2%. 


Mais M. Brendel?) a démontré, que dans le cas d'une pe- 
tite planéte caractéristique nous avons la relation 


d, — 20, = 0— 1. 


Si la planète considérée est critique (le cas de Thule) y ne 
peut étre zéro et dans ce cas l'équation suivante: 


a, — 206, =C,+1 


1) Brendel, 1. c. pg. 92, 100. 
2) Brendel, 1. с. kap. V, sparsim. 


http://rcin.org.pl 


SUR LE MOUVEMENT DES PETITES PLANETES. 117 


a lieu aussi— c'est pourquoi nous trouvons: 
(17) %=0; TEG 

On voit donc que dans le cas qui nous intéresse, l'équation 
53 pg. 81 texte pol. et l'équat. (17) vont nous déterminer la con- 


stante ay. 
La fonction W nous pouvons écrire: 


(18) W = (0,4-т- то) v + les termes périodiques 
où 1=1+% 


mais у, est au moins du second degré, alors += 7. 
On a trouvé la relation suivante: 


w, =(1—p)v—B—yuW, 


et 

(19) 4p = 10 —uK 

donc w=w, +y K —(1—g)v — B-pW-+yK, 
mais W=W.+1v 


et w—(1—y,) + les termes period. [в =p (1 4- 1)]. 
De ces relations nous trouvons ainsi la partie séculaire de W: 
pars sec w= (1-Е pa) v. 


Rappellons— nous le définition des quantités y et p. 


On avait les relations suivantes: . 
n' n 
=; e — = 
Ж УЕ? n, № 


(n c'est une constante d'intégration „la constante du mouvement“ 
selon M. Brendel, n, c'est le mouvement diurne vrai). Dans 
les relations précédentes introduisons à et à, (qui sont des petites 
quantités) alors nous trouvons: 


E V E 
р = 4 м = 4 


On voit que dans nos expressions, dont nous avons besoin, in- 
tervient maintenant à, seul, qui dans le cas même de commensu- 
rabilité exacte est toujours beaucoup plus grand que la masse per- 
turbatrice, tandis que la quantité 2, dans le cas d'une commensu- 
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rabilité exacte, peut devenir 0. L'idée de M. Brendel est que 
dans les expressions des perturbations apparaissent, comme petits 
diviseurs, non les à mais le à. M. Brendel a démontré aussi 
que le rapport de commensurabilité exacte peut avoir lieu seule- 
ment entre n et n’ et non entre л, et »'. Dans ce cas aussi nous 
aurons à — 0. 

Nous aurons maintenant les relations suivantes !): 


t= (1 + Co 3-1) 


à,— 0,04 -21--6 — 36, 31. 
Comme nous étudions le mouvement d'une petite planète 
critique pour laquelle: 


et 


nous aurons: 
9,=05—47p. 


D'une manière tout à fait analogue nous pouvons écrire: 


n' 3—8. 
“э Мз. a qe et 2,=02, — 467. 


Quand à, devient très petit ou presque zéro, nous aurons 
alors à faire à la libration qui, dans la méthode de Gyldén n'est 
pas difficile à étudier. La quantité n, nous donne le mouvement 
moyen diurne dans le cas de la libration; nous devons aussi con- 
naître exactement cette quantité pour pouvoir déterminer la pé- 
riode de revolution de la petite planète. 

Nous devons remarquer encore que y, et y, se distinguent 
seulement par la présence de la constante y, Cette quantité est 
de second degré et comme telle pour le moment elle n'entre pas 
dans nos calculs; += est toujours positive et, dans le cas des 
planétes caractéristiques de premiére classe (auxquelles appar- 
tient Thule), étant méme du degré zéro, elle doit étre conservée 
dans nos calculs. 

5. Passons maintenant aux applications numériques de nos 
formules. La petite planète Thule (279) fut découverte par M. Pa- 
lisa à Vienne le 25 octobre 1888; au moment de la découverte 
elle était de grandeur 13.5. Les premiers éléments calculés par 


) Brendel, A, N., 3346. 
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M. M. Lange et Halm ont montré que cette nouvelle planète 
est un astre très lointain, qui se meut à peu près à la distance de 4 
unités astronomiques du Soleil. Une orbite plus éxacte a été cal- 
culée à nouveau par M. Bidschoff!), qui a utilisé dans ses cal- 
culs toutes les observations de Thule de 25/X 1888 à 10/II 1891. 
M. Bfdschof, dans son travail, trouve que Thule se rapproche 
beaucoup de Jupiter (en moyenne trois fois pendant un siécle). 

Les éléments calculés par M. Bitdschof sont donnés pg. 86 
texte pol. 

Aprés M. Bidschof, M. Wedemeyer?) à calculé une nou- 
velle orbite, en utilisant toutes les observations et en leur appli- 
quant les perturbations, calculées par la méthode d'O ppolzer. 
Les éléments de M. Wedemeyer se trouvent pg. 87 texte pol. 
Les observations de Thule, effectuées aprés la correction de son 
orbite par M. Wedemeyer, sont trés peu nombreuses, notam- 
ment en 1906 Thule füt observée le 7/X Vienne et en 1911 le 20 
et 30 avril photographiée à Heidelberg. 

La petite planéte s'approchait de plus en plus de Jupiter, ses 
perturbations augmentaient, et sans doute l'accord entre les po- 
sitions de Thule effectives et calculées devenait de moins en 
moins satisfaisant. La proximité avec Jupiter arrivait en 1912. 
M. Viljev?) à effectué à nouveau les calculs des perturbations 
et les a appliqué aux éléments de Thule de M. Wedemeyer; il 
a donné ainsi des nouveaux éléments osculateurs (pg. 88 tex. pol.) 
et une éphémeride, déduite de ces éléments. Thule ne fut point 
retrouvée en 1913. 

Pour nos calculs nous avons utilisé les élements de M. We- 
demeyer (B. J. für 1914) et les éléments de Jupiter donné par la 
,Conn. des Temps“ pour 1914, en prenant comme valeur de la 
— Dans ce travail nous nous limitons pour 
le moment à la recherche de la lacune (dans le sens de M. Bren- 


masse m’ = 
^ $9 2 

del) pour les petites planètes du type —, c'est pour celà que nous 

donnerons ici seulement les calculs qui s'y rattachent. Dans le $ 2 


|) Sitzungsberichte d. Ak. d. Wiss. Wien, Bd. 100, pg. 937. 
+) Wedemeyer, Archiv. d. deutschen Scewarte XXXI Jhgg. 
3) Astr. Nachr. 4661. 
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nous avons donné, conformément à la théorie, que nous appli- 
quons à notre cas spécial, les expressions analytiques du dévelop- 
pement de la fonction perturbatrice. Dans ces développements 
entre la quantité z, qui s'exprime au moyen du rapport du demi- 
axe de l'orbite de Thule à celui de Jupiter. Les coéfficients nu- 
mériques de ces développements sont exprimés au moyen des quan- 
tités В et y. 

Calculons ces coéfficients. Avec les éléments de Jupiter et 
de Thule nous trouvons: 


log a = 9.9156572. 


Comme cette quantité differe peu de l'unité, nous voyons 
que les coéfficients $, т ne se trouvent point ni dans les tables de 
M. Masal') ni dans celle de Gyldén?). 

Nous devons donc les calculer spécialement en se basant sur 
leurs expressions analytiques. Nous calculons d'abord f, au mo- 
yen des formules de la moyenne arithmético-geométrique; nous 
trouvons: 

log B® = 0.1117634. 


Aprés avoir calculé £,, nous calculons les autres coéfficients 
B. pour 2? —1, 2...12. Ces coéfficients nous obtiendrons le plus 
facilement en employant la méthode de M. Ludendorff3). 

La table des valeurs auxiliaires, donnée par M. Luden- 
dorff fut augmentée et au moyen de cette table ainsi transformée, 
jai obtenu les $,. Comme contrôle de nos calculs, j'ai obtenu 
les mémes coéfficients en se servant des quadratures mécaniques. 

L'accord des valeurs des coéfficients ainsi obtenus à été com- 
plet. Les coćfficients 6%, f2, f12, &*,, ont été calculés au mo- 
yen des quadratures mécani ques, et controlés par les tables de M. 
Ludendorff. Pour les coéfficients qui nous restent encore à de- 
terminer, j'ai appliqué la relation suivante: 

Ba = emn 

') Masal H. Astronomiska iakttagelser och undersökningar änstälda pá 
Stockholms Observatorium. Fjerde Bandet. Stockholm, 1891. 

2, Gyldén H. Astr. Ges. Publ. 21. 

2) Ludendorff H. А. N. Bd. 160. 
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Les résultats obtenus furent controlés aussi par la quadrature 
mécanique. Ainsi nous avons obtenu les logarithmes de tous les 
B„* qui sont données pg. 93 texte pol. 


Maintenant en appliquant la relation: 


1:8,9. (20 —1) oid gs) 

Bec dr а 

nous avons obtenu les logarithmes de q„., qui sont donnés pg. 94 
tex. pol. Nous pouvons déjà calculer les À et B (avec des indices 
différents) qui entrent dans nos formules. Ce calcul peut étre 
effectué de deux manières différentes. On peut appliquer les re- 
lations (20) et (24) pg. 27 et 29 tex. pol. et faire le calcul de À et B 
au moyen des quantités 9, ,.,, qui sont des fonctions de q,.. et des 
Pump et Qua, , fonctions a leur tour des 9,,,,. 

Cette facon de procéder est recommandée par M. Brendel, 
elle est trés pratique, quand on doit calculer tous les coéfficients 
nécessaires pour le calcul des perturbations, mais dans notre cas 
spécial où nous nous limitons aux termes du degré zéro, c'est 
beaucoup plus pratique de calculer les A et B directement en fonc- 
tion de т en se servant des tables très pratiques de M. Masal’). 

Le résultat de notre calcul est donné pg. 95 tex. pol. 

Nous avons donc trouvé tous les coéfficients nécessaires 
dans le développement de la fonction perturbatrice, nous pouvons 
maintenant calculer pour Thule la lacune correspondante dans 
l'anneau planétoidal. 

Nous avons eu l'équation suivante, oit B, est l'inconnue: 


(28, 4- 8,58, =P! +p" 8, - p" B° 


(où 2, s'exprime au moyen de f, (67) pg. 84 texte pol.). 

M. Brendel a montré que pour trouver les limites de la la- 
cune il faut déterminer f, de cette manière, qu'en assignant une 
certaine valeur pour à, cette équation nous donne seulement deux 
différentes valeurs pour 2,, qui correspondent aux deux valeurs 
de, Ces valeurs particulières déterminent les limites de la la- 
cune. En se basant sur la théorie de M. Brendel, nous pouvons 
déterminer à, de telle manière que cette fonction devienne discon- 


1) Masal Н. Formeln und Tafeln zur Berechnung der absoluten Stó- 
rungen der Planeten. Kongl. Sv. Vet. Ak. Handl. Bd. 23. 
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tinue, alors suivant M Brendel à ces valeurs de à, qprresnon 
dent les 8,, qui déterminent la lacune. 


Puisque 2, ne peut étre zéro donc m ne peut être égale 
1 


à une fraction simple, tandis que à peut être égale à zéro, quand y. 
devient une fraction simple. 

C'est pourquoi dans toutes nos formules où nous trouvons 
en dénominateur à,, nous ne pouvons le remplacer par à, même 
dans le cas lorsque ces deux quantités diffèrent très peu l’une de 
l'autre 

Nous avons maintenant à calculer 8,. Si à, ne peut être de 
l'ordre 0, alors 2,28, est d'un ordre plus élévé que les autres ter- 
mes qui se trouvent dans notre formule. Au lieu de à, écrivons 


ò — 4p y =ò — бр", 
nous aurons: 


285,— 12p p, — p! b ds +p" pr? 

g. 2 28 
6! 4-28? JB, +k=0 3 
pu DR Loup 

— — o 12y. . 
Résolvons l'équation (1) en admettant qu'elle doit avoir une 
racine double; nous trouvons pour f, les valeurs suivantes: 
log(8,)) = 8.759360;  log(8,)u — 9.154128, 


ou bien: 


où 


Les deux valeurs de 8, ainsi obtenues sont justement les va- 
leurs critiques. Nous pouvons aussi sans difficulté calculer 2. 
Nous trouvons: 

log à — 8.559884 ; 


Connaissant (2,); et (8,)и nous obtiendrons deux valeurs pour 
т en appliquant la formule suivante, donnée par M. Brendel!): 
1 


"ure 


1) Brendel A. N. 3346. 
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Les valeurs de y sont: 
log yı = 7.696539; log yı = 8.495577 
et les valeurs correspondantes de 2,. 
102(8,), = 8.334772;  log(2,)i = 8.750883, 
Maintenant de la relation 
TUE. 
1 3—à, 
nous trouvons les deux valeurs du moyen mouvement diurne vrai, 
qui limitent la lacune pour le type $. Ces valeurs sont: 
(n,)i — 401". 732; (n,)n=391". 485 
La constante n pour les valeurs critiques à, sera: 
_ 4" 
go m 
d'où: 
n= 403". 723. 
Les résultats, que nous avons trouvés peuvent s'exprimer de 
la maniére suivante: 
Pour les petites planètes du type $ la quan- 
tité В doit satisfaire aux conditions suivantes: 


log 8,<8.759360 ou log f, 29.154128, 
log 8, <8.334772 ou log 2,=8 750883, 


d'où le vrai mouvement diurne ж pour le type Thule 
doit être 


et 


т, > 401". 732 ou n, < 391". 485. 


Entre ces deux valeurs de п, se trouve la lacune, dans le 
sens de M. Brendel, dans l'anneau planétoidal; une petite pla- 
néte avec un moyen mouvement diurne vrai (comme il est défini 
par M. Brendel) contenu entre ces limites de n, aurait une or- 
bite instable. 

(Il faut remarquer qu'on ne peut comparer n, immédiatement 
avec le moyen mouvement diurne osculateur, calculé avec les élé- 
ments de l'orbite de la petite planéte). 
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Nous avons donc obt enu les valeurs critiques du vrai moyen 
mouvement diurne qui nous déterminent la lacune dans le type 1. 
Cette lacune jusqu’à présent n’a point été déterminée. 

Cet élément ne fut déterminé par M. Brendel que dans le 
cas d'Hécube, par M. Buchholz dans le cas de Hilda et par 
M. M. Schwarzschild et Dziewulski dans le cas de Hécube 
et Hilda. Les deux derniers auteurs employèrent la méthode de 
Poincaré et recherchaient les solutions périodiques; les résul- 
tats qu'ils ont obtenus, sont en parfait accord avec ceux de M. 
Brendel et confirment ainsi la haute valeur de la nouvelle mé- 
thode. 
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DO NABYCIA WE WSZYSTKICH KSIEGARNIACH 


NASTEPUJACE DZIELA, 
WYDANE 7 ZAPOMOGI 


KASY POMOCY DLA OSÓB PRACUJACYCH NA POLU NAUKOWEM 
imienia D-RA MED. JÓZEFA MIANOWSKIEGO, 


LUB OFIAROWANE NA RZECZ KASY. 


MATEMATYKA. Rb. 


Autenrieth Ed. Mechanika Techniczna. Podrecznik nauki statyki i dy- 
namiki dla mech. inż. budow. Przeł. z upow. autora Stanisław 
Patschke 1900, XI + 613, rys. 327 . . 2 
Baraniecki M. A. Początkowy wykład syntetyczny własności przecięć 
stożkowych na podstawie ich pokrewieństwa harmonicznego z ko- 
łem 1885, XVI + 131, rys. 63 . - 
Danielewicz A. Metoda najmniejszych kwadratów 1904, X, sb 185+X 1 
Podstawy matematyczne ubezpieczeñ Zyciowych 1895, 335, tabl. X — 
Danielewicz Boleslaw. Z dziedziny statystyki matematycznej 1884, 30 — 
Dedekind Ryszard. Ciagtosé a liczby niewymierne czwatego nie- 
zmienionego wydania przełożył St. Straszewicz 1914, VIII + 17 — 
Enriques F. Zagadnienia dotyczace geometryi elementarnej. Tom I. 
Krytyka podstaw. Z drugiego wydania ое przełożyli: de 
Kwietniewski i Wł. Wojtowicz 1914, IV + 331 2 2 
Feldblum M. Geometrya wykrešlna 1902. XVI ys 395, rys. 172 ; 
Folkierski WI. Zasady rachunku rózniczkowego i calkowego, wyd. 2, 
znacznie zmien. Tom 1, 1904, XII + 574 + 55 fig. . . . . 

— Tom Il, 1909, XI + 630 + 101 fig. . 

Franke J. М. Mechanika teoretyczna 1887—1889, XXXI + 645 +72 drzew. 

Gauss K. F. Rozważania ogólne o powierzchniach krzywych, z orygi- 
nału łacińskiego pyi Stefan Finkelkraut, objaśnieniami i uwa- 
gami opatrzył Juliusz Rudnicki z przedmową Kazimierza с 
skiego, wyd. Wladystaw Wojtowicz 1913, 71. .— 
Gosiewski Wł. Rachunek prawdopodobieństwa 1906, X + 965. ` 2 
Hardy G. H. M. A. Wyklady elementarne z soy ШАШУ, przelozyl 

Wi. Wojtowicz 1916, XV + 492. . 3 
Min O sposobach obliczania przeszkód ` biegu ciat niebieskich 

1901, XII, 624 + 3 tabl. 7 
Loria Gino. Przeszłość i stan obecny najważniejszych teoryi geome- 
trycznych. Przekł. uzup. dod. S. Dicksteina 1889, VII + 112. . — 
Majewski W. Geometrya praktyczna. Podręcznik dla rzemieślników. 
1903, 301 (karton . . Е 
Malanowicz J. Kreślenie geometryczne i jego praktyczne zastosowania. 

1907, XI +- 176 + 45 tabl. i 346 rys. (karton) . . _ 

— Kreslenie i zdobienie geometryczne. Kurs pierwszy, wyd. 2, 1912, 

УШ + 110 + II + 45 tab. i 225 rys. (karton) . . . а 
— Rzuty geometryczne. (Geometrya wykreślna dostosowana do po- 

trzeb przemysłu i rzemiosł). ascia dla szkôt techn. i zawod. 

z 269 rys. 1913. . . E 
Pieri Mario. Gouna elementarna oparta | na "pojęciach „Punktu i Kuli“ 
z oryginału włoskiego przełożył Stefan Kwietniewski 1915, 152 . 1 
Poradnik dla samouków. Wskazöwki metodyczne dla studyujacych 
poszczególne nauki. Wydawnictwo A. Heflicha i St. Michalskie- 

o. Wydanie nowe. Tom l-szy w opracowaniu J. Lukasiewicza, 

. Janiszewskiego, St. by 1915, XXXIX St. Mazurkiewicza, Wt. 
Sierpiñskiego i St. Zaremby 1915, XXXIX + 618, z 314 fig. w tek- 
ście i 1 tablicą. . a 2 
Prace matematyczno - fizyczne, wydawane | przy "współudziale "Wład. 
Natansona, J. Puzyny, M. Smoluchowskiego, S. Zaremby, K. Zó- 
rawskiego przez S. Dicksteina. Tom XXVil, 1916, 277, rys. 10 . — 


NN t€ 


http://rcin.org.pl 


k. 


75 


40 


Rozmarynowicz T. Matematyczne podstawy ubezpleczeñ na wypadek 
niezdolności do pracy w zastosowaniu do urządzenia kas emeryt. 
Wyd. Bolesław Danielewicz, 1886, IV + 53.. . . . + + +. . 

Schür Е. Podręcznik geometryi analitycznej. Przelo2. T. Łopuszański 
1901,25 296. nib. Че де 

Sierpiñski WI. dr. Teorya liczb niewymiernych. Wyktady uniwers. 
1910. она 

— Teorya liczb (kurs uniwersytecki) 1914, XIV + 412 . . . 

Sochocki. Rozwiazywanie równañ fczebnych. 1884, XII + 212 

Sprawozdania z pi$miennictwa naukowego polskiego w dziedzinie nauk 
matematyczno - przyrodniczych: 

I — (1882) 1883, III + 187, 
H — (1883) 1885, 111 + 222, 
HI — (1884) 1886, 294 + 10, 
IV — (1885) 1887, 300 + VI, kazdy tom tps 

Steiner J. Konstrukcye geometryczne wykonane za pomocą linii prostej 
i stałego koła; przełożył Stefan Kwietniewski 1915, VIII + 69, 
а Wa CAS AES оо анаа ве а ДА 

Szczepański Józef. Stopień wyższy matematyki elementarnej i poczat- 
ki rachunku nieskończonostkowego. Podręcznik dla klas wyższych 
szkół średnich i dla samouków. 1906, 450 . . . . . . $ 

Wektor Czasopismo matem. fiz. Red. Władysław Wojtowicz w War- 
szawie. 5,— z przesyłką 6, — zeszyt pojedyńczy . E 2446 

Witwiñski Romuald. O ukladach odwracalnych powierzchni potrójnie 
ortogonalnych. Odbitka z tomu XXVII-go Prac mat.-fizycznych. 
JARO A bA u Z e. ON RCE Жыз 

Zaremba Stanislaw. Wstep do analizy. Czesé I-sza. Pojecie dowodu 
matematycznego oraz inne wiadomości pomocnicze 1915, 124. 


NAUKI PRZYRODNICZE. 
ii Ely dr. Flora Tatr, Pienin i Beskidu Zachodniego, 1890, VI + 
Braun Juljan. Badania w dziedzinie azotowych związków organicznych 
i ich pochodnych (1900 — 1908). 1908, VII + 238 . : RME 
Chmielewski Z. Podrecznik analizy chemiczno-rolniczej 1905, 169 
Dyakowski B. Zarys metodyki elementarnego kursu historyi naturalnej. 
Wyd. W. Jezicrski. 1909, 38. bL АИ NUM SEC. 
Dziedzicki H. Atlas organów rozrodczych (hypopygium) typów Winnertz'a 
i gatunków znajdujacych sie w jego zbiorze. Mycetophilidów, 
130, ДОЛЛАР e CPC E IS LPS he EMIR 
Dzieje myśli. Tom I zesz. 1. О rozwoju metod badań naukowych. Wie- 
dza ludów pierwotnych. Dzieje astronomii. Rys rozwoju fizyki. 
W opracow. Wł. Heinricha, Ludwika Krzywickiego. Stanisława 
Kramsztyka i Ludwika Brunera, 1907. XXXI + 296, rys. 82 
— Tom | zesz. 2. Rozwój historyczny pojęć chemicznych. Szkic 
ewolucyi pojęć w mineralogii. Zarys rozwoju matematyki: a) roz- 
wój matematyki do końca XVI w., b) zarys rozwoju geometryi 
w starożytności, wiekach średnich i w epoce odrodzenia, c) rozwój 
matematyki od początku w XVII. W oprac. Leona Marchlewskiego, 
Józefa Siomy, M. Feldbluma, Wł. Smosarskiego i Stef. Kwietniew- 
skiego, 1911, 279. z 33 ilustracyami . . . . . « . + + so 
— Tom Il zesz. 1. Historya ogólnej nauki o ziemi (geografii — geolo- 
gii). Dzieje nauk biologicznych. Dzieje antropologii. Dopełnienie 
do historyi fizyki. W oprac. Wacława Nałkowskiego, Józefa Nus- 
bauma, Ludwika Krzywickiego i L. Brunera. 1907, 471, 40 illu- 
stracyi w tekście, 2 tablice . . . . . РАЙ Ce SAA 
— Tom И zesz. 2. Dzieje psychologii. Dzieje jezykoznawstwa. 
oś и $. Lorii i J. Baudouina de Courtenay. Warszawa, 1909, 
Faraday. Dzieje świecy, przeł. M. i St. Kalinowscy. Warszawa, 1914, 
SUL IA SD 203) mika ir m x de „AEO e 
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Filipowicz Kazimierz dr. Wiadomosci p z botaniki (podiu 
dziela d-ra Le Maout: ,Leçons élémentaires de botanique“) z 19 
drzew. w tekście, 1884, Ш + 225 + Il. (karton) . . . . . . 

Grzybowski J. prof. Przeglądowa mapa geologiczna ziem polskich 
z tekstem objaśniającym z trzema przekrojami, pod redakcyą prof. 
J. Morozewicza, wydał Z. Weyberg, 1912, 139, 1 mapa kolorowa . 

Guenther Konrad. Zagadnienia życia w świetle darwinizmu. Z upoważ- 
nienia autora spolszczyli Ad. Kudelski i Kazimierz Kulwieć. 1906, 
А OES жшт PELO (m dema во PO > 

Holleman A. F. prof. Podręcznik chemii сы z 3 niem. wyd. 

rzel, według 7 wyd. niem. poprawił К. Jabłczyński wyd. 2. 1910, 
HOHEN LR SNS то с al Ie RR 
ledrzejewicz J. Kosmografia. Wyd. 2, oprac. przez d-ra M. Ernsta, 
z 246 fig. w tekście i 11 tabl. 1907, XVI- 442 . . . . . . . 

Joubert J. Zasady elektryczno$ci. Z czwartego wydania francuskiego 
przełożył Maryan Grotowski. 1915, XV - 507 z 354 rys. w tekście 

Klein P. Metereologia ogólna. Przełożył R. Merecki. Warszawa. 1916, 
92 eV tapiti; ое РЕТ КЕЯ RERS 

кашек January. Stosunki florystyczne jeziora Switezi. 1915, 
VII + 437. Sprostowania str. 7. . . . . . . . . SR 

Kontkiewicz S. Krótki podręcznik mineralogii, 1907, V + 226 + 3 tabl. 
(КОЙЫ 5 ro wie WIZY ар 2899 к! NEC) TAPETA 

Kozłowski Wł. M. Zasady przyrodoznawstwa w świetle teoryi pozna- 
EN E EAS an set S KQ N et eie lr A 

Kulwieó Kazimierz. Chrząszcze polskie. Klucz do określania owadów 
tę cars M dla użytku młodzieży, amatorów i ogrodników, 
1907, NE 4 


Loth E. Wskazówki do badañ antropol. na czlowieku 2ywym. 1914 
Malinowski Edmund dr. Świat roślin. O kształtach roślin, powstawaniu 
gov. krążeniu soków w roślinach. 1912, VI + 2 nlb 145 + 
nib., 108 rys., 2 tabl. barwne. . "NM TS БЫ с 
Męczkowska T., Rychterówna St. Zbiór ćwiczeń i doświadczeń z przy- 
rody martwej (202 doświadczeń z 112 rysunkami), 1915, 156 . `. 
wide | G. UNUM nad mieszańcami roślin, przełoż. W. Wolska. 1915, 
Merczyng H. Teorya prądu elektrycznego. Zarys zasadniczych praw 
ustalonego i nieustalonego prądu elektrycznego i towarzyszących 
mu zakłóceń magnetycznych. Podstawy elektromagnetycznej teoryi 
AA 1905 DA SO uo а жу» 
Merecki R. Klimatologia ziem polskich. Warszawa. 1914, 313. . . 
Milobedzki Tadeusz. Szkola analizy jakosciowej. 1910, VIII—271 (Karton) 
Mohn H. Zasady meteorologii, przełożył St. Kramsztyk. 1888, XVI + 
sap UN E Ла › лш! «ч» wi cha IN OY: з 
Neumayr M. prof. Dzieje ziemi, w opr. prof. d-ra Wiktora Uhliga: 

l. Geologia ogólna. Wyd. 2 pod red. J. Morozewicza, opra- 
cował K. Koziorowski, z dopełn. M. Limanowskiego, 1912, XX+-837, 
mapa barwna, 16 tablic, 300 rys. w tekście. A LĄ CV вА 

|. Geologia opisowa, przeł. z 2 niem. wyd. J. Lewiński 
i K. Koziorowski; dopełnienia poczynili: K. Bohdanowicz i J. Grzy- 
bowski. Wydał J. Morozewicz. 1908. XVI + 674 + 343 rys. 
w tekście, 2 mapy barwne, 9 tabl. (1 kolor.) . Rua zę vi D 

Nusbaum Józef dr. prof. Zasady anatomii porównawczej. 

l. Wiadomości wstępne i anatomia porównawcza zwierząt 
sis ii 1899, Ш + 744 + XXI, rys 212 + tabl. 5. 

I. Anatomia porównawcza zwierząt kręgowych. 1903, X + 
552, 134 rys. E E в UN E 

Nusbaum J. dr. Zootomia praktyczna. Wyd. staraniem d-ra Jana Tura, 
ze 100 drzeworytami. 1908, VIII + 263 . SEA CLP IS 
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Nusbaum Hilarowicz Józef. Szlakami nauki ojczystej. Wyd. Jan Tur. 
1916; ALIS LE POON, « vds узе Salal š 
Pamiętnik Fizyograficzny, wydany staraniem E. Dziewulskiego i B. Zna- 
towicza: 
Tom Ш. Dział I Meteorologia i hydrografia. 11. Geologia 
z chemią. Ш. Botanika i zoologia. IV. Antropologia. V. Miscelańea. 
1883 536 + 2 + 2 + 13 tab., rys lit., 21 drzeworytow w tekście; 
V. Dział I, II, III, IV, V. 1885, 4 nib. 113 + 76 + 223 + 
744111 4- 4. 
VIII. Dział I, П, Ш. IV, V. 1888, 2 nlb. + XIX + 19 + 
155 + 389 + 17 + 33 +4 nlb.; 27 tabl rys. lit. i drzew. w tekście; 
WINE A. Slósarski i Br. Znatowicz. 
IX. Dział I. Il, Ш, IV. 1889 2 nlb.4+XIX+235445+11+ 
295+ 77-Н\У, 24 tabl. гуз. litogr. i drzewor. 
X. Dział I, H, Ш, IV. 1890. 2 nib. + XXI + 202 + 75 + 
437 + 2 nlb. + 20 +11 -+ ll. 29 tabl. rys. lit. i drzew. w tekście. 
ХІ. Dział I. Il, Ш. 1891, 8 + 18 + 186 + 162 + 133 + II + 
H 14 tabl. гуз. litogr. i drzewor. w tekscie. 
XI. Dział I, ЇЇ, Ш, IV. 1892. 17 + 214 + 235 + 23 + II + 
II + 12 tabl. rys. litogr. i drzewor. w tekscie. 
ХШ. Dział 1, 11, Ш. 1895, 19 +- 152 +- 231 + 1+ 1 +7 ta- 
blic rys. lit. 
XIV. Dział I, Il, Ш. 1896, 23 + 151 + 30 + 228 +1 + I + 
7 tabl. rys. lit. 
Wydawcy: W. Wróblewski i Br Znatowicz. 
XV. Dział I, Il, Ш. 1898, 19 + 183 + 285 + 39 +1 + I + 
4 mapy + 3 tabl. lit. 
XVI. Dział І, II, Ш. 1900, 13+139+ 13 + 44-+-208. 
ХУП. Dział I, II, Ш, IV. 1902, 16 + 134 + 144 + 104 + 
22--]--1--1 mapa i tabl. litogr. 
XVIII. Dział 1, II, Ш, IV, V. 1904, 61 + 193 + 147 + 
104 + 24 4- 42 + + I. 
XIX. Dział I, П, Ш, IV. 1907, 79 + 183 + 59 -- 82 +- 7 -- I+ I. 
Pamietnik Fizyograficzny, wyd star. E. voi TS à i B. Znatowicza. 
XX. Meteorologia i Miscelanea 1910, XLI + 203 + 46, 
Wydawcy: K. Kulwieé i K. Stolyhwo. 
ХХІ. Dział I, II, Ш, IV, V. 1913, IX + XV + 155 + 30 + 
25 + 117 + 48 + 41 + 4 mapy + 19 rys. + 24 tabl. fot. 
XXII. 1, Il, Ш, IV, V. 1914 IX + XV + 155 + 30 + 25 + 
117 + 48 + 41 + 4 тару + 19 rys. + 24 tabl. fot. 
XXII. Dział I, П, HI, IV, V. 1916, 262 + XXVIII tabl. + 
142-1- 1 tabl. + 246 + 1 mapa + 16 rys. + 16 + 8, każdy tom 
Pogorzelski W. Badania teoretyczne ilo$ci ciepla, otrzymywanych na 
a ziemskiej, z uwzglednieniem strat promieniowania w atmos- 
Ee Bes w. ee Se E Nera Жы MEER S PA BEAC T SN EY 
Pol G. Słownik łacińsko - polski nazw gatunkowych roślin (12 + 17), 
Pawłowski Stanisław. Ze studyów nad zlodowaceniem Czarnohory. 
IO MOI ZARABIA 2 Met ce e ła vi^ Е 
Pożaryski Mieczysław. Podstawy naukowe elektrotechniki łącznie z za- 
sadami pomiarów. 1915 X +-415, z 427 rysunkami w tekście . 
Roszkowski W. i Żebrowska A. O budowie pochewek prącia u błotnia- 
ER: PIO IO UBI 4012 PRM EU EE Aoc c m 
Routh E. |. Statyka teoretyczna z licznymi przykladami z drugiego wy- 
dania Sup cH go przełożył Zygmunt Straszewicz. 1916, X + 
TS TRIE de za rc ea ЫЕ «жы 
Rydzewski Bronistaw. Próba charakterystyki paleobotanicznej Dabrow- 
skiego Zagłębia Węglowego. 1915, 86, tabl. 5 . . . . . . 
Siemiradzki 1. Gąbczaki jurajskie ziem polskich (Paleontologia ziem 
polskich pod red. J. Lewiñskiego № 1), 1913, 49 + tabl. VIII . 
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Silberstein Ludwik. ia irr zm. |. 1908, NN 366. . 350 


Ш. cz 1, 1913, 173 . . 1 80 
Slownik ficzny Król, Pol. i innych krajów słowiańskich. Komplet . 60 — 
нна iatopelk Zygmunt, Atlas ziem polskich Тот 1, część } ME 
kie Ksi ęstwo Poznańskie, 46 map i planów. 1915 . ‚ei = 
. Stra E. dr, Jost L. dr., Schenk K. dr., Karsten d. dr Podre 
aniki día szkół wyższych. Z E wyd niem. przelo2yli Ja 
i Karol Szteinbokowie. Zeszyt 1. 1913, 160. Zeszyt И. 19 4, 
161 — 320. Zeszyt Ш. 1915... 3 — 
Świat i człowiek. Zeszyt I. wyd. 2. Pojęcie rozwoju. Wszechświat 
еро rozwój. Rozwój ziemi, оргас, 1. Waserberg, S. Kapusi; 
alkowski, 1908, XVI + 215 -+82 ilustr. -+ 3 tablice kolor. . 135 
Zeszyt ll, wyd. 2. Rozwój życia оте N Parry, 
roślin. pargar zwierząt. Pochodzenie człowieka. Roz 
wieka, орг usbaum, Z. Wóycicki, J. Eismont, K. мо 
KO ki, 1912, 821 + 73 ilustr. +1 tabl. Ta 1 60 
ts na md © W PLK кем e е А 
w sun o arczyc 0 rzyw i ppe- 
ANSE 191 sir 350% 65 hr. 'R 'ój hic Jw 
е52у wy 'OZWÓ, spoleczny. 02W psyc zny. 
Rozwój w dziejach sztuki. Znaczenie rozwoju. W opr. L. Krzy- 


wickiego, M. wskiego, Wł. Tatarkiewicza i F. ARCH 
Rios e 1913, str. 355-+5 ilustr. . . 2= 
kac: W. W. Aa Poczatek i rozwój umysłowości w przyrodzie. ` 1885, id 
Tenenbaum Szymon: Fauna koleopterologiczna wysp Balearskich. "1915. 2 

Tombeck D. i Gouard E. Chemia przemyslowa, przełożył 7. Haraba- 
M,  szewski. 1915, XI + 422. . 1 80 

Tür Jan. Nowe badania nad rożwojem układu nerwowego potworów 
platyneurycznych. 1915, 128. — 50 

w E. Zbiorowiska roślinne. Zarys ekologicznej grafii roślin. 


wydania niem. E. Knoblaucha przeł. z "OF: autora E. Strumpf 

i Trzebiński. 1900, XV + 450 1 50 
Witkawaki EN prof. uniwers. Jagiellonskiego. Zasady fizyki. Tom L 
d. 3. (Fizyka ogólna. Dynamiczne w $ci quer Ayka, 

V -+ 536 + 205 fig. . 2 — 

| oo 1, wydanie 4-te, 1915. XX 4 535. ry s. 205 . 2 4 
AE И, wyd. 2. (Cieplo. Fizyka costes owe. Promienio- 

wanie.). + GSC ig + 2 tabl. kol. 2 40 
rom ill (сикы TORBAY ZAM) 1914, IX + 1 nbl.+ 

Mox Bal) ba k ki do hydrograf. dí we. 
2 ї ога, өз st 0 aśniając ma r awne 

Ve qu eiie pótnocno - йе Чо à 2 у 14195 +1 nbl. — 5 
куч Zygmunt. Obrazy roślinności Królestwa Polekieżo. Zeszyt 1. 
Roślinność ШАУ. SE гин 1911, 12 nlb + tabl. 10 + 20 

str. nlb. objasnie (cb 
lin IL Roślinność wyżyny Kielecko - Sandomierskiej. 

1912, 36 +- 10 tabl. 1 — 
Zeszyt III. „Roślinność: „wyżyny Kielecko - Sandomierskiej 

1912, 32 + Yo tabl. . 
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Zeszyt VI. Roślinność Ojcowa. 1913, 26 + 10 КЫ 
су! VIL Roślinność okolic Częstochowy i Olsztyna. 1914. 
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„Ziemia” Tygodnik Ilustrowany Redaktor i wyd. K Kulwieé. Rocznie . 6 — 
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